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detonation  parameters  such  as  the  Chapman-Jouguet  velocity, 
pressure,  etc.  which  are  derived  from  one-dimensional  steady 
state  equilibrium  theory  (hence  static).  Experience  shows 
that  the  dynamic  parameters  reflect  more  intimately  the 
detonation  properties  which  are  of  great  practical  interest 
(e.g.  initiation,  transmission,  detonability , etc.-) .  In 
terms  of  critical  tube  diameter,'  the  universality  of  the 
important  empirical  law  of  dc  ■  13X  has  been  conclusively 
verified  in  the  current  work.  Moreover,  the  study  of  the 
critical  tube  diameter  which  is  essentially  the  problem  of. 
the  transmission  of  detonation  wave  from  a  confined  circular 
tube  to  the  unconfined  three-dimensional  geometry  has  .been 
extended  to  other  transmission  geometries  and  the  results 
correlated  succesfully  again  in  terms  of  a  simple  empirical 
model.  The  validity  of  the  surface  energy  model  to  correlate 
the  detonation  cell  size  and  the  critical  initiation  energy 
has  been  further  confirmed  for  a  broad  range  of  explosive 
mixtures.  This  simple  phenomenological  model  has  been  used 
to  succesfully  link,  on  a  quantative  basis,  the  detonation 
cell  size  X  with  the  critical  initiation  energy  not  only 
for  hydrogen-air  mixtures  but  for  most  of  the  common  gaseous 
hydrocarbon  fuels  ranging  from  ethylene  to  the  alkanes. 

Finally,  detonability  limits  have  also  been  correlated  in 
terms  of  thecell  size  X.  Apart  from  the  direct  importance 
of  being  able  to  deduce  detonability  limits  again  from  a 
relatively  straight-forward  procedure  of  measuring  the  cell 
size,  there  arises  the  further  important  conclusion  that  one 
cannot  define  detonability  limits  in  terms  of  thermochemical 
properties  of  the  mixture  alone  but  that  one  must  also  define 
the  confining  geometrical  boundary  conditions  to  unambiguously 
define  the  true  detonability  limits  for  a  given  mixture  in  a 
configuration. 

In  the  current  program  year,  emphasis  was  also  placed 
on  the  study  of  the  problem  of  intense  turbulent  mixing 
between  the  burnt  and  the  unburnt  explosive  mixture. to  explore 
specifically  the  environment  under  which  rapid  amplification 
to  detonation  can  be  achieved.  The  turbulent  mixing  was  '! 
induced  by  having  a  weakly  initiated  flame  propagate  into 
an  obstacle  environment  along  a  linear  combustion  tube  and 
thereby  create  the  intense  turbulent  flow  field  ahead  of 
itself  by  virtue  of  fluid  displacement  as  a  result  of  the. 
increase  of  £he  specific  volume . in  an  exothermic  reaction. 

The  conditions  of  turbulent  mixing  between  the  burnt  and  unburnt 
FAE  mixture,  the  resulting  rapid  flame  acceleration  and  the 
minimum  necessary  conditions  for  transition  to  detonation 
have  been  studied  in  the  present  program. 
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Introduction 


1 . 

The  underlying  task  in  FAE  research  as  it  pertains  to 
AFOSR  needs  can  be  summarized  in  the  following  way:  For  a  given 
fuel-air  mixture  at  prescribed  initial  and  boundary  conditions, 
the  problem  is  to  determine  the  essential  requirements  for 
direct  initiation  and  sustained  propagation  of  a  detonation  wave 
within  the  FAE  cloud.  This  implies  establishing  a  priori  on  a 
quantitative  basis  the  minimum  required  spatial  extent  of  the 
activation  region  for  initiation,  the  characteristic  chemical 
time  scale  and  the  energetics  required  for  the  subsequent  evo¬ 
lution  of  the  detonation  wave  in  the  cloud.  Our  approach  in 
the  AFOSR  sponsored  program  at  McGill  has  been  to  look  at  the 
fundamental  physical  and  thermo-chemical  phenomena  underlying 
these  events  from  a  sufficiently  broad  range  of  perspectives 
such  that  the  important  controlling  mechanisms  can  be  identified 
and  thereby  the  events  can  be  anticipated  and  controlled  on 
a  quantitative  parametric  basis.  This  approach  has  entailed 
a  broad  range  of  in-house  exploratory  laboratory  scale  experiments 
followed  up  by  theoretical  modelling  on  a  phenomenological  basis 
to  lay  the  groundwork  and  establish  quidelines  for  large  scale 
field  experiments  done  in  collaboration  with  government  and  quasi¬ 
government  research  agencies  which  have  large  scale  field  test 
facilities.  The  present  progress  report  describes  the  results 
that  have  evolved  from  such  an  approach  where  the  activity  perti¬ 
nent  to  AFOSR  needs  has  gone  beyond  the  immediate  in-house  research 
work  sponsored  at  McGill  and  benefited  in  scope,  scale,  level  of 
activity  and  productivity  by  having  access  to  and  interaction 
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with  external  agencies  in  the  manner  described.  Specifically, 
in  the  current  program  year,  the  results  of  which  are  being 
reported  herein,  large  scale  tests  have  been  conducted  in 
collaboration  with  Sandia  National  Laboratories,  Albuquerque, 

New  Mexico  and  the  Defense  Research  Establishment  Suffield, 

Alberta  and  the  Norwegian  Defence  Construction  Services,  Oslo, 
Norway.  The  results  of  these  tests  form  an  integral  part  of 
the  progress  pertinent  to  the  AFOSR  program  and  are  included 
herein  as  discrete,  self-contained  but  intimately  related 
sections  of  the  progress  report.  The  whole  report  is  composed 
of  individual  Appendices  describing  the  progress  achieved  in 
the  various  aspects  of  the  research.  All  of  the  Appendices 
form  part  of  a  series  of  publications  which  have  appeared  recently 
in  archive  journals  in  the  field  or  are  being  presented  at  inter¬ 
national  conferences  as  indicated  on  the  introductory  page  of 
each  study.  In  the  next  several  paragraphs  a  summary  description 
of  the  progress  report  is  presented  to  assist  in  placing  the 
progress  reported  in  each  Appendix  in  the  context  of  the  AFOSR 
program. 


2 .  Measurement  of  Detonation  Cell  Size  (Chemical  Length  Scale) 
and  Correlation  with  the  Dynamic  Detonation  Parameters 

Over  the  span  of  the  past  several  years,  extensive 
studies  largely  carried  out  in  our  laboratory  have  shown  that 
the  detonation  cell  size  X  is  of  great  fundamental  significance 
to  describe  an  overall  chemical  time  scale  that  truly  characterizes 


the  chemical  reactions  in  the  transient  three-dimensional 
reaction  zone  of  the  complex  multi-ieaded  (cellular)  structure 
of  a  real  detonation  wave  itself.  It  is  concluded  that  the 
detonation  cell  size  X  (or  the  equivalent  chemical  times  de¬ 
rived  from  it)  provide  a  much  more  useful  parameter  as  input 
to  the  required  theories  and  empirical  laws  for  the  prediction 
of  the  desired  requirements  for  initiation  and  propagation 
than  the  currently  used  induction  time  t  itself.  It  is  the 
first  and  perhaps  the  most  basic  in  the  hierarchy  of  so-called 
"dynamic"  detonation  parameters  which  are  parameters  that  depend 
on  the  rates  of  chemical  reactions  (hence  "dynamic")  as  opposed 
to  other  detonation  parameters  such  as  the  Chapma n- Jouguet 
velocity,  pressure,  etc.  which  are  derived  from  one-dimensional 
steady  state  equilibrium  theory  (hence  static).  Experience 
shows  that  the  dynamic  parameters  reflect  more  intimately  the 
detonation  properties  which  are  of  great  practical  interest 
(eg.  initiation,  transmission,  detonability,  etc.).  In  last 
year’s  progress  report,  a  limited  range  of  cell  size  measurements 
were  reported  and  correlated  with  the  so-called  critical  tube 
diameter  and  initiation  energy  largely  in  hydrogen-air  mixtures. 
The  critical  tube  diameter  refers  to  the  minimum  tube  diameter 
for  which  a  planar  detonation  in  a  particular  explosive  mixture 
will  transmit  without  failure  from  a  confining  tube  into  an  uncon 
fined  environment.  It  has  been  shown  to  be  a  very  important 
dynamic  detonation  parameter  because  it  can  potentially  be  used 
in  a  number  of  important  practical  estimates  in  the  context  of 
the  FAE  problem.  In  the  first  instance,  it  can  be  used  to  infer 


the  minimum  size  of  a  detonable  cloud  and  hence  the  detonability 
limits  for  a  given  set  of  boundary  conditions.  In  the  context 
of  initiation  via  chemical  sensitizers  (FAE  III  concept)  it 
can  serve  as  a  quantitative  measure  of  the  minimum  spatial  length 
scale  of  the  activation  region  required  for  detonation  initiation. 
Finally,  the  critical  tube  diameter  has  also  been  proposed  as  an 
alternate,  more  easily  measurable  parameter  to  assess  shock 
sensitivity  of  a  given  explosive  mixture  ins  ead  of  initiation 
energy  which  from  the  practical  point  of  view  is  prohibitively 
difficult  to  measure  for  other  than  relatively  sensitive  mixtures. 
Nevertheless,  initiation  energy  in  itself  still  remains  an  impor¬ 
tant  and  basic  dynamic  detonation  parameter;  its  limitation 
rests  with  the  serious  difficulty  of  its  quantitative  experimental 
measurement. 

The  current  effort  that  is  being  reported  extends 
the  preliminary  success  over  a  limited  range  of  conditions  described 
in  the  last  progress  report  of  quantitatively  correlating  the 
cell  size  X  with  the  critical  tube  diameter,  the  initiation 
energy  and  the  detonability  limits.  In  terms  of  the  critical 
tube  diameter,  the  universality  of  the  important  empirical  law 
of  dc  =  13X  has  been  conclusively  verified  in  the  current  work. 
Moreover,  the  study  of  the  critical  tube  diameter  which  is  essen¬ 
tially  the  problem  of  the  transmission  of  a  detonation  wave  from 
a  confined  circular  tube  to  the  unconfined  three-dimensional 
geometry  has  been  extended  to  other  transmission  geometries  and 
the  results  correlated  successfully  again  in  terms  of  a  simple 
empirical  model.  The  validity  of  the  surface  energy  model  to 
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correlate  the  detonation  cell  size  and  the  critical  initiation 
energy  has  been  further  confirmed  for  a  broad  range  of  explo¬ 
sive  mixtures.  This  simple  phenomenological  model  has  been 
used  to  successfully  link,on  a  quantitative  basis,the  detonation 
cell  size  X  with  the  critical  initiation  energy  not  only  for 
hydrogen-air  mixtures  but  for  most  of  the  common  gaseous  hydro¬ 
carbon  fuels  ranging  from  ethylene  to  the  alkanes.  Further 
correlation  requires  additional  data  on  the  quantitative  measure¬ 
ment  of  initiation  energy,  the  bulk  of  which  so  far  has  been 
produced  by  Ellsworth  at  the  Shell  Thornton  Research  Center  in 
Great  Britain.  The  significant  importance  of  the  cell  size  - 
initiation  energy  correlation  cannot  be  overemphasized  since 
now  via  a  relatively  simple  laboratory  scale  measurement  of 
cell  size  X,  one  can  infer  the  initiation  energy  and  hence  shock 
sensitivity  of  just  about  any  explosive  gaseous  medium.  Finally, 
detonability  limits  have  also  been  correlated  in  terms  of  the 
cell  size  X.  Apart  from  the  direct  importance  of  being  able 
to  deduce  detonability  limits  again  from  a  relatively  straight¬ 
forward  procedure  of  measuring  the  cell  size,  there  arises  the 
further  important  conclusion  that  one  cannot  define  detonability 
limits  in  terms  of  thermochemical  properties  of  the  mixture 
alone  but  that  one  must  also  define  the  confining  geometrical 
boundary  conditions  to  unambiguously  define  the  true  detonability 
limits  for  a  given  mixture  in  a  particular  configuration. 

Thus  in  terms  of  the  progress  report,  Appendix  I  describes 
in  great  detail  the  recent  successful  universal  quantitative  cor¬ 
relation  between  the  cell  size  X  and  the  other  practically 
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important  "dynamic"  detonation  parameters  such  as  the  critical 
tube  diameter,  initiation  energy  and  detonability  limits. 

Appendix  II  followed-up  by  large  scale  studies  described  in 
Appendix  III  and  Appendix  IV  pertain  to  the  assessment  of  the 
transmission  properties  of  detonation  waves  from  confined 
to  unconfined  geometries.  Appendix  V  constitutes  a  summary 
review  or  compendium  of  essentially  the  current  state  of  the 
art  in  FAE  research  derived  largely  from  the  studies  described 
herein. 

3 .  Intense  Turbulent  Mixing,  Flame  Acceleration  and  Transition 
to  Detonation 

In  the  FAE  III  concept,  one  relies  on  the  tailored 
intense  turbulent  mixing  between  the  chemical  sensitizer  and 
the  explosive  fuel-air  cloud  to  set  up  suitable  conditions  for 
’shockless"  initiation  of  detonation.  By  "shockless"  initiation 
is  meant  that  one  does  not  supply  a  priori  a  detonation  level 
shock  wave  in  the  initiation  process  as  one  does,  for  example, 
in  the  blast  initiation.  In  fact,  the  medium  initially  is 
entirely  shock  free.  Rather  the  shock  wave,  which  ultimately 
must  reach  detonation  strength,  must  evolve  spontaneously  within 
the  chemically  sensitized  activation  region  purely  by  virtue 
of  the  intensity  of  the  exothermic  chemical  reactions  relying 
exclusively  on  the  energy  content  of  the  sensitized  portion  of 
the  fuel-air  cloud  itself.  Clearly,  phenomenologically  this 
is  not  a  simple  process  and  to  work  successfully  requires  a 
very  judicious  tailoring  of  a  number  of  time  scales,  length  scales. 
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reaction  rates  and  energetics.  Suitable  conditions  must  be 
generated  such  that  once  reactions  are  intiated  rapid  ampli¬ 
fication  of  pressure  waves  that  are  produced  in  the  exothermic 
chemical  reactions  must  progress  in  a  coherent  fashion  to 
lead  to  detonation  level  shock  waves  within  the  chemically 
sensitized  region  which  must  be  of  sufficient  size,  the  scale 
of  which  has  been  alluded  to  in  the  discussion  on  the  critical 
tube  diameter  in  the  previous  section.  The  "detonation"  wave 
thus  formed  within  the  chemically  sensitized  portion  of  the 
fuel-air  cloud  will  then  transmit  into  the  rest  of  the  cloud 
leading  to  global  detonation  of  the  cloud.  The  most  difficult 
aspect  of  the  FAE  III  concept  is  to  know  a  priori  and  be  able 
to  create  experimentally  the  kind  of  conditions  whereby  the 
chemical  reactions  in  the  sensitized  portion  of  the  cloud 
progresses  in  step  (i.e.  coherently)  with  the  evolving  pressure 
field  and  thereby  provide  the  adequate  shock  amplification 
environment  which  can  ultimately  lead  to  detonating  the  ini  tally 
shock  free  medium. 

In  the  current  program  year,  emphasis  was  placed  on 
the  study  of  the  problem  of  intense  turbulent  mixing  between 
the  burnt  and  the  unburnt  explosive  mixture  to  explore  specifically 
the  environment  under  which  rapid  amplification  to  detonation 
can  be  achieved.  The  turbulent  mixing  was  induced  by  having  a 
weakly  initiated  flame  propagate  into  an  obstacle  environment 
along  a  linear  combustion  tube  and  thereby  create  the  intense 
turbulent  flow  field  ahead  of  itself  by  virtue  of  the  fluid 


displacement  as  a  result  of  the  increase  of  the  specific  volume 
in  an  exothermic  reaction.  The  conditions  of  intense  turbulent 
mixing  between  the  burnt  and  unburnt  FAE  mixture,  the  resulting 
rapid  flame  acceleration  and  the  minimum  necessary  conditions 
for  transition  to  detonation  have  been  studied  in  the  present 
program  and  are  reported  in  Appendix  VI. 
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for  initiation  and  propagation  than  the  currently  used 
induction  time  xitself.  It  is  the  first  and  perhaps  the 
most  basic  in  the  hierarchy  of  so-called  "dynamic"  detonation 
parameters  which  are  parameters  that  depend  on  the  rates 
of  chemical  reactions  (hence  "dynamic’’'*  as  opposed  to  other 
detonation  parameters  such  as  the  Chapman-Jouguet  velocity, 
pressure,  etc.  which  are  derived  from  one-dimensional  steady 
state  equilibrium  theory  (hence  static).  Experience  shows 
that  the  dynamic  parameters  reflect  more  intimately  the 
detonation  properties  which  are  of  great  practical  interest 
(e.g.  initiation,  transmission,  detonability , etc . ) .  In 
terms  of  critical  tube  diameter,'  the  universality  of  the 
important  empirical  law  of  dc  =  13X  has  been  conclusively 
verified  in  the  current  work.  Moreover,  the  study  of  the 
critical  tube  diameter  which  is  essentially  the  problem  of 
the  transmission  of  detonation  wave-  from  a  confined  circular 
tube  to  the  unconfined  three-dimensional  geometry  has  been 
extended  to  other  transmission  geometries  and  the  results 
correlated  succesfully  again  in  terms  of  a  simple  empirical 
model.  The  validity  of  the  surface  energy  model  to  correlate 
the  detonation  cell  size  and  the  critical  initiation  energy 
has  been  further  confirmed  for  a  broad  range  of  explosive 
mixtures.  This  simple  phenomenological  model  has  been  used 
to  succesfully  link,  on  a  quantative  basis,  the  detonation 
cell  size  X  with  the  critical  initiation  energy  not  only 
for  hydrogen-air  mixtures  but  for  most  of  the  common  gaseous 
hydrocarbon  fuels  ranging  from  ethylene  to  the  alkanes. 

Finally,  detonability  limits  have  also  been  correlated  in 
terms  of  thecell  size  X.  Apart  from  the  direct  importance 
of  being  able  to  deduce  detonability  limits  again  from  a, 
relatively  straight-forward  procedure  of  measuring  the  cell 
size,  there  arises  the  further  important  conclusion  that  one 
cannot  define  detonability  limits  in  terms  of  thermochemical 
properties  of  the  mixture  alone  but  that  one  must  also  define 
the  confining  geometrical  boundary  conditions  to  unambiguously  | 
define  the  true  detonability  limits  for  a  given  mixture  in  a 
configuration. 

In  the  current  program  year,  emphasis  was  also  placed 
on  the  study  of  the  problem  of  intense  turbulent  mixing 
between  the  burnt  and  the  unburnt  explosivemixture.  to  explore 
specifically  the  environment  under  which  rapid  amplification 
to  detonation  can  be  achieved.  The  turbulent  mixing  was 
induced  by  having  a  weakly  initiated  flame  propagate  into 
an  obstacle  environment  along  a  linear  combustion  tube  and 
thereby  create  the  intense  turbulent  flow  field  ahead  of 
itself  by  virtue  of  fluid  displacement  as  a  result  of  the 
increase  of  the  specific  volume  in  an  exothermic  reaction. 

The  conditions  of  turbulent  mixing  between  the  burnt  and  unburnt 
FAE  mixture,  the  resulting  rapid  flame  accelaration  and  the 
minimum  necessary  conditions  for  transition  to  detonation 
have  been  studied  in  the  present  program. 
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ABSTRACT 


Experimental  measurements  of  the  detonation  cell  size  In  mixtures  of  H~ 

^2h2*  ^4,  C'2h6*  *-3h8  and  ^4^10  a"*r  over  a  range  of  fuel  concentrations 
have  been  carried  out  In  three  cylindrical  tubes  of  diameters  5  cm,  15  cm 
and  30cm.  The  cell  size  has  been  determined  from  the  signatures  on  smoked 
aluminum  foils  placed  inside  the  tube  as  well  as  from  the  frequency  of  the 
pressure  fluctuations  recorded  by  piezoelectric  transducers.  Based  on  the 
cell  size  data  obtained,  estimates  of  the  critical  tube  diameter  using  the 
empirical  law  of  Soloukhin  and  Mitrofanov  (dc  *  13X)  have  been  found  to  be 
in  agreement  with  experimental  data  from  direct  measurement  of  the  critical 
tube  itself.  Hence,  the  Important  empirical  law  dc  »  13X  Is  thus  verified. 
Estimates  of  the  critical  charge  weight  from  the  cell  size  data  using  the 
surface  energy  theory  proposed  by  Lee  have  been  found  to  agree  reasonably 
well  with  the  experimental  results  of  Elsworth.  Based  on  the  criteria  for 
stable  prooagation  In  tubes  (d*  *  a/it)  and  in  two-dimensional  channels  (W*  *3 
A),  detonabil i ty  limits  can  also  be  predicted  from  a  knowledge  of  the  cell 
size  x.  Based  on  Westbrook's  kinetic  calculations.  It  Is  found  that  the  cell 
size  data  are  directly  proportional  to  the  Induction  time  of  the  oxidation 
process  which  confirms  qualitatively  Schelkhin's  model.  However,  on  a  quan¬ 
titative  basis,  Schelkhin's  model  predicts  cell  sizes  an  order  of  magnitude 
larger  than  the  present  experimental  datd. 


Presented  at  the  9th  International  Colloquium  on  Dynamics  of 
Explosions  and  Reactive  Systems ,  .Pol  tiers ,  France,  July  1983. 
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I .  Introduction 

An  assessment  of  the  relative  detonation  sensitivity  of  various  fuel - 
uxygen-ni trogen  mixtures  was  carried  out  by  Matsui  and  Lee1  using  the  criti¬ 
cal  initiation  energy  for  spherical  detonations  as  a  basis  for  comparison. 

The  critical  initiation  energy  itself  depends  on  the  energy-time  character¬ 
istics  as  well  as  the  geometry  of  the  source2’3  and  is  thus  a  rather  complex 
parameter.  Furthermore ,  for  the  relatively  insensitive  fuel-air  mixtures, 
the  critical  initiation  energy  is  usually  measured  in  terms  of  an  equivalent 
weight  of  a  high  explosive  charge  typically  of  the  order  of  100  g  or  more. 

This  requires  rather  large  volumes  of  gases  for  a  long  enough  detonation  tra¬ 
vel  to  conclusively  determine  the  stability  of  the  detonation  in  the  case  of 
successful  initiation.  The  use  of  long  rectangular  bags  rather  than  spheri¬ 
cal  balloons  simplifies  the  experiments  somewhat  but  suffers  from  a  certain 
degree  of  lateral  confinement  of  the  plastic  walls  as  shown  by  Moen  et  al." 
Thus,  a  rectangular  bag  does  not  simulate  a  truly  unconfined  spherical  deto¬ 
nation.  For  insensitive  fuels  or  for  near-limit  mixtures,  the  critical  cross- 
sectional  dimension  of  the  rectangular  bag  required  for  stable  propagation  is 
also  very  large  and  is  not  known  a  priori. 

In  recent  years,  the  critical  tube  diameter  required  for  the  successful 
transformation  of  a  confined  planar  detonation  into  an  unconfined  spherical 
wave  has  been  shown  to  be  easily  measurable.  Since  the  re-initiation  pro¬ 
cess  for  successful  transformation  of  a  planar  into  a  spherical  wave  occurs 
within  a  distance  of  the  order  of  one  tube  durti-cer  from  the  exit  of  the  tube, 
even  the  large-scale  exoeriments  required  are  of  rather  modest  dimensions.  Al¬ 
though  the  critical  tube  diameter  has  been  related  to  the  critical  energy  by 
Lee  and  Matsui  using  a  "work-done"  concept,  a  far  more  important  fundamental 
relationship  between  the  critical  tube  diameter  "dc"  and  the  cell  size  or 
transverse  wave  spacing  "A"  has  been  established  in  recent  years.  First  ob¬ 
served  by  Mitrofanov  and  Soloukhin6  and  confirmed  later  by  Edwards  et  al.,7 
it  was  found  that  dc  -  13A  for  C2H2-O2  mixtures  at  sub-atmospheric  pressures 
of  about  80  torr.  Following  the  suggestion  of  Edwards  that  this  empirical 
law  should  have  a  wider  applicability,  Knystautas  et  al.8  carried  out  exten¬ 
sive  experiments  in  fuel -oxygen-ni trogen  mixtures  to  measure  both  the  criti¬ 
cal  tube  diameter  "dc"  as  well  as  the  cell  size  "A"  over  a  wide  range  of  ini¬ 
tial  pressures  and  confirmed  the  validity  of  this  empirical  law  of  dc  =  13A 
for  the  fuels  tested  (H2,  C2H2,  C2H4,  C3H5,  CjHg,  CH^,  C2H5,  C4H10  and  MAPP). 
With  the  critical  tube  diameter  linked  to  the  cell  size  through  this  simple 
empirical  law, it  is  then  possible  to  perform  laboratory-scale  experiments  in 
confined  tubes  to  determine  the  cell  size  "A"  and  deduce  the  critical  tube 
diameter  from  it.  A  program  for  the  systematic  measurement  of  the  detonation 
cell  size  in  atmospheric  fuel-air  mixtures  has  been  carried  out.  The  present 
paper  reports  the  results  obtained  ar,d  the  correlation  of  the  cell  size  to 
existing  data  from  actual  measurements  of  critical  tube  diameter  as  well  as 
-riticai  initiation  energies. 


2.  Experimental  Details 

Three  steel  tubes  (5,  15  and  30  cm  in  diameter)  of  length  £  20  m  were 
used  in  the  present  study.  The  tubes  were  equipped  with  ionization  probes 
for  the  measurement  of  the  detonation  velocity  and  piezoelectric  transducers 
for  pressure  measurements.  Commercial  grades  of  C0H2,  H2,  C2H4,  CoHg,  C3H8 
and  C4H10  and  bottled  compressed  air  were  used.  In  the  5  cm  and  T5  cm  dia¬ 
meter  tubes,  the  mixture  of  the  desired  composition  was  prepared  in  a  conti¬ 
nuous  flow  system  with  flow  rates  monitored  through  standard  calibrated  rota¬ 
meters.  The  tube  was  purged  with  the  premixed  gases  for  at  least  5  tube  vol¬ 
umes  prior  to  the  experiment.  In  the  larger  diameter  (30  cm)  tube,  the  mix¬ 
ing  procedure  was  as  follows:  the  tube  was  first  evacuated  and  the  desired 
volumes  of  air  and  fuel  were  then  introduced  into  the  tube  via  the  method  of 
partial  pressures.  A  bellows-type  pump  was  then  used  to  recirculate  the 
gases  from  one  end  to  the  other  to  permit  thorough  mixing  of  the  components. 
Initiation  of  the  detonation  was  via  an  exploding  wire  or  a  blasting  cap  de¬ 
pending  on  the  sensitivity  of  the  mixture.  In  the  30  cm  diameter  tube,  it 
was  often  required  to  use  an  additional  booster  charge  with  the  blasting  cap 
for  direct  initiation.  A  short  length  of  wire  spiral  was  usually  placed  at 
the  initiation  end  to  guarantee  the  rapid  formation  of  the  detonation  wave. 

To  record  the  detonation  cell  signatures,  smoked  aluminum  foils  were  used. 

The  aluminum  foils  were  usually  of  width  -rd  and  as  long  a  length  (in  the 
direction  of  propagation)  as  experimentally  possible  for  a  more  easier  esti¬ 
mation  of  the  averaged  cell  size.  For  cell  sizes  large  compared  to  the  dia¬ 
meter  of  the  pressure  transducer,  periodic  pressure  fluctuations  could  be 
seen  superimposed  on  the  main  pressure  trace.  The  periods  of  these  pressure 
fluctuations  "tc"  can  be  used  to  estimate  the  cell  length  "Lc"  (Lc  *  tr  Vrj) 
which  can  then  be  converted  into  cell  size  via  the  approximate  geometrical 
relationship  A  =  0.6  Lc.  Thus  the  pressure  trace  as  well  as  the  smoked  foil 
record  provide  two  independent  means  for  estimating  cell  sizes.  In  general, 
a  few  experiments  had  to  be  carried  out  for  each  mixture  composition  in  order 
to  conclusively  establish  the  cell  size  for  the  particular  mixture. 

2.  Resul ts 

The  averaged  detonation  cell  diameter  "A"  for  six  gaseous  fuels  (H2, 
C2H2,  C2H4,  C2Hg,  C3H0  and  04^0)  is  plotted  against  the  equivalence  ratio 
<t>  in  Fig.  1.  Except  for  hydrogen,  the  minimum  cell  size  (or  the  most  sensi¬ 
tive  composition)  corresponds  to  a  mixture  composition  slightly  on  the  fuel- 
rich  side  (<$)>!)  rather  than  stoichiometric  (<J>  =  1).  As  can  be  observed, 
estimates  of  the  cell  size  from  pressure  fluctuations  are  in  good  agreement 
with  those  obtained  from  smoked  foil  records.  Except  for  methane  (CH4),  all 
the  alkanes  (C3H8,  C2Hg,  C^n)  appear  to  have  the  same  sensitivity  in  that 
their  cell  sizes  are  practically  identical.  For  methane,  the  use  of  50  g  of 
explosive  charge  failed  to  cause  direct  initiation  for  the  stoichiometric 
composition  although  Kogarko9  reported  detonations  in  fuel -lean  (from  6.9  to 
8.2%  CH4)  and  fuel-rich  (from  11.1  to  13.5%  CH4)  C^-air  mixtures  in  a  simi¬ 
lar  tube  using  50  g  and  70  g  of  explosive  charges,  respectively.  However, 
v/e  feel  that  a  50  g  charge  is  the  upper  limit  for  safe  operation  in  a  uni¬ 
versity  laboratory.  Furthermore,  Moen10  recently  performed  experiments  in  a 
6  foot  diameter  tube  and  estimated  from  his  smoked  foil  record  a  cell  size 
of  about  33  cm  for  stoichiometric  methane-air  mixtures.  This  would  corres¬ 
pond  to  about  6  times  the  minimum  cell  size  of  about  5.35  cm  for  the  other 
alkanes  (i.e.,  C3H8,  C2HR  and  04^0)-  In  increasing  order  of  sensitivity 
(or  decreasing  cell  size),  ethylene  (C2H4)  follows  the  alkanes  and  hydrogen 
(H2)  is  slightly  more  sensitive  than  C^H^,  As  expected,  acetylene  is  found 
to  be  the  most  sensitive  fuel  with  a  minimum  cell  size  A  =  0.565  cm  as  com¬ 
pared  to  A  =  1 . 5  cm  for  H2,  A  =  2.6  cm  for  C2H4  and  A  =  5.35  cm  for  the  al¬ 
kane  group  (C3Hg,  C2Hg  and  C4H-|q).  For  the  stoichiometric  compositions  of 
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the  mixtures  studied  in  the  present  work,  cell  sizes  were  determined  by  Bull 
et  al.11  in  a  detonation  tube  of  rectangular  cross  section  (7.6  cm  x  3.8  cm) 
using  smoked  stainless  steel  foils.  Bull's  results  are  X  =  0.92  cm,  1.  cm, 
2.4  cm,  5.4  cm  and  5.6  cm  for  C0H2.  Hz,  C2H4,  C2  Hg  and  C3H0,  respectively 
and  are  in  good  agreement  with  the  present  observations.  It  is  also  of  im¬ 
portance  to  note  that  the  width  of  the  typical  U-shaped  curves  for  cell  size 
X  versus  fuel  composition  <j>  decreases  with  decreasing  sensitivity  of  the 
fuel.  Also,  the  increase  in  cell  size  as  the  equivalence  ratio  approaches 
the  rich  limit  (<p  >  1 )  is  in  general  much  slower  than  towards  the  lean  limit 
(<J>  <  1).  This  universal  behavior  is  reflected  in  practically  all  dynamic 
detonation  parameters  that  are  dependent  on  the  induction  time  of  the  mix¬ 
ture  itself. 

Based  on  cell  size  data,  the  critical  tube  diameter  dc  can  be  estimated 
from  the  empirical  relationship  dc  =  13X.  The  solid  curves  shown  in  Fig.  2 
represent  critical  tube  diameters  estimated  from  the  present  data  obtained 
for  X.  Although  an  extensive  experimental  measurement  of  the  critical  tube 
diameter  for  fuel -oxygen-nitrogen  mixtures  was  carried  out  by  Knystautas  et 
al.9,  relatively  few  data  exist  for  "dc"  in  fuel-air  mixtures,  in  particular 
for  off-stoichiometric  compositions.  Some  of  the  data  that  are  currently 
available  from  large-scale  field  tests  are  also  shown  in  Fig.  2.  Because  of 
the  interest  in  H2-air  mixtures  in  connection  with  nuclear  reactor  safety, 
large-scale  experiments  on  the  measurements  of  dc  for  both  lean  and  rich  Hz- 
air  mixtures  were  carried  out  recently  by  Benedick  at  Sandia  National  Labora¬ 
tory  (New  Mexico).12  As  can  be  observed  from  Fig.  2,  the  agreement  with  the 
estimates  of  dc  from  cell  size  data  via  the  empirical  law  dc  =  13X  is  ex¬ 
tremely  good  for  the  case  of  H2.  For  ethylene  (C9H4),  a  limited  number  of 
experiments  were  carried  out  at  Raufoss  (Norway)13  and  at  DRES14  for  lean 
Cz^-air  mixtures.  The  results  are  also  in  accord  with  the  present  esti¬ 
mates  from  cell  size  data.  For  lean  acetylene-air  mixtures,  some  direct  mea¬ 
surements  of  "dc"  were  also  carried  out  at  Raufoss.13  Although  no  special 
attempt  has  been  made  to  narrow  down  the  range  of  mixture  compositions  bet¬ 
ween  success  and  failure  for  transmission  in  a  given  tube,  the  results  are 
also  found  to  be  in  good  agreement  with  the  predictions  from  cell  size  data. 
For  stoichiometric  C2H2~air  and  H2-air  mixtures,  direct  measurements  of  "dc" 
were  carried  out  by  Knystautas  et  al.8  For  stoichiometric  Hz-air  mixtures, 
the  agreement  is  almost  perfect.  For  C2H2-air  mixtures,  the  measurements  of 
"dc"  are  slightly  higher  than  those  estimated  from  "X".  However,  "bottle" 
to  "bottle"  variations  of  the  purity  of  the  commerical  C^Hg  used  can  account 
for  the  observed  deviation  in  view  of  the  extreme  sensitivity  of  the  C9H2- 
air  mixture  itself.  It  would  be  of  interest  to  extend  the  critical  tube 
diameter  measurements  to  fuel-rich  hydrocarbon-air  mixtures  as  in  the  case 
of  hydrogen-air  mixtures  in  order  to  verify  the  validity  of  the  dc  -  13X 
correlation  over  the  entire  range  of  fuel  compositions  between  the  two 
1 imi ts . 

Although  the  critical  initiation  energy  has  been  deduced  from  a  know¬ 
ledge  of  the  critical  tube  diameter  from  energy  considerations  usina  a  "work 
done"  concept,5’15  a  better  model  was  recently  developed  by  Lee.12*16  The 
model  postulates  that  there  exists  a  minimum  surface  energy  per  unit  area  of 
the  wave  front  for  successful  transformation  of  a  planar  detonation  into  an 
unconfined  spherical  wave.  This  minimum  "surface  energy"  is  then  equated 
directly  to  the  area  of  the  critical  tube  (i.e.,  ndg2/4).  Assuming  the  re¬ 
quirement  of  a  minimum  size  of  a  detonation  kernel  for  direct  initiation, 
the  surface  area  of  the  critical  detonation  kernel  is  then  equated  to  the 
area  of  the  critical  tube  (i.e.,  4ttR£2  =  ndc74).  Using  strong  blast  theory 
to  relate  the  blast  energy  Ec  to  the  kernel  radius  R£  where  the  detonation 
is  assumed  to  be  of  Chapman-Jouguet  strength,  the  following  simple  expression 


was  obtained  by  Lee  et  al.16 


Ec  ■  4’V0Mc/I(T> 
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where  use  has  been  made  of  the  empirical  relationship  discussed  previously 
that  d<-  =  13A.  In  Eqs.  1  and  2,  y0  and  p,-,  denote  the  specific  heat  ratio  and 
the  initial  pressure  of  the  mixture,  M^j  is  the  Chapman-Oouguet  detonation 
Mach  No.  and  I  is  a  dimensionless  constant  representing  the  energy  integral 
in  strong  blast  theory.  Based  on  the  cell  size  data  obtained,  the  critical 
initiation  energies  for  the  various  mixtures  studied  can  be  evaluated  from 
Eq.  2.  Figure  3  shows  a  comparison  of  the  predicted  results  using  Eq.  2  and 
those  obtained  by  Elsworth  from  direct  experimental  measurements  of  the  cri¬ 
tical  weight.  For  C2H4,  Murray  et  al.4  have  also  obtained  a  few  critical  ini¬ 
tiation  charge  weights  for  lean  mixtures.  The  critical  energy  in  Fig.  3  has 
been  expressed  as  an  equivalent  charge  weight  of  tetryl  (1  g  tetryl  is  equi¬ 
valent  to  4270  joules)  for  easy  comparison  with  Elsworth's  experimental  data. 
The  agreement  in  general  is  reasonably  good  in  view  of  the  simplicity  of  the 
model.  No  experimental  data  of  critical  energies  have  been  reported  for 
C2H2-air  mixtures  to  permit  a  comparison  with  the  present  prediction  from  Eq. 
2.  For  most  practical  situations,  Eq.  2  can  be  used  to  predict  critical  ini¬ 
tiation  energies  (or  charge  weights)  with  quite  acceptable  accuracy. 


The  detonability  limits  in  tubes  have  also  been  related  directly  to  the 
cell  size  by  Lee  et  al.16  Based  on  the  experimental  work  of  Donato,  8»19  it 
was  found  that  the  onset  of  single-headed  spinning  detonations  correspond  to 
the  limit  for  stable  propagation  in  the  .qiven  tube.  Donato  showed  that  for 
mixtures  outside  the  single-head  limit,  a  finite  perturbation  would  cause  the 
detonation  to  fail  and  the  combustion  wave  will  continue  to  propagate  as  a 
deflagration  thereafter.  Only  for  mixtures  within  the  single-head  limit 
would  regeneration  into  the  detonation  mode  occur  shortly  downstream  of  the 
finite  perturbation  that  causes  the  wave  to  fail.  It  should  be  emphasized 
that  using  a  strong  enough  initiator,  spinning  detonations  can  always  be  ini¬ 
tiated  in  a  given  tube  over  a  wide  range  of  composition  beyond  the  onset 
limit  when  spinning  first  occurs.  However,  it  has  been  shown  that  these 
spinning  waves  will  decay  when  encountering  finite  perturbations  and  after 
failure  will  not  transit  again  to  the  detonation  mode  without  the  support  of 
a  strong  initiator.  However,  without  a  finite  perturbation  (such  as  a  couple 
of  turns  of  a  Schelkhin  spiral),  these  spinning  waves  may  persist  for  as  much 
as  100  tube  diameters  without  showing  any  sign  of  decay. 


Adopting  the  criterion  that  the  onset  of  single-head  spinning  waves  in  a 
given  tube  should  correspond  to  the  detonability  limits  for  that  tube,  the 
limits  in  circular  tubes  can  be  estimated  when  cell  size  data  are  available. 
Since  the  onset  of  single-head  spin  corresponds  to  a  cell  size  of  the  order 
of  the  tube  circumference  (i.e.,  A  =  ird),  the  composition  limits  for  a  given 
circular  tube  can  be  estimated  if  A  is  known.  Figure  4  gives  the  limiting 
tube  diameter  "d*"  as  function  of  equivalence  ratio  for  the  various  fuels  stu¬ 
died.  For  a  given  fuel  composition  (i.e.,  a  given  <f>),  stable  detonations  can¬ 
not  be  propagated  in  tubes  with  a  diameter  d  <  d*.  Thus,  for  the  case  of 
stoichiometric  Hp-air  mixtures,  d*  -  0.5  cm  when  compared  to  d*  -  0.2  cm  for 
C2H2  and  d*  =  1.7  cm  for  the  alkanes.  It  is  of  interest  to  note  that  accord¬ 
ing  to  the  cell  size  reported  by  Moen  for  methane  (i.e.,  X  -  33  cm),  a  mini¬ 
mum  tube  diameter  for  stable  detonation  propagation  in  stoichiometric  methane- 
air  mixtures  would  be  of  the  order  of  at  least  10  cm.  Thus  the  results 


reported  by  Wolanski  et  al.20  in  a  5  cm  tube  should  correspond  to  overdriven 
transient  waves  only.  This  is  supported  by  the  fact  that  they  observed 
single-head  spinning  waves  over  the  whole  range  of  fuel  concentrations 
studied.  For  two-dimensional  channels,  the  limit  should  correspond  to  a 
composition  where  the  channel  height  "W"  is  of  the  order  of  the  cell  size. 
This  is  supported  by  the  recent  experiments  of  Vasil iev21  who  found  that 
stable  detonation  propagation  corresponds  to  W  >,X. 

From  a  practical  point  of  view,  the  detonation  limits  for  a  thin  hori¬ 
zontal  layer  of  explosive  mixture  bounded  by  one  solid  surface  (i.e.,?  ground)  is 
of  interest.  Recent  experiments  by  Liu22  for  transmission  of  planar  detona¬ 
tions  through  rectangular  slits  into  unconfined  media  as  well  as  large  scale 
transmission  experiments  from  two-dimensional  channels  by  Benedick23  showed 
that  the  critical  channel  width  W*  =  3X.  Thus  it  is  reasonable  to  assume 
that  the  minimum  thickness  "h"  of  a  layer  of  explosive  gas  that  can  support 
a  detonation  should  be  of  the  order  of  one  and  a  half  cell  size  (i.e.,  h*  = 
1.5X).  Thus,  for  stoichiometric  H2-air  mixtures,  h*  ^  3  cm  as  compared  to 
the  alkanes  C7H3,  C^Hg  and  C4H1Q  for  which  h*  =  8  cm.  For  CH4,  h*  *  0.5  m 
based  on  Moenfs  estimate  of  X  =  0.33  cm,  which  is  much  smaller  than  commonly 
expected  for  methane.  Direct  experiments  must  be  performed  to  verify  the 
limit  postulate  of  an  unconfined  thin  layer  of  explosive. 

As  suggested  by  Schelkhin,24  the  cell  size  should  be  directly  propor¬ 
tional  to  the  induction  zone  thickness.  However,  no  theory  exists  at  pre¬ 
sent  whereby  the  cell  size  X  can  be  predicted  a  priori  from  a  knowledge  of 
the  induction  kinetics.  However,  Westbrook2 5-2 7  has  made  extensive  correla¬ 
tions  between  experimentally  measured  X  (or  equivalently  dc)  and  induction 
times  computed  from  detailed  kinetics  of  the  oxidation  processes.  In  gene¬ 
ral,  it  can  be  shown  that  X  'v  £  where  l  is  the  induction  length  but  the  con¬ 
stant  of  proportional i ty  has  to  be  determined  by  matching  with  an  experimen¬ 
tal  datum.  The  proportional i ty  constant^depends  on  how  the  induction  time 
is  computed.  For  example,  Westbrook  based  his  induction  time  on  a  constant 
volume  explosion  process,  while  Shepherd28  integrated  the  kinetic  rate  equa¬ 
tions  along  a  Rayleigh  line  as  in  the  classical  ZDN  model  for  the  detonation 
structure.  There  is  no  special  justification  for  either  method  since  both 
have  to  use  an  experimental  point  to  find  the  proportionality  constant  it¬ 
self.  The  solid  curves  in  Fig.  1  represent  the  estimated  cell  sizes  from 
kinetic  data  matching  at  the  stoichiometric  composition  point.  The  results 
are  relatively  in  good  agreement  with  the  data  despite  some  discrepancies 
as  the  mixture  composition  departs  from  the  stoichiometric  composition  which 
may  result  from  the  choice  of  the  matching  point.  Nevertheless,  the  reason¬ 
able  agreement  confirms  Shchelkhin's  postulate.  However,  useful  models  for 
linking  the  cell  size  X  to  the  kinetics  of  the  reaction  must  somehow  take 
into  consideration  the  complex  transient  gasdynamics  processes  of  the  cell 
interactions. 

4.  Conclusions 

From  the  cell  size  data  obtained  in  the  present  investigation,  critical 
tube  diameter,  critical  initiation  energy  (charge  weight)  as  well  as  detona- 
bility  limits  have  been  estimated.  Comparison  with  available  data  for  the 
critical  tube  from  direct  experimental  measurements  confirms  the  general 
validity  of  the  empirical  law  dc  =  13X  first  observed  by  Mitrofanov  and 
Soloukhin.  However,  more  extensive  experiments  should  be  carried  out  parti¬ 
cularly  in  fuel-rich  mixtures  to  consolidate  the  validity  of  this  empirical 
expression.  Due  to  its  fundamental  significance,  it  is  also  of  importance 
to  carry  out  experiments  in  fuels  other  than  the  hydrocarbons  to  verify  the 
universality  of  the  dc  ~  13X  law.  The  critical  initiation  ene,  lies  predicted 
bv  the  minimum  surface  energy  model  using  experimental  cell  size  data  agree 


reasonably  well  with  direct  experimental  measurements.  Thus,  apart  from  the 
practical  usefulness  of  Eq.  2,  the  linkage  of  the  cell  size  X  to  the  critical 
initiation  energy  Ec  should  elucidate  the  fundamental  mechanisms  of  the  fail¬ 
ure  and  reinitiating  processes  in  the  critical  tube  experiment,  the  blast 
wave  initiation  processes  as  well  as  the  detailed  cell  dynamics  of  the  trans¬ 
ient  shock  interactions.  This  may  lead  to  a  theory  for  the  empirical  law 
dc  -  13X  itself.  Prediction  of  detonability  limits  (or  minimum  tube  diame¬ 
ters  for  stable  propagation)  from  cell  size  data  is  an  important  contribution. 
However,  extensive  experiments  must  be  carried  out  further  to  verify  the  pos¬ 
tulates  for  stable  propagation  (i.e.,  d*  *  X/tt  for  circular  tubes,  W*  =5X  for 
two-dimensional  channels  and  h*  =  1.5X  for  a  thin  layer  bounded  by  one  solid 
surface).  Thus  far,  the  experiments  in  circular  tubes  by  Donato  and  two-di¬ 
mensional  channels  by  Vasil iev  are  of  limited  scope.  Although  it  has  been 
shown  that  all  the  dynamic  detonation  parameters  (i.e.,  critical  tube  diame¬ 
ter,  critical  initiation  energy  and  detonability  limits)  can  be  linked  to  the 
cell  size,  a  theory  for  the  prediction  of  the  cell  size  X  from  basic  informa¬ 
tion  on  the  physical  chemical  properties  of  the  explosive  mixture  is  still 
lacking. 
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Effect  of  Geometry  on  the  Transmission  of  Detonation 
through  an  Orifice 

Y.  K.  LIU*,  J.  H.  LEE,  and  R.  KNYSTAUTAS 

Department  of  Mechanical  Engineering.  McGill  University,  Montreal.  Canada 


An  experimental  study  has  been  carried  out  to  investigate  the  transmission  of  a  planar  detonation  wave  through  an  orifice 
into  an  unconfined  medium.  Mixtures  of  2H2  +  02  +  0N2  and  C2H4  +  3(02  +  0N2)  for  a  range  of  nitrogen 
concentrations  corresponding  to  1  s  /}  *  3.76  and  at  an  initial  presiue  of  I  atm  were  used  in  the  experiments.  It  is  found 
that  the  critical  diameter  for  the  transmission  through  an  orifice  is  identical  to  that  for  a  straight  tube  and  both  follow  the 
empirical  correlation  of  dc  —  13X.  The  transmission  through  square,  triangular,  elliptical,  and  rectangular  orifices  has 
led  to  the  development  of  a  correlation  based  on  the  effective  diameter  similar  to  the  case  of  circular  geometry  (i.e. ,  de„ 
21  1 3X).  The  effective  diameter  is  defined  as  the  mean  value  of  the  longest  and  shortest  dimensions  of  the  orifice  shape. 
The  effective  diameter  correlation  suggests  that  the  criterion  for  transmission  may  be  based  on  the  mean  curvature  of  the 
wave  front,  the  implication  being  that  it  is  not  to  exceed  a  certain  critical  value.  The  results  suggest  that  local  properties 
in  the  immediate  vicinity  of  the  wave  front  are  the  controlling  parameters  for  reinitiation  rather  than  the  properties  of  the 
gas  dynamic  flow  structure  in  the  wake.  Expressions  are  developed  to  provide  estimates  for  the  critical  transmission 
dimensions  for  arbitrarily  shaped  openings.  For  the  two-dimensional  limit  when  one  of  the  characteristic  linear 
dimensions  becomes  very  large  compared  to  the  other,  it  is  found  that  the  transmission  is  based  on  a  critical  value  for  the 
shorter  dimension  of  the  order  of  3  times  the  cell  diameter.  This  result  is  in  accordance  with  the  recent  large  scale 
experiments  of  Benedick  in  two-dimensional  channels  of  an  aspect  ratio  L  .W  as  large  as  35.  These  observed  results  can 
also  be  successfully  explained  in  terms  of  the  critical  wave  curvature  criterion. 


1.  INTRODUCTION 

l 

When  a  planar  detonation  wave  emerges  from  a 
circular  tube  into  an  unconfined  environment 
filled  with  the  same  explosive  gas  mixture,  the 
detonation  may  continuously  transform  to  propa¬ 
gate  as  a  spherical  wave  without  failure  if  the  tube 
diameter  is  greater  than  a  critical  value  (d  >  dc )  or 
may  fail  to  a  deflagration  wave  i i  d  <de.  Zeldo- 
vich  et  al.  [1]  were  the  first  to  suggest  that  the 
critical  diameter  is  related  to  the  induction  zone 
length  of  the  detonation  wave.  Mitrofanov  and 
Soloukhin  [2]  later  made  the  important  observa¬ 
tion  that  the  critical  tube  diameter  dc  a  13X  for 
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a  circular  tube  and  Wc  =  10X  for  a  rectangular 
channel  (Wc  is  the  channel  width  and  X  is  the  cell 
size  of  the  detonation  wave).  This  observation  was 
later  confirmed  by  Edwards  at  al.  [3] ,  who  also 
suggested  that  the  correlation  may  be  universal 
and  valid  for  all  mixtures.  The  empirical  correla¬ 
tion  dc  =  13X  gives  a  direct  link  between  the  de¬ 
tonation  transmission  phenomena  and  the  micro¬ 
structure  of  the  detonation  front  itself.  By  using 
this  correlation,  other  dynamic  detonation  proper¬ 
ties  of  mixtures  (such  as  the  critical  initiation 
energy  and  detonability  limits)  may  be  deduced 
[4-7] ,  On  the  practical  side,  dc  =  13X  is  also  use¬ 
ful  in  the  design  of  detonation  wave  traps  and  in 
estimating  the  critical  dimensions  of  pipes  and 
openings  between  interconnecting  chambers  to 
prevent  transmission  of  detonation  from  one  to 
the  other.  Recently,  considerable  effort  has  been 
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devoted  to  the  experimental  measurement  of  the 
critical  tube  diameter  for  a  wide  range  of  fuels  and 
initial  conditions.  It  was  found  that  the  empirical 
correlation  dc  =  1 3  A  indeed  holds  for  many  hydro¬ 
carbon  gas  mixtures  and  under  various  initial  con¬ 
ditions  [4,  5] . 

All  the  experiments  to  date  have  been  con¬ 
cerned  mainly  with  the  circular  geometry.  It  is  of 
great  interest  to  know  what  the  criterion  is  for  the 
transmission  of  a  detonation  through  an  orifice 
whose  geometry  is  other  than  circular.  Some 
preliminary  studies  on  the  detonation  transmission 
through  rectangular  orifices  have  already  been 
reported  [7] .  However,  systematic  investigations 
of  the  geometric  effect  of  the  opening  on  transmis¬ 
sion  of  detonation  have  not  been  made.  Results 
reported  by  Mitrofanov  [2]  and  Edwards  [3]  are 
based  on  square  channels.  For  true  two-dimen¬ 
sional  geometries,  rectangular  orifices  having  an 
aspect  ratio  of  infinity  (length/width)  must  be 
used.  For  L/W  -*  °°,  the  transmission  process  will 
be  two-dimensional  resulting  in  the  formation  of  a 
cylindrical  detonation.  The  present  paper  describes 
the  recent  results  of  an  experimental  study  on  the 
effect  of  geometry  of  the  opening  on  the  transmis¬ 
sion  of  detonation  waves. 

2.  EXPERIMENTAL  DETAILS 

The  apparatus  used  in  the  present  experiments  is 
essentially  similar  to  that  described  in  great  detail 


in  a  previous  paper  on  critical  tube  diameter 
measurements  [4],  Briefly,  as  shown  in  Fig.  1, 
the  apparatus  consists  of  a  tandem  arrangement  of 
detonation  tubes  extending  a  total  length  of  3.4  m 
and  culminating  in  a  maximum  tube  diameter  of 
20  cm  diam.  The  20  cm  diam  tube  is  joined  to  a 
cylindrical  detonation  chamber  (56  cm  diam  X 
78  cm  long).  At  the  interfacing  flanges  between 
the  detonation  tube  and  chamber,  provision  is 
made  for  inserting  orifice  plates  which  are  typi¬ 
cally  1  cm  thick  X  22  cm  diam  with  a  differently 
shaped  orifice  opening  milled  into  each  plate. 
However,  the  maximum  linear  opening  dimension 
for  any  shape  of  orifice  is  limited  to  the  inner 
diameter  of  the  detonation  tube,  i.e.,  20  cm.  The 
orifice  openings  have  sharp  edges  formed  by  cham¬ 
fering  at  45°.  A  planar  detonation  wave  is  initiated 
at  the  extreme  end  of  the  detonation  tube  arrange¬ 
ment  by  a  powerful  ignition  source  using  an  elec¬ 
trical  condenser  bank  connected  to  an  ■exploding 
wire  igniter.  The  planar  detonation  wave  propa¬ 
gates  the  3.4  m  toward  the  orifice  and  emerges 
into  the  detonation  chamber  by  diffracting  through 
the  orifice.  Before  each  experiment,  the  tube- 
chamber  combination  is  first  evacuated  (p0  <  1 
Torr)  and  then  filled  to  1  atm  and  flushed  five 
volumes  by  the  flow  displacement  technique.  The 
flow  composition  is  controlled  by  calibrated  ro- 
tometer-type  flow  meters,  and  typical  flow  rates 
are  of  the  order  of  30-40  1/min.  The  fuels  studied 
were  acetylene,  hydrogen,  and  ethylene  of  com- 
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mercial  purity  (CP  grade).  The  explosive  mixture 
consists  of  stoichiometric  proportions  of  fuel- 
oxygen  with  nitrogen  dilution  as  specified  by  the 
parameter  (3  3  N2/02.  The  range  for  (3  is  1  <  (3  < 
3.76,  the  upper  bound  corresponding  to  the 
fuel-air  composition.  All  experiments  were  carried 
out  at  1  atm  initial  pressure. 

To  study  the  geometric  effect  of  the  detonation 
transmission  phenomenon,  three  groups  of  orifices, 
each  group  corresponding  to  the  same  orifice  area, 
were  used.  Figure  2  illustrates  the  orifice  configu¬ 
rations  for  each  group  and  consists  of  circular, 
square,  triangular,  and  elliptical  (2 : 1  axis  ratio) 
orifices.  The  three  groups  of  orifices  have  areas  of 
20,  64,  and  127  cm2,  corresponding  to  circular 
orifice  dimensions  of  5,  9,  and  12.6  cm  diam.  In 
order  to  study  the  detonation  transmission  process 
for  a  two-dimensional  expansion,  a  series  of 
rectangular  orifices  was  used.  The  range  of  aspect 
ratio  (L/W)  for  the  rectangular  orifices  is  3.75  < 
L/W  <  15,  with  the  range  of  actual  dimensions 
being  8  <  L  <  19  cm  and  1  <  W  <  5  cm.  To  assess 
the  detonation  transmission  from  the  detonation 
tube  through  the  orifice  and  into  tire  detonation 
chamber,  two  calibrated  piezoelectric  pressure 
transducers  (PCB113A24,  5  mv/psi)  were  used, 
one  in  the  detonation  tube  and  the  other  at  the 
center  of  the  back  flange  (opposite  to  the  detona¬ 
tion  tube,  exit)  of  the  detonation  chamber.  Suc¬ 
cessful  transmission  into  the  chamber  or  failure  of 
the  detonation  wave  in  the  chamber  was  established 
by  the  time  of  arrival  and  the  magnitude  of  the 


Fig.  2.  Comparison  of  different  orifice  geometries  for  the 
same  area. 


reflected  pressure  as  monitored  by  the  transducer 
at  the  back  flange  of  the  chamber.  A  smoked  foil, 
placed  along  the  inner  wall  of  the  20  cm  diam  de¬ 
tonation  tube,  was  used  to  determine  the  charac¬ 
teristic  cell  diameter  X  for  each  mixture. 

3.  RESULTS  AND  DISCUSSIONS 

The  preliminary  results  reported  earlier  [7) 
have  indicated  that  the  critical  diameter  for 
detonation  transmission  from  a  tube  and  an  orifice 
is  the  same  under  similar  conditions.  However,  due 
to  a  higher  upstream  pressure  from  the  reflection 
of  the  detonation  wave  from  the  orifice  plate,  one 
may  expect  that  the  critical  orifico  diameter 
should  be  somewhat  less  than  that  for  a  straight 
tube.  Since  the  results  reported  earlier  are  based 
on  a  limited  number  of  large  scale  experiments 
under  field  conditions  where  experimental  param¬ 
eters  could  not  be  controlled  accurately,  a  series  of 
test  were  carried  out  in  H2-02-N2  and  C2H4- 
02-N2  mixtures  at  an  initial  pressure  of  1  atm  to 
verify  this  important  result.  Figures  3  and  4  show 
the  critical  orifice  diameter  for  a  range  of  nitrogen 
concentration  in  stoichiometric  mixtures  [i.e., 
2H2  +  02  +  0N2  with  1  <  /3  <  3.76  and  C2H4  + 
3(02  +  0N2)  for  1  <  0  <  3.0).  The  results  are 
found  to  be  in  perfect  accord  with  those  for 
straight  circular  tubes  obtained  previously  by 
Knystautas  et  al.  [4] .  Cell  size  measurements  have 
also  been  carried  out,  and  the  cell  size  data  are  also 
in  agreement  with  those  of  Knystautas  et  al.  Based 
on  the  cell  size  data,  the  empirical  correlation  of 
dc  2:  13X  is  also  found  to  be  valid  for  circular 
orifices  as  it  is  for  straight  circular  tubes.  Thus  the 
present  results  conclusively  demonstrate  the 
equivalence  of  the  critical  diameter  for  orifices  and 
for  tubes,  even  though  in  certain  experiments,  the 
ratio  of  the  orifice  diameter  to  the  tube  diameter 
in  which  the  orifice  plate  is  mounted  is  quite  small 
(d/D  ^  0.1).  One  would  expect  significant  in¬ 
fluence  of  the  reflection  of  the  detonation  wave 
on  the  flow  structure  of  the  product  gases  emitting 
from  the  orifice  under  this  condition.  Thus  the 
present  experiments  indicate  that  the  reinitiation 
process  as  the  detonation  wave  exits  from  the 
orifice  (or  tube)  is  a  very  local  phenomenon 
limited  to  a  small  region  of  the  order  of  a  cell 
length  in  thickness  at  the  detonation  front  itself. 
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Dilution  Ratio  6  ■  N^/O^ 


Fig.  3.  Comparison  of  the  critical  tube  and  orifice  diameter*  as  a  function  of  dilution 
ratio  in  H2-02-N2  mixtures. 

The  flow  behind  the  front  seems  to  play  a  minor  being  equated  to  the  work  done  by  the  product 

role  in  the  reinitiation  process.  Hence,  the  so-  gases  emitted  from  the  tube  (or  orifice),  does  not 

called  work-done  model  suggested  by  Lee  and  appear  to  give  the  correct  description.  However, 

Matsui  [8]  in  which  the  reinitiation  is  equated  to  the  later  surface  energy  model  proposed  by  Lee 

the  blast  initiation  process,  with  the  blast  energy  [5] ,  where  a  minimum  surface  energy  is  postu- 


Dilution  R* tie  . 

Fig.  4.  Comparison  of  the  critical  tube  and  orifice  diameters  as  a  function  of  dilution 
ratio  in  C2H4-02-N2  mixtutes. 
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lated  to  be  required  for  the  evolution  of  a  planar 
to  a  spherical  wave,  appears  to  be  more  compatible 
with  the  present  findings. 

For  the  transmission  of  detonation  through  ori¬ 
fices  with  shapes  other  than  circular,  a  series  of 
experiments  has  been  carried  out  for  square, 
triangular,  and  elliptical  orifices  for  the  two  mix¬ 
tures  2H2  +  02  +  (3N2  for  1  <  (J  £  3.76  and 
C2H«  +  3(02  +  /3N2)  for  1  <  0  $  3  at  initial 
pressures  of  1  atm.  Attempts  to  use  the  hydraulic 
diameter  ( dH  =  4  X  area/wetted  perimeter)  to 
characterize  the  different  geometrical  openings 
and  seek  a  similar  correlation  between  the  critical 
hydraulic  diameter  and  the  cell  size  of  the  mixture 
as  was  done  in  the  case  of  the  circular  geometry 
(i.e.,  dc  -  13X)  failed  to  yield  meaningful  result*. 
Figures  5  and  6  give  the  critical  hydraulic  diameter 
dH  for  the  H2  and  C2H4  mixtures  for  various 
values  of  nitrogen  concentration  for  the  square, 
triangular,  and  elliptical  orifices.  Each  group  of 
points  corresponds  to  the  same  orifice  area.  As  can 
be  observed,  for  the  same  orifice  area,  the  circular 
orifice  gives  the  largest  critical  hydraulic  diameter, 
or  for  the  same  hydraulic  diameter,  the  circular 
orifice  requires  the  most  sensitive  mixture  (i.e., 
the  smallest  value  of  0).  The  triangular  orifice 
appears  to  be  the  most  effective  in  the  transmis- 
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sion  of  the  detonation  for  the  same  value  of  the 
hydraulic  diameter,  being  distinctly  more  effective 
than  the  square  or  the  elliptical  orifices.  This 
suggests  that  the  critical  hydraulic  diameter  is 
geometry  dependent.  No  universal  correlation 
is  observed  for  the  differently  shaped  orifices, 
and  the  hydraulic  diameter  does  not  appear  to  be 
the  appropriate  correlation  parameter  to  describe 
the  detonation  transmission  process.  However,  it  is 
found  that  one  can  define  an  effective  diameter 
del t  based  on  the  mean  value  of  diameter  of  the 
inscribed  d2  and  the  circumscribed  d2  circles  for 
each  of  the  geometrical  openings  [i.e.,  dtU  = 
(d i  +  d2)l 2] ,  and  this  effective  diameter  appears 
to  follow  the  same  universal  empirical  law  of 
d,tt  =  13X  as  for  the  case  of  the  circular  orifice 
geometry.  The  results  for  d, ff/X  plotted  against 
the  hydraulic  diameter  for  the  three  geometries 
(i.e.,  square,  triangular,  and  elliptical  orifices)  for 
both  H2  and  C2H4  mixtures  are  given  in  Figs. 
7-9.  As  indicated  in  these  figures,  the  experi¬ 
mental  results  show  that  d,tt  a  13X  for  all  the 
geometries.  That  the  effective  diameter  d,tt  =* 
(di  +  d2)/2  is  found  to  be  equivalent  to  the 
diameter  of  a  circular  orifice  suggests  that  the 
reinitiation  process  depends  on  the  curvature  of 
the  wave  front.  Upon  emerging  from  the  orifice 


Fig.  5.  Variation  of  the  critical  hydraulic  diameter  with  dilution  ratio  in  H2-02-N2 
mixtures. 


Y.  K.  LIU  ET  AL. 


Dilution  Rotio  0  •  Nj/Oj 

Fig.  6.  Variation  of  the  critical  hydraulic  diameter  with  dilution  ratio  in  C2H4-O2-N2 
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Fig.  7.  Variation  of  the  effective  diameter  with  hydraulic  diameter  for  the  square  orifice. 
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Fig.  9.  Variation  of  the  effective  diameter  with  hydraulic  diameter  for  the  elliptical 
orifice. 


(or  tube),  the  planar  detonation  wave  front  is  sub¬ 
jected  to  an  expansion  process.  The  penetration 
of  the  rarefaction  waves  toward  the  axis  of  the 
detonation  gives  rise  to  a  curved  front  as  the 
transformation  to  an  eventual  unconfined  spheri¬ 
cal  detonation  takes  place.  The  radius  of  curvature 
depends  on  the  linear  dimension  of  the  front  as 
well  as  on  the  geometry.  For  a  circular  orifice, 
there  is  only  one  linear  dimension  (i.e.,  the  tube  or 
orifice  diameter),  and  the  existence  of  a  critical 
diameter  indicates  a  corresponding  critical  radius 
of  curvature  for  successful  propagation  of  the 
curved  front.  For  noncircular  geometries ,  there 
exist  two  linear  dimensions  characterizing  the 
initial  wave  front  shape.  The  effective  diameter 
being  the  average  of  the  inscribed  and  circum¬ 
scribed  circles  gives  the  appropriate  length  scale 
on  which  the  mean  radius  of  curvature  of  the 
curved  detonation  front  depends.  Thus  one  may 
expect  a  similar  correlation  of  dtlt  =  13X  if  a 
certain  critical  radius  of  curvature  is  required  for 
the  successful  transformation  of  an  eventual 
spherical  wave.  The  results  for  different  geometries 
and  the  successful  correlation  of  datt  =  13X  as 
explained  above  are  compatible  with  the  conclu¬ 
sions  stated  earlier  that  the  reinitiation  process  is 
a  very  local  phenomenon,  limited  to  the  region  in 
the  immediate  vicinity  of  the  detonation  front, 
with  little  dependence  on  the  flow  structure  of  the 
product  gases  emitting  from  the  opening  or  tube. 
The  key  parameter  that  characterizes  the  wave 
front  appears  to  be  the  mean  curvature,  which 
depends  strongly  on  the  initial  geometry  of  the 
detonation  front  as  it  exits  the  orifice  and  is 


subjected  to  transverse  penetration  by  the  rare¬ 
faction  waves. 

Generalizing  the  above  discussions  to  arbitrary 
geometries,  one  may  consider,  for  example,  a 
polygonaily  shaped  orifice  of  n  sides,  with  /  being 
the  length  of  one  side.  The  hydraulic  diameter  dH 
for  the  polygon  is  thus  (Fig.  10) 

4A  TT 

dH=-=lcot--  (1) 

nl  n 

The  relation  between  dtf1  and  the  hydraulic  diam- 


Rg.  10.  Schematic  illustration  of  polygon  geometry. 
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eter  dH  is  therefore 

1  +  cos  (ir/n)  , 

deft  ~  .  7  .  dH.  (2) 

2  cos  (ir/n) 

Thus  if  we  use  the  present  experimental  results 
obtained  for  dttt  =  13X,  the  critical  hydraulic 
diameter  for  the  polygonal  shaped  orifice  becomes 

2  cos  (n/n) 

=  ~ - 7T  13X.  (3) 

1  +  cos (n/n) 

For  n  -*  “  the  above  expression  yields  dn  ■*  13X 
for  the  circular  orifice.  The  general  expression 
represented  by  Eq.  (3)  can  now  be  used  to  corre¬ 
late  the  experimental  data  for  the  square  ( n  “  4) 
and  the  triangular  (n  =  3)  orifices.  Table  1  gives  a 
comparison  between  the  prediction  given  by  Eq. 
(3)  and  the  experimental  results.  The  agreement  is 
quite  good,  and  it  would  be  of  interest  to  extend 
the  experiments  to  pentagonally  and  hexagonally 
shaped  orifices  to  test  the  general  validity  of  Eq. 

(3). 

For  the  ellipse,  the  effective  diameter  is  related 
to  the  hydraulic  diameter  by  the  expression 


64  -  3 B* 
64  -  16B2 


d,n  =  iU+*)  (1  + 


where  K  =  a/b  is  the  ratio  of  the  major  (a)  to  the 
muior  (6)  axes  and  B  =  (K  —  1)I(K  +  1).  Thus 
using  the  correlation  d,tl  =  13X,  we  obtain  for  the 
critical  hydraulic  diameter  for  the  ellipse 


(1  +/Q(I  +  1/AT)[(64  -  35*)/(64  -  16 B3)] 


For  a  *  b  ( a  circular  orifice),  K  -  1 ,  and  B  ■  0,  the 
above  expression  yields  dH  =  de  ■  13X.  Table  2 
compares  the  predicted  results  using  Eq.  (5)  with 
experimental  data.  The  agreement  is  fair  and 
shows  at  least  the  correct  qualitative  trend  of  a 
decreasing  da  as  K  increases.  Again  the  general 
validity  of  Eq.  (5)  requires  further  experimental 
work  using  elliptical  orifices  with  different  values 
of  K. 

Rectangular  orifices  of  different  aspect  ratios 
(L/W  ■  length  of  ilot/width  of  slot)  are  of  partic¬ 
ular  interest  since  the  variation  of  L/W  from  one 
to  infinity  gives  the  continuous  transition  from 
three-  to  two-dimensional  geometry,  correspond¬ 
ing  to  the  detonation  evolving  to  a  cylindrical 
wave.  Previous  experiments  by  both  Mitrofanov 
and  Soloukhin  [2]  and  Edwards  et  al.  [3]  have 
reported  that  for  square  channels  Wc  at  10X.  This 
agrees  with  the  present  results  for  the  square  ori¬ 
fice,  where  the  critical  width  We  is  found  to  be  in 
the  range  10X  <  Wc  S  1 IX  for  both  the  H2  and 
C2H4  mixtures  tested.  Results  for  H'./X  for  values 
of  L/W  ranging  from  1  to  15  are  shown  in  Fig.  ■ 
11.  Also  plotted  in  Fig.  11  are  the  earlier  results 
of  Bjerketvedt  et  al.  [7] ,  who  have  performed  a 
few  experiments  with  both  C2H2  and  C2H4 
mixtures  for  orifices  with  1  S  L/W  5  5.  The 
results  show  that  Wc/\  decreases  from  a  value 
of  the  order  of  10  for  L/W  as  1  quite  rapidly  to  an 
asymptotic  value  of  about  3  for  L/W  £  7.  In 
other  words,  as  one  approaches  the  two-dimen¬ 
sional  geometry  (L/W  -*  «>),  the  critical  channel 
width  -*  3X.  This  rather  surprising  observation 
has  prompted  Benedick  et  al.  [9]  to  perform  some 
large  scaled  field  tests  in  which  two-dimensional 
channels  (Instead  of  rectangular  slots)  of  very  large 
aspect  ratios  ( L/W- *  35)  could  be  used.  The  recent 
results  obtained  by  Benedick  for  H2-air  mixtures 
are  also  plotted  in  Fig.  1 1 ,  which  confirms  the 
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Fig.  1 1 .  Variation  of  the  critical  slot  width  with  aspect  tatio  foi  tectang-ulai  orifices. 


asymptotic  limit  value  of  WJ X  a  3  that  has  been 
obtained  in  the  present  study  using  slot  orifices. 

To  attempt  to  explain  all  these  observed 
detonation  transmission  phenomena,  and  espe¬ 
cially  the  surprising  result  that  the  critical  trans¬ 
mission  criterion  for  rectangular  orifices  changes 
from  WC[K  a  10  for  an  aspect  ratio  LjW  =  1 
(square  orifice)  to  WC[K  a  3  for  aspect  ratios 
L/W  >  7,  one  can  invoke  the  recently  proposed 
wave  curvature  concept  of  Lee  [10].  In  the  gen¬ 
eral  transmission  situation,  Lee  suggests  that  when 
a  planar  detonation  wave  propagating  in  a  confining 
linear  tube  or  channel  emerges  suddenly  into  an 
unconfined  volume,  it  is  subject  to  lateral  relief 
due  to  the  centered  expansion  wave  at  the  edge  of 
the  opening.  The  rarefaction  waves  that  are  thus 
generated  penetrate  into  the  wave  front,  retard  the 
affected  portions,  and  thereby  essentially  impose  a 
curvature  on  the  wave  front.  He  then  defines  the 
critical  condition  for  transmission  in  terms  of  a 
minimum  local  radius  of  curvature  of  the  dif¬ 
fracted  wave,  namely:  If  the  rarefaction  waves 
give  rise  to  a  curvature  k  of  the  diffracted  detona¬ 
tion  wave  exceeding  a  certain  critical  value,  then 
failure  of  the  wave  results. 

For  example,  this  critical  curvature  can  be 


determined  in  the  three-dimensional  case— say,  for 
the  critical  tube  diameter  situation— and  the 
de  =  13X  correlation  can  then  be  explained.  Con¬ 
sider  a  planar  detonation  wave  characterized  by  a 
hydrodynamic  thickness  AH  emerging  through  the 
tube  opening  into  the  unconfined  environment. 
The  hydrodynamic  thickness  AH  represents  the 
equivalent  thickness  of  a  "one-dimensional” 
detonation  wave  where  equilibrium  conditions 
prevail  at  the  C- J  plane  situated  at  a  distance  of 
one  hydrodynamic  thickness  from  the  leading 
front.  The  hydrodynamic  thickness,  as  reported 
by  Vasiliev  et  al.  [11]  and  Edwards  et  al.  [12], 
ranges  from  about  3  to  5  cell  lengths  in  extent. 
Taking  the  observed  relation  between  the  cell 
length  I  and  the  cell  width  X  as  I  «  1.6X  gives  AH 
equal  to  about  5-8  cell  widths  X.  Taking  an  aver¬ 
age  value,  one  can  relate  AH  and  X  as  AH  ^  6.5X. 
With  regard  to  transmission,  it  is  plausible  to 
postulate  that  if  the  rarefaction  waves  do  not  reach 
the  core  of  the  fully  exposed  wave  (i.e.,  AH  thick) 
by  the  time  the  wave  has  propagated  a  distance 
Ah  ,  then  the  detonation  wave  will  survive.  In 
other  words,  for  successful  transmission  the 
detonation  must  have  an  unattenuated  portion  at 
least  of  the  order  of  one  hydrodynamic  thick- 
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ness.  This  is  compatible  with  the  concept  advanced 
by  Urtiew  and  Tarver  [13]  of  an  unattenuated 
core.  For  the  unattenuated  detonation  core  to  be 
of  the  order  of  AH  s  6.5X  (as  contrasted  with 
X/2 ,  proposed  by  Urtiew  on  the  basis  that  the  core 
represents  the  minimum  constant  area  tube  diam¬ 
eter),  the  detonation  front  would  have  to  propa¬ 
gate  a  distance  of  2AH.  In  other  words,  this  im¬ 
plies  that  the  wave  AH  thick  having  traveled  the 
distance  AH  is  now  2AH  beyond  the  tube  exit  and 
that  it  has  traversed  this  distance  by  traveling  at 
the  C-J  velocity  D.  Therefore,  the  time  of  detona¬ 
tion  travel  rD  beyond  the  tube  exit  can  be  written 
as  rD  =  2A n/D.  On  the  other  hand,  the  time  for 
the  radial  penetration  of  the  rarefaction  waves  is 
fR  -  ^c/(2/ci),  where  ^  is  the  acoustic  speed  in 
the  product  gases.  Thus,  the  criterion  for  success¬ 
ful  transmission  is  that  tq  must  be  just  less  than 
tr,  or  in  the  limit  rD  =  rR,  which  leads  to 


2Ah 

D 


2c , 


For  most  detonable  mixtures  D  ^  2^,  hence  dc  Ss 
2Ah,  and  therefore  dc  =  13X,  which  is  in  perfect 
agreement  with  experiment.  It  follows  that  the 
local  radius  of  curvature  at  the  forward  tip  of  the 
diffracted  wave  is  R  a  AH  =  dj 2  -  6.5X  and  can 
be  taken  as  representative  of  the  curvature  im¬ 
posed  by  a  three-dimensional  expansion. 

Analogous  reasoning  would  be  applicable  to, 
say,  a  square  orifice  with  some  equivalent  radius  of 
curvature  proportional  to  the  average  of  the  mini¬ 
mum  and  maximum  radial  dimensions  of  the 
square  opening, 


1 


R  =  - 
2 


W  +  s/2W 


2 


6.5X, 


which  immediately  yields  (Vc/X  =  10.8,  which  is  in 
accord  with  the  present  experimental  results  an 
and  the  results  of  Mitrofanov  [2]  and  Edwards 
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One  can  extend  the  above  arguments  to  the 
two-dimensional  case.  In  the  case  of  the  rectang¬ 
ular  orifice  of  large  aspect  ratio  ( L/W  >  7),  the 


planar  wave  diffracting  through  the  opening  trans¬ 
forms  into  a  cylindrical  wave.  For  the  same  critical 
curvature,  by  simple  geometric  considerations, 
namely, 


1  <14 


'cyl 


'sph 


A  dR 


one  can  easily  show  that  the  radius  of  curvature 
for  the  equivalent  (i.e.,  same  curvature  x)  cylindri¬ 
cal  wave  would  be  one-half  that  for  the  spherical 
case,  so  that 


Rcyi  =  =  —  =  3.25X  =  Wc, 

2  2  c 


which  again  is  in  agreement  with  the  present  re¬ 
sults.  For  the  intermediate  values  of  L/W,  it  is  of 
interest  to  see  whether  the  correlation  dm  tt  x  13X 
as  given  by  Eq.  (3)  for  the  polygon  applies  as  well. 
For  a  rectangular  slot,  the  effective  diameter  is 
related  to  the  length  L  and  the  width  W  of  the  slot 
as  follows: 


d.rr  =  i(fc'  +  v/h'4+l2). 


(6) 


Replacing  dt  tt  by  13X,  the  ratio  Wc/X  can  be  ob¬ 
tained: 


W e 


26 


X  1  +  sJ{L/W)2  + 


(7) 


The  above  expression  gives  the  limit  as  L/W  1 
of  Wc/X  =  10.8  as  already  discussed  for  the  square 
orifice.  For  the  value  of  L/W  -  7,  Eq.  (7)  yields 
WJ\  =  3.2,  which  is  in  accord  with  the  experi¬ 
mental  result  of  L/W  =  7,  at  which  point  the 
approach  to  the  asymptotic  value  is  observed  to 
occur.  Predictions  from  Eq.  (7)  are  shown  by  the 
dotted  line  in  Fig.  1 1  and  appear  to  give  the  cor¬ 
rect  trend  for  the  intermediate  values  of  L/W  in 
the  approach  to  the  asymptotic  limit.  Beyond 
L/W  m  7,  the  curve  based  on  Eq.  (7)  begins  to 
deviate.  This  is  not  surprising  as  now  two-dimen¬ 
sionality  is  being  approached  and  Eqs.  (6)  and  (7) 
are  only  valid  in  the  context  of  the  three-dimen¬ 
sional  geometry. 
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4.  CONCLUSIONS 

The  present  study  of  the  transmission  of  a  planar 
detonation  from  orifices  of  different  geometrical 
shapes  into  an  unconfined  medium  has  yielded  a 
number  of  interesting  results.  The  fact  that  the 
transmission  through  an  orifice  (where  the  reflected 
detonation  from  the  orifice  plate  gives  rise  to 
higher  stagnation  pressures  and  hence  influences 
the  structure  of  the  downstream  flow)  has  the 
same  requirement  as  that  from  a  straight  tube 
indicates  that  the  reinitiation  phenomenon  is  a 
very  local  process,  occurring  in  the  immediate 
vicinity  of  the  detonation  front  itself.  Thus, 
properties  at  the  wave  front  (such  as  the  wave 
curvature)  play  a  more  important  role  than  the 
gas  dynamic  processes  of  the  product  gases  in  the 
wake.  The  results  for  the  transmission  through 
different  geometrically  shaped  orifices  lead  to  the 
concept  of  an  effective  diameter  dett,  which  is 
the  mean  value  of  the  longest  and  shortest  charac¬ 
teristic  dimensions  of  the  orifice.  The  universal 
correlation  based  on  this  effective  diameter 
d»u  —  13X  suggests  that  the  ability  of  a  detona¬ 
tion  to  evolve  without  a  failure  requires  a  certain 
minimum  limit  for  the  wave  front  curvature.  In 
other  words,  as  the  wave  emerges  from  the  tube  or 
orifice  and  is  subjected  to  lateral  expansion,  the 
initially  planar  front  will  become  curved  as  the 
rarefaction  waves  penetrate  into  the  core  of  the 
detonation.  Failure  occurs  when  the  mean  curva¬ 
ture  of  the  attenuating  detonation  exceeds  a  cer¬ 
tain  limit.  The  mean  curvature  being  proportional 
to  the  average  linear  dimensions  of  the  wave  re¬ 
sults  in  the  successful  correlation  of  dett  s*  13X. 
Based  on  the  effective  diameter  correlation, 
formulas  have  been  developed  to  provide  the  trans¬ 
mission  criteria  for  arbitrarily  shaped  orifices  (i.e., 
polygonal,  elliptical,  rectangular),  and  experi¬ 
mental  results  agree  with  these  general  expressions, 
Of  particular  interest  in  the  present  study  are  the 
results  of  the  two-dimensional  limit  when  one 
characteristic  dimension  becomes  much  larger  than 
the  other.  Experiments  conclusively  indicate  that 
the  two-dimensional  limit  gives  the  shorter  charac¬ 
teristic  dimension  a  critical  limiting  value  of  about 
three  cell  diameters  (i.e.,  Wc  —  3X).  This  limiting 


value  becomes  independent  of  the  longer  dimen¬ 
sion  when  the  ratio  L/WiZ  7.  The  observed  change 
in  the  transmission  criterion  from  We/X  st  10  for 
the  square  orifice  to  Wc/X  at  3  for  a  rectangular 
orifice  of  large  aspect  ratio  ( L/W  >  7)  can  be  ex¬ 
plained  in  terms  of  the  recently  proposed  wave 
curvature  model  of  Lee.  Lee's  model  postulates 
a  minimum  critical  wave  curvature  criterion  for 
transmission.  Since  for  the  same  curvature  the 
associated  radius  of  curvature  is  less  for  a  cylindri¬ 
cal  wave  than  a  spherical  wave,  such  direct  geo¬ 
metric  considerations  clearly  confirm  the  experi¬ 
mentally  observed  effects. 
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ABSTRACT 

Results  of  large-scale  experiments  on  the  transmission  of 
hydrogen-air  and  ethylene-air  detonations  from  a  two-dimensional 
channel  into  a  large  volume  are  reported.  These  data  cohfirm  the 
surprising  observations  of  Liu  et  al.  who  reported  recent  laboratory- 

scale  results  which  indicate  that  successful  transmission  of  a 
detonation  emerging  from  a  rectangular  orifice  with  a  large  aspect 
ratio  L/w  >  5,  requires  that  the  minimum  width  of  the  channel  need 
be  only  about  three  detonation  cell  diameters,  i.e.  Wc  «  3i.  These 
results  are  surprising  in  view  of  the  previous  findings  of  SolouXhin 
and  Edwards,  where  for  square  channels  (L/W  -  1),  they  observed  that 
Wc  -  10X. 

The  apparatus  for  the  large-scale  experiments  consisted  of  a 
two-dimensional  confined  channel  (1.83  m  wide  x  2.44  m  long)  with 
a  variable  width  orifice,  to  the  end  of  which  was  attached  an  in¬ 
flated  polyethylene  bag  of  much  larger  dimensions  to  provide  an  un¬ 
confined  environment  for  detonation  transmission.  Observations  were 
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made  with  high  speed  cinematography  and  with  smoked  foils  to  measure 
the  detonation  cell  size.  Our  results  indicate  that  a  minimum  chan¬ 
nel  width  of  3X.  is  required  for  successful  transmission  from  channels 
whose  aspect  ratios  varied  from  34  to  8.  For  channel  aspect  ratios 
of  less  than  about  8,  the  minimum  channel  width  for  successful 
detonation  transmission  increases  continuously  from  3^  to  10\  as 
the  aspect  ratio  approaches  1. 


INTRODUCTION 


During  the  past  few  years,  there  have  been  many  investigations 
to  determine  the  critical  tube  diameter  for  the  transmission  of 
confined  planar  detonations  to  unconfined  spherical  detonations. 

Both  laboratory  experiments  (Matsui  and  Lee  1978;  Knystautas  et  al. 
1982)  as  well  as  large  scale  field  tests  have  been  carried  out  for 
a  variety  of  hydrocarbon-air  and  oxygen  enriched  mixtures  over  a 
range  of  initial  pressures*.  The  importance  of  the  critical  tube 
diameter  is  that  it  is  directly  related  to  the  dynamic  detonation 
parameters  such  as  critical  initiation  energy  and  detonability 
limits  (Lee  et  al.  1982).  Knowledge  of  the  critical  tube  diameter 
permits  these  important  dynamic  parameters  of  the  explosive  to  be 
estimated.  The  critical  tube  diameter,  dc,  is  also  found  to  be 
directly  related  to  the  cell  size,  \,  of  the  detonation  front  via 
the  simple  empirical  formula,  dc  »  13  This  observation  was  first 

made  by  Mitrofanov  and  Soloukhin  (1965)  who  observed  that  for  low 
pressure  acetylene-oxygen  detonations  in  circular  tube,  dc  =  13\. 

The  significance  of  this  went  unnoticed  until  Edwards  et  al. 

(1979)  repeated  the  experiments,  confirmed  the  results  and  suggested 
that  the  law  should  be  valid  for  other  gas  mixtures.  Subsequently, 
there  have  been  a  number  of  experimental  demonstrations  of  the 
validity  of  this  simple  empirical  law  for  many  different  fuel- 
oxidizer  combinations  at  various  initial  pressures.  Despite  the 
fundamental  significance  and  simplicity  of  this  empirical  law,  no 
theoretical  justification  has  appeared.  However,  extensive 

♦See  Guirao  (1982),  Moen  (1982)  and  Rinnan  (1982). 


calculations  by  Westbrook  and  Urtiew  show  a  strong  correlation 
between  the  critical  tube  diameter  (or  equivalently  the  cell  size) 
and  the  reaction  zone  thickness  as  determined  by  the  detailed 
kinetics  of  the  oxidation  process*. 

To  test  the  range  of  validity  of  this  simple  empirical  law, 
experiments  have  been  carried  out  recently  by  Liu  et  al.  (1983) 
for  detonation  transmission  through  orifices  of  different  geometries 
They  found  that  the  transmission  through  an  orifice  is  identical 
to  that  through  a  tube  of  the  same  cross-sectional  area.  For  the 
two-dimensional  case  of  a  rectangular  orifice  of  large  aspect 
ratio  (i.e.,  L/w  >  5),  Liu  obtained  the  rather  surprising  result 
that  the  critical  channel  width  Wc  for  transmission  is  only 
about  three  cell  diameters.  This  result  is  in  disagreement  with 
that  reported  by  Mitrofanov  and  Soloukhin  (1965),  as  well  as 
Edwards  et  al.  (1979)  who  gave  Wc  =*  10  X  for  the  two-dimensional 
case.  However,  both  Soloukhin' s  and  Edwards'  experiments  were 
based  on  square  channels,  whereas  Liu’s  two-dimensional  orifice 
experiments  covered  a  range  of  aspect  ratios  1  <  L/w  <  15.  For 
L/w  =  1,  Liu's  results  agree  with  Mitrofanov  and  Soloukhin  (1965) 
and  Edwards  et  al.  (1979),  Wc  »  10  X. 

Due  to  the  fundamental  significance  of  this  two-dimensional 

result  reported  by  Liu,  we  felt  it  was  desirable  to  perform  large- 

scale  experiments  to  confirm  that  Wc  =  3  X.  Whereas  Liu’s 

experiments  were  based  on  rectangular  orifice  plates  of  a  maximum 

L/W  =  15,  two-dimensional  channels  with  large  aspect  ratios  (e.g., 

L/w  =  30)  can  readily  be  used  in  field  experiments  to  give  a  truly 

*See  Westbrook,  (1982),  Westbrook  and  Urtiew  (1982),  Westbrook  and 
Urtiew  (1982),  Westbrook  (1980). 


two-dimensional  geometry.  The  present  paper  reports  the  results 
of  a  series  of  such  field  tests  recently  carried  out  at  Sandia 
National  Laboratories  (New  Mexico)  to  verify  the  empirical  law 
Wc  =  3X. 

EXPERIMENTAL  METHOD 

The  experimental  detonation  channel  consisted  of  two  25  mm 
thick  steel  plates  1.83  m  wide  and  2.44  m  long.  The  channel  width 
W  (the  spacing  between  the  plates)  was  varied  by  using  four  spacer 
blocks  placed  at  each  of  the  four  corners  of  the  plates  (Fig.  la). 
One  of  the  side  walls  of  the  channel  is  covered  by  a  sheet  of 
clear  plastic  to  permit  photographic  observations  of  the  detonation 
propagation  in  the  channel,  while  the  opposite  side  wall  of  the 
channel  is  covered  by  a  plywood  sheet.  The  closed  end  of  the  chan¬ 
nel  was  also  covered  by  a  sheet  of  plywood  onto  which  are  taped 
thin  strips  of  sheet  explosive  for  initiating  the  gas  detonation. 

The  quantity  of  strip  explosive  ranged  from  50  to  140g  depending 
upon  the  width  of  the  channel.  The  high  explosive  strips  are 
centrally  initiated  by  an  exploding  bridgewire  detonator.  The 
plastic  and  plywood  side  walls  and  the  initiation  end  of  the  channel 
are  blown  off  during  the  experiment  and  have  to  be  replaced  for 
each  test.  The  open  end  of  the  channel  is  connected  to  a  large 
plastic  bag  fabricated  from  100  ym  thick  polyethylene  sheets  which 
are  taped  to  the  steel  plates.  The  plastic  bag  was  about  3.5  m  in 
length.  The  filling  procedure  is  to  first  retract  the  entire 
plastic  bag  into  the  channel  by  removing  almost  all  of  the  air  in 
the  system  with  a  vacuum  pump;  then  a  premixed  hydrogen:  air  mixture 
is  introduced  from  a  large  underground  storage  tank.  The  plastic 
bag  is  extruded  from  the  channel  as  the  system  is  being  filled 


with  the  explosive  gas  mixture.  The  final  pressure  in  the  system 
is  slightly  above  ambient  so  that  the  plastic  bag  remains  taut. 

On  windy  days  or  when  the  hydrogen  concentration  is  rich  enough 
such  that  the  buoyancy  of  the  bag  is  significant,  thin  nylon  ropes 
are  used  to  anchor  the  plastic  bag  in  the  horizontal  position. 

The  detonation  channel  is  located  on  a  thick  concrete  pad  with  the 
plates  extended  0.55  m  over  the  edge  of  the  slab  to  allow  room  for 
the  inflated  plastic  bag.  The  only  diagnostic  was  high-speed 
cinematography.  Detonation  velocity  as  well  as  successful  transmis 
sion  or  failure  can  be  discerned  quite  readily  from  the  movie 
record.  Hydrogen-air  and  ethylene-air  mixtures  have  been  used  and 
the  desired  composition  was  established  in  the  mixing  tank  by  tne 
method  of  partial  pressures.  To  measure  the  detonation  cell  size, 
soot  was  deposited  on  a  thin  aluminum  sheet  by  a  welding  torch 
burning  rich  mixtures  of  MAPP  (methyl  acetylene,  propadiene,  pro¬ 
pane,  butane  and  propylene)  and  oxygen.  The  smoked  aluminum  shee 
was  then  taped  to  the  bottom  steel  plate  of  the  channel.  We  found 
that  the  smoked  sheet  remains  in  place  during  the  test  and  the  cell 
size  could  be  determined  from  it  after  the  top  plate  of  the  channel 
was  removed  following  the  detonation. 

RESULTS 

For  the  case  of  hydrogen-air  mixtures,  the  cell  size  \  as 
determined  from  the  smoked  aluminum  sheets  placed  in  the  two- 
dimensional  detonation  channel  is  plotted  against  hydrogen-concen¬ 
tration  in  Fig.  2.  Also  shown  for  comparison  are  the  results  for 
\  obtained  at  McGill  by  Knystautas  et  al.  (1982)  a  circular  tube 


of  300  mm  diameter.  For  the  same  hydrogen-air  mixture,  the  present 
data  for  X  are  generally  slightly  larger  than  the  corresponding 
values  obtained  at  McGill  (solid  line).  The  slight  discrepancy 
may  be  accounted  for  by  the  difference  in  the  ambient  pressures 
between  Albuquerque  (~  1600  m  altitude)  and  Montreal  (~  100 
m) .  Since  the  cell  size  is  directly  proportional  to  the  induction 
time  which  is  inversely  proportional  to  the  initial  pressure,  the 
McGill  results  can  be  scaled  to  the  lower  ambient  pressure  of 
Albuquerque,  i.e., 

^Montreal 

^Albuquerque  =  p  ^Montreal 

Albuquerque 

Using  this  correction,  there  is  better  agreement  between  measured 
cell  sizes  and  the  rescaled  McGill  data  (shown  as  the  dotted  curve). 

Rich  mixtures  of  hydrogen-air  were  not  used  in  the  present  series 
of  experiments  since  the  objective  was  only  to  verify  the  empirical 
law  Wc  =  3  X  for  two-dimensional  channels.  Extensive  experiments 
have  already  been  carried  out  for  the  rich  hydrogen-air  mixtures  for 
circular  tubes  and  have  been  reported  previously  Guirao  (1982). 

It  should  be  noted  that  the  cell  size  is  indicative  of  the  rate 
processes  in  the  complex  detonation  structure  and  is  a  very  useful 
parameter.  Detonation  velocities  deduced  from  the  high  speed 
movies  were  in  agreement  with  the  Chapman-Jouguet  values  within 
our  estimated  experimental  error. 

The  normalized  critical  channel  height  Wc/ X  measured  in 
the  present  experiments  is  plotted  against  the  aspect  ratio  of  the 
channel  L/W  in  Fig.  3.  Note  that  for  large  aspect  ratios, 

L/w  >  10,  the  present  results  give  a  value  Wc/ X  =  3  in  approximate 
agreement  with  the  results  obtained  by  Liu  using  rectangular 


orifices  instead  of  a  full  two-dimensional  channel.  Experiments 
with  aspect  ratios  up  to  35  yield  the  same  value  for  Wc  2  3x. 

Using  different  hydrogen  concentrations  and  hence  different  X's 
and  different  Chapman- Jouguet  states,  also  gives  the  same  results, 

Wc/ X  =  3.  This  observation  indicates  that  the  critical  channel 
height  Wc  scales  according  to  the  cell  size  X,  and  within  the 
present  experimental  accuracy  does  not  depend  explicitly  on  the 
Chapman-Jouguet  states.  Several  experiments  were  also  carried 
out  for  6%  ethylene-air  mixtures  with  aspect  ratio  of  14.5  -  18.5. 

The  result  Wc/x  =  3  was  obtained  for  ethylene  as  well,  reinforcing 
our  conjecture  that  the  scaling  law  for  channels  is  based  only  on 
the  cell  size  X  as  has  been  demonstrated  in  the  case  of  circular 
tubes  (Knystautas  et  al.  1982). 

In  the  present  two-dimensional  geometry,  it  was  also  possible 
to  test  the  symmetry  of  the  transmission  process  about  the  centerline 
of  the  channel.  By  retracting  the  top  plate  of  the  channel  relative 
to  the  bottom  plate  so  that  the  emergent  detonation  is  supported  by 
one  plane,  we  constructed  a  half  channel  (Fig.  lb).  By  symmetry, 

transmission  should  now  be  possible  for  a  spacing  between  the  plates 
corresponding  to  only  1.5  x.  Several  experiments  were  conducted 
using  20  -  20.5%  hydrogen-air  mixtures  with  channel  widths  from  10 
to  10.5  cm  and  a  channel  length  of  1.82  m,  yielding  17.5  <  L/W  <,  18.2. 
The  go  -  no  go  results  bracketed  a  critical  spacing  between  the 
plates  of  close  to  1.5  X,  confirming  that  the  two-dimensional 
transmission  mechanism  is  symmetrical  about  the  channel  centerline. 

For  aspect  ratios  L/W  <  8,  the  present  results  indicate  that 
the  value  for  Wc/X  rises  sharply.  A  similar  trend  was  also  obtained 
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toy  Liu  for  the  case  of  transmission  through  rectangular  orifices, 
but  the  increase  does  not  occur  until  L/W  <  5.  Why  three-dimensional 
effects  begin  to  become  important  at  larger  aspect  ratios  for  the 
channel  than  for  the  orifice  plates  has  not  been  resolved.  Other 
than  this  discrepancy,  the  present  channel  experiments  are  in 
qualitative  agreement  with  the  rectangular  orifice  plate  experiments 
of  Liu. 

To  explain  the  surprising  result  that  the  critical  condition 
for  the  successful  transmission  of  detonation  changes  from  Wc  =  10X. 
for  a  square  channel  to  Wc  =  3\  for  a  rectangular  channel  of  large 
aspect  ratio  (L/w  >  8),  one  can  invoke  the  wave  curvature  concept 
proposed  recently  by  Lee  et  al.  (1984).  Lee’s  arguments  suggest  that 
when  a  planar  detonation  wave  propagating  in  a  confining  linear 
tube  or  channel  emerges  suddenly  into  an  unconfined  geometry,  it 
is  subject  to  a  gas  dynamic  expansion  process  by  virtue  of  the 
rarefaction  waves  that  are  generated  at  the  edge  which  then  pene¬ 
trate  into  the  wave  front  and  impose  a  curvature  on  the  wave.  He 
then  defines  the  critical  condition  for  transmission  in  terms  of 
the  minimum  radius  of  curvature  of  the  diffracted  wave,  namely; 
if  the  rarefaction  waves  give  rise  to  a  curvature  of  the  diffracted 
detonation  exceeding  a  certain  critical  value,  then  failure  of  the 
wave  results.  For  a  three-dimensional  diffraction  process  the 
minimum  radius  of  curvature  would  be  of  the  order  of  the  hydro- 
dynamic  thickness  A^  so  that  R  =  A^  =  6.5X..  Thus  for  the  dif¬ 
fraction  of  detonation  through  a  square  opening,  the  expansion 
process  is  three-dimensional  and  the  curvature  imposed  on  the  dif¬ 
fracting  wave  would  be  spherical  in  nature  corresponding  to  some 


effective  radius  of  curvature  related  to  the  average  of  the  mini- 
mun  and  maximum  dimensions  of  the  square  opening  R<*(W+/2W)/2 . 

Liu  et  al.  (1983)  has  shown  that  for  this  case  Wc  »  10\  in  agree¬ 
ment  with  the  experimental  results  of  Soloukhin  (1965)  and  Edwards 
(1979)  . 

For  rectangular  channels  of  high  aspect  ratio  (L/W>8);  the 
planar  wave  diffracting  through  the  opening  transforms  into  a 
cylindrical  wave.  For  the  same  critical  curvature,  by  simple 
geometric  considerations,  namely  that  xCyi  =  xSph  =  (l/A)dA/dR 
one  can  show  that  the  radius  of  curvature  for  the  cylindrical  wave 
would  be  one  half  that  for  the  spherical  case  so  that  RCyq  = 
RSph/2=AH/2  =  3.25X  which  is  in  agreement  with  the  present  experi¬ 
mental  results. 

CONCLUSION 

The  present  series  of  large  scale  tests  confirm  the  smaller 
laboratory  experiments  of  Liu  et  al.  that  transmission  from  confined 
planar  detonations  to  unconfined  cylindrical  detonations  requires 
a  critical  channel  height  of  only  about  three  cell  diameters. 
Furthermore,  the  channel  results  are  almost  identical  to  the  results 
from  a  rectangular  orifice.  This  suggests  that  the  re-initiation 
mechanism  is  very  local  and  does  not  depend  strongly  on  the  flow 
structure  of  the  product  gas  emerging  into  confined  space  behind 
the  leading  front.  The  onset  of  three-dimensional  effects  at 
surprisingly  large  aspect  ratios  (i.e.,  L/W  =  5)  suggests  that 
quenching  by  the  rarefaction  waves  may  not  be  a  dominant  mechanism. 
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FIGURE  CAPTIONS 

Fig.  1  Two-dimensional  channel  apparatus  with: 

(a)  unconfined  exit  plane 

(b )  one  degree  of  confinement  at  the  exit  plane. 

Fig.  2  Detonation  cell  size  measured  in  Montreal, 

normalized  for  lower  ambient  pressure  in 
Albuquerque,  and  as  measured  in  channel 
experiments  in  Albuquerque,  all  as  a  function 
of  H2:Air  mixture. 

Fig.  3  Critical  channel  width  normalized  to  detonation 

cell  size,  Wc/x,  as  a  function  of  the  channel  as¬ 
pect  ratio,  L/W,  for  detonations  emerging  from  the 
exit  plane  of  the  channel  into  an  unconfined 
volume.  The  data  point  3\/2  is  for  detonation 
supported  by  one  plane  after  exiting  the  channel. 
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DIFFRACTION  OF  DETONATION  FROM  TUBES  INTO  A  LARGE  FUEL- 

AIR  EXPLOSIVE  CLOUD 

I  O.  MOEN  and  S  B.  MURRAY 
Defence  Research  Establishment  Suf field 
Ralston,  Alberta,  Canada 

D.  BJERKETVEDT  and  A.  RINNAN 
Norwegian  Defence  Construction  Services 
Akershtis,  Oslo,  Norway 

AND 

R.  KNYSTAUTAS  and  J.  H.  LEE 
McCill  University 
Montreal,  Canada 


This  paper  reports  on  a  series  of  field  tests  performed  at  the  Canadian  Defence  Research 
Establishment  Suffield  and  at  Raufoss,  Norway,  to  obtain  the  critical  tube  diameter,  dc,  for 
ethylene-air  mixtures  by  investigating  the  diffraction  of  detonations  from  tubes  into  large 
plastic  bags  simulating  an  unconfined  fuel-air  cloud.  The  critical  ethylene-air  compositions 
for  successful  re-establishment  of  detonation  upon  emerging  from  tubes  of  diameters  0.31 
m,  0.45  m,  0.89  m  and  1.36  m  were  determined  by  monitoring  the  diffracted  detonation 
wave  in  the  bag.  High-speed  cinematography  of  the  diffracted  wave  shows  that  re-establish- 
ment  of  detonation  is  via  one  or  more  re-ignition  centres  at  sites  along  the  head  of  the 
expansion  wave  which  originates  at  the  area  change.  The  characteristic  transverse  wave  spae- 
ings,  S,  associated  with  the  detonations  were  measured  from  smoked  foils  mounted  in  the 
tubes  These  measurements  demonstrate  that  the  empirical  relation  dr  =  13  S  provides  a 
good  correlation  between  the  critical  tube  diameter  and  the  cell  size  over  a  wide  range  of 
ethylene-air  compositions.  However,  a  better  understanding  of  the  coupled  chemical-gas- 
dynamic  processes  within  the  cell  is  required  in  order  to  clarify  the  link  between  cell  size 
and  chemical  kinetics. 

In  addition  to  the  open  tube  tests,  transmission  of  detonations  through  circular  orifice 
holes  was  investigated  in  both  ethylene-air  and  acetvlene-air  mixtures.  Based  on  these  tests 
it  was  concluded  that  the  critical  tube  and  critical  orifice  diameters  are  equal. 

1.  Introduction  The  mechanism  by  which  the  transverse  waves 

are  excited  and  maintained  is  not  completely 
Two  key  properties  which  characterize  the  deto-  understood,  However,  the  recent  detailed  obser- 

nability  of  an  explosive  mixture  are  the  critical  en-  vations  on  the  gasdvnamics  and  chemistry  within 

ergy,  Ef,  required  to  initiate  detonation  in  the  a  detonation  cell  by  Libouton  et  al.2  have  provided 

mixture  and  the  critical  tube  diameter,  dc,  re-  new  insight  into  the  rc-initiation  process  at  the  end 

quired  tor  a  planar  detonation  emerging  from  a  of  the  cell.  It  is  this  re-initiation  which  ensures  the 

tube  to  transmit  to  an  unconfined  spherical  deto-  maintenance  of  the  cell  structure  characteristic  of 

nation  Matsui  and  Lee1  have  proposed  that  cither  a  self-sustained  detonation.  The  delay  to  re-initia- 

Er  or  d,  can  be  used  to  assess  the  relative  sensi-  tion  and  thus  the  length  of  a  detonation  "ell  is  cun- 

tivity  of  explosive  mixtures  to  detonation.  In  fact,  trolled  by  the  chemical  kinetic  properties  of  the 

they  propose  a  simple  relation  between  these  two  explosive  mixture.  The  critical  tube  dian,  'or,  on 

parameters  based  on  a  work-done  model.  On  a  the  other  hand,  is  a  measure  of  minimum  deto- 

more  fundamental  level,  these  properties  arc  re-  nation  size  required  for  the  regeneration  or  cxci- 

lated  to  the  detailed  coupling  between  chemical  tation  of  transverse  waves  as  the  planar  detonation 

energy  release  and  gasdvnamics  which  is  also  re-  wave  emerging  from  the  tube  is  being  quenched 

sponsible  for  the  three-dimensional  transverse  wave  by  the  transversely  advancing  expansion  waves, 

or  cellular  structure  of  detonation  waves.  The  first  observation  that  a  minimum  tube  di- 
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amcter  is  required  for  an  established  detonation  to 
emerge  from  a  tube  and  become  a  spherical  deto¬ 
nation  in  an  unconfined  cloud  was  made  by  Zcl- 
do\ ich  et  al.  ’  Mitrofanov  and  Soloukhin4  later  ob- 
served  that  the  critical  diameter  for  oxy-aeetylenc 
systems  is  10  to  13  transverse  wave  spacing*.  Ed- 
wards  et  al.’6  have  provided  further  evidence  that 
a  minimum  number  of  transverse  waves  is  required 
for  successful  transmission  They  also  suggest  that 
the  number  of  transverse  waves  required  depends 
only  on  the  geometry  of  the  tube  (i.e.,  dc  =  13  S 
for  a  circular  tube  and  r/(  =  10  S  for  a  planar  chan¬ 
nel  and  not  on  the  detonative  system.  The  relation 
dr  =  1.3  S  for  critical  transmission  from  circular 
tubes  has  recently  been  confirmed  for  a  variety  of 
oxygen-enriched  fuel-air  mixtures  by  Knystautas  et 
ah’  However,  the  direct  tests  of  this  empirical  re¬ 
lation  haxe  been  limited  to  mixtures  having  rela¬ 
tively  small  cell  structure  (S  —  16  mm).  Edwards 
et  al.6  present  some  evidence  that  detonations  in 
systems  with  large  and  less  regular  cell  structure 
require  more  transverse  waves  to  ensure  transmis¬ 
sion  of  the  diffracted  wave.  It  is  therefore  impor¬ 
tant  to  examine  the  relation  between  cell  size  and 
criti  il  tube  diameter  for  larger  diameter  tubes. 

In  the  present  tests,  the  critical  C2H4-air  com¬ 
positions  for  successful  re-establishment  of  deto¬ 
nation  upon  emerging  from  tubes  of  diameters  0.31 
m,  0.35  m,  0.89  m  and  1.36  m  were  determined 
by  monitoring  the  diffracted  wave  in  the  bag  using 
ionization-gap  probes  and  pressure  transducers.  In 
addition,  high-speed  cinematographic  records  of  the 


diffracting  detonation  wave  were  obtained  using 
two  high-speed  cameras.  The  detonation  cell  struc¬ 
ture  was  monitored  by  the  smoked-foil  technique. 


2.  Experimental  Details 

The  experimental  field  facility  at  DRES  has  been 
described  in  detail  elsewhere.8,9  The  facility  is 
centred  around  an  18.3  in  x  7.6  m  concrete  test 
pad  onto  which  the  experimental  apparatus  can  be 
mounted.  A  photograph  of  the  test  pad  with  a  typ¬ 
ical  test  section  used  in  the  critical  tube  diameter 
tests  is  shown  in  Fig.  1.  The  test  sections  for  these 
tests  consisted  of  steel  tubes  connected  to  large 
plastic  bags  constructed  from  0. 13  mm  thick,  seam¬ 
less,  extruded  polyctnylene  tubing.  The  dimensions 
of  tube-bag  configurations  used  in  the  present  tests 
are  also  shown  in  Fig.  1. 

At  DRES,  the  test  gas  (CP  grade,  99.5%  pure 
ethylene)  was  mixed  with  the  initial  air  in  the  test 
volume  by  a  multipath  re-circulation  system  using 
a  high-capacity  centrifugal  blower.  The  composition 
and  mixture  homogeneity  in  the  test  volume  were 
verified  by  continuously  analyzing  samples  taken 
from  four  ports  positioned  strategically  along  the 
apparatus  using  a  calibrated  “Wilks  Miran  80”  in¬ 
frared  gas  analyzer.  This  sampling  system  is  ade¬ 
quate  to  guarantee  the  concentration  and  homo¬ 
geneity  of  the  mixture  to  within  ±0.05%  C2H4. 
Once  the  desired  fuel-air  composition  had  been 
attained,  detonation  was  initiated  by  using  a  De- 
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0.45 
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Fu.  1  Photograph  showing  the  DRES  testing  facility  with  a  detonation  tube  (0.45  m  diameter)  and 
polyethylene  hag  l  75  m  diameter)  test  section.  Dimensions  for  each  of  the  four  critical  tube  configu¬ 
rations  investigated  appear  in  the  table. 
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tasheet  initiator  charge  mounted  on  the  plywood 
sheet  which  covered  the  ignition  end  of  the  tube. 

Four  types  of  diagnostics  were  employed  in  the 
DRES  tests.  Eight  piezo-electric  pressure  trans¬ 
ducers  (PCB  113A24)  were  positioned  along  the 
apparatus.  Pressure-time  signatures  and  velocities 
obtained  from  the  pressure  records  were  used  to 
determine  success  or  failure  of  transmission  from 
the  tube  to  the  bag.  In  addition,  a  20-channel  ion¬ 
ization-gap  probe  system,  together  with  a  time-of- 
arrival  counter,  was  used  to  independently  measure 
wave  velocities.  Cinematographic  records  were  ob¬ 
tained  by  using  a  “Hycam"  camera  (15,000-18,000 
halt  frames  per  second)  looking  normal  to  the  di¬ 
rection  of  propagation  and  a  "Fastax"  camera 
i  -5,000  full  frames  per  second!  looking  at  a  45° 
angle  to  the  direction  of  propagation.  Finally,  350 
•mil  x  500  mm  smoked  foils  mounted  in  the  tube 
near  the  exit  to  the  bag  were  used  to  record  the 
detonation  wave  structure.  Each  foil  was  made  by 
applying  a  thin  layer  of  carbon  soot  to  a  polished 
tin  sheet. 

The  experimental  ce  iguration,  diagnostics  and 
procedure  used  in  tl  Raufoss  tests  were  similar. 
In  addition  to  open  tube  tests,  transmission  through 
circular  orifice  holes  was  investigated.  Each  orifice 
hole  was  cut  in  a  20  mm  thick  plywood  disc  mounted 
oxer  the  exit  end  of  the  tube.  These  transmission 
tests  were  performed  using  both  ethylene-air  (C2H4- 
air'  and  acetylene-air  (C2H2-air)  mixtures. 

The  gas  flow  system  used  at  Raufoss  was  based 
on  the  flow  displacement  method.  Five  test  vol¬ 
umes  were  displaced  to  guarantee  the  desired  mix¬ 
ture  composition.  Infrared  analysis  was  employed 
as  in  the  DRES  tests.  However,  only  the  input  gas 
mixture  was  sampled.  The  estimated  uncertainty  in 
mixture  composition  in  the  test  section  was  ±0. 1% 
fuel.  Pressure  transducers,  ionization-gap  probes 
and  high-speed  cinematography  were  the  principal 
diagnostics.  Smoked  foils  were  not  used. 


3.  Results  and  Discussion 

i'  Critical  Tube  Diameter  and  Transverse  Waves 

Selected  frames  from  a  high-speed  (—15,000  frames/ 
sec)  cinematographic  record  showing  the  successful 
transmission  of  detonation  from  a  0.89  m  diameter 
tube  into  a  large  1.75  m  diameter  bag  simulating 
an  unconfined  gas  cloud  are  shown  in  Fig.  2.  Frame 
1  shows  the  planar  detonation  wave  just  as  it 
emerges  from  the  tube  Frame  2  shows  the  shrink¬ 
ing  planar  detonation  core  as  the  expansion  waves 
which  originate  at  the  tube  exit  sweep  in  towards 
the  centre  quenching  the  detonation  in  their  wake. 
The  head  of  the  expansion  fan  reaches  the  centre 
at  a  distance  of  approximately  one  tube  diameter 
from  the  tube  exit  Unless  re-ignition,  capable  of 


Fig.  2.  Selected  frames  from  a  high-speed  cine¬ 
matographic  record  (—15,000  frames/sec)  showing 
successful  transmission  of  detonation  from  a  con¬ 
fined  tube  of  0.89  m  diameter  to  a  large  gas  bag 
of  1.75  m  diameter  simulating  an  unconfined  en¬ 
vironment  (0.2  msec  between  frames). 


re-establishing  the  detonation  in  the  adverse  ex¬ 
pansion,  occurs  somewhere  behind  the  leading  shock 
wave  prior  to  this,  the  detonation  fails  to  re-estab¬ 
lish  itself  and  the  explosive  mixture  in  the  bag  is 
consumed  by  a  deflagration.  For  the  test  shown  in 
Fig.  2,  re-ignition  is  seen  to  occur  in  Frame  3  at 
an  ignition  nucleus  in  the  centre  of  the  bag.  The 
subsequent  formation  of  detonation  bubbles  which 
sweep  across  the  entire  cross  section  of  the  bag  are 
clearly  seen  in  Frames  4-6  of  the  sequence.  The 
final  two  frames  of  the  sequence  show  a  fully  es¬ 
tablished  detonation  in  the  bag. 

Fig,  2  represents  a  typical  sequence  showing 
close  to  critical  re-establishment  of  detonation  after 
diffraction  from  a  tube.  If  the  sensitivity  of  the 
mixture  is  reduced  for  the  same  diameter  tube,  no 
re-ignition  centres  are  observed  and  the  detonation 
fails  to  re-establish  itself  in  the  bag,  while  with 
more  sensitive  mixtures  re-ignition  occurs  closer  to 
the  tube  exit,  usually  at  more  than  one  ignition 
centre  at  the  edge  of  the  shrinking  detonation  core. 
Only  for  mixtures  much  more  sensitive  than  critical 
composition  does  the  detonation  wave  transmit  di¬ 
rectly  from  a  planar  wave  to  a  hemispherical  wave 
without  the  appearance  of  re-ignition  centres. 

The  above  observations  are  in  agreement  with 
the  laboratory  Schlieren  and  smoked-foil  observa¬ 
tions  on  the  diffraction  of  detonations  in  low  pres¬ 
sure  oxy-acetylene  reported  by  Edwards  et  al.5,6 
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Their  records  also  show  that  the  re-establishment 
of  detonation  after  an  abrupt  area  expansion  is  via 
the  formation  of  re-ignition  centres  at  sites  along 
the  head  of  the  expansion  which  originates  at  the 
area  change.  In  spite  of  the  apparent  randomness 
in  the  appearance  of  these  re-ignition  centres,  the 
critical  mixture  composition  for  re-establishment  of 
detonation  for  a  given  tube  diameter  is  well  defined 
provided;  i)  the  detonation  emerging  from  the  tube 
is  a  planar  detonation  wave  no  longer  influenced 
by  the  blast  wave  from  the  initiating  charge,  and 
ii)  the  plastic  bag  is  sufficiently  large  so  that  it  does 
not  influence  the  transmission  process. 

In  the  present  tests,  care  was  taken  to  ensure 
that  a  steady-state  planar  detonation  had  been  es¬ 
tablished  and  that  the  initiating  charge  did  not  in¬ 
fluence  the  critical  fuel-air  composition.  This  was 
accomplished  by  varying  the  initiation  energy  once 
the  critical  composition  had  been  identified  and  by 
carefully  monitoring  both  the  pressure  and  velocity 
of  the  detonation  wave  in  the  tube. 

In  order  to  ensure  that  wave  interactions  with 
the  plastic  bag  did  not  influence  the  transmission, 
tests  were  done  using  plastic  bags  of  different  sizes. 
The  re-establishment  of  detonation  due  to  wave  in¬ 
teractions  with  the  plastic  bag  was  observed  on  the 
high-speed  film  records  from  only  one  test  (0.45  m 
tube  into  a  0.89  m  plastic  bag  at  6.3%  C2H4).  The 
apparent  Go  outcome  of  till  test  has  therefore 
been  discounted. 

The  results  from  the  critical  tube  diameter  tests 
are  summarized  in  Fig.  3,  where  the  Go  -  No-Go 
outcome  of  tests  with  different  diameter  tubes  is 
shown  as  a  function  of  %C2H4  by  volume.  Also 
shown  in  this  figure  are  the  results  from  the  1980 
Raufoss  tests.10  Except  for  the  0.45  m  tube,  the 
results  from  the  two  test  series  are  in  good  agree¬ 
ment  considering  the  different  experimental  pro¬ 
cedures  and  conditions  at  DRES  and  Raufoss.  The 
Raufoss  results  indicate  successful  re-establishment 
of  detonation  in  leaner  mixtures  for  a  given  tube 
diameter.  This  could  be  due  to  the  fact  that  air 
containing  impurities  in  the  form  of  oil  mist  from 
an  air  compressor  was  used  in  the  Raufoss  tests. 
These  impurities  could  have  an  effect  on  the  mix¬ 
ture  sensitivity  and  also  on  the  monitoring  of  fuel 
concentration.  At  DRES,  atmospheric  air  was  used 
in  all  the  tests. 

It  is  important  to  note  that  no  detonation  re-cs- 
tablishment  was  observed  with  a  0.31  m  tube  up 
to  an  ethylene  concentration  of  10.5%.  Since  it  is 
very  unlikely  that  re-cstablishment  will  occur  in 
even  richer  mixtures,  we  can  conclude  that  the 
critical  tube  diameter  is  larger  than  0.31  m  for  any 
C2H4-air  mixture. 

The  solid  and  chain-dotted  curves  shown  in  Fig. 
3  are  theoretical  predictions  based  on  a  correlation 
proposed  by  Moen  et  al.u  for  initial  pressures  of 
101.3  kPa  and  92.5  kPa,  respectively;  the  latter 


pressure  being  the  mean  atmospheric  pressure  on 
the  DRES  range  where  most  of  the  tests  were 
done.  The  dashed  curve  is  a  similar  theoretical  pre¬ 
diction  using  a  more  detailed  chemical  reaction 
scheme  proposed  by  Westbrook.12,13 

All  of  these  predictions  are  based  on  empirical 
relations  between  the  critical  tube  diameter,  dc, 
the  transverse  wave  spacing,  S,  and  the  induction 
zone  length,  A,  calculated  for  the  hypothetical  one¬ 
dimensional  C-J  detonation.  It  is  assumed  that  the 
detonation  cell  size  is  determined  by  the  chemical 
induction  time,  t,  behind  a  shock  wave  propagating 
at  the  C-J  velocity,  with  S  being  directly  propor¬ 
tional  to  the  induction-zone  length  A  =  t u,  where 
u  is  the  post-shock  relative  particle  velocity.  When 
combined  with  the  empirical  13  S  formula  for  dc, 
this  gives  dc  =  13Atu,  where  A  is  the  proportion¬ 
ality  constant  relating  S  and  A.  Thus,  once  the  C-J 
shock  parameters  and  the  corresponding  induction 
times  are  known  as  functions  of  mixture  composi¬ 
tion,  the  critical  tube  diameters  are  determined  by 
normalizing  at  one  composition.  For  the  ethylene 
system  Moen  et  al.11  use  the  induction-time  for¬ 
mula  of  Hidaka  et  al.14 

log10T[O2]  =  -11.45  +  27.5  x  103/4.58T, 

where  t  is  the  induction  time  in  seconds  and  [02] 
is  the  oxygen  concentration  in  moles/liter.  They 
then  obtain  A  =  480  by  monitoring  the  onset  of 
near-limit  phenomena  in  circular  tubes  and  find 
that  the  predicted  critical  tube  diameters  for 
C2H4 — 02 — N2  mixtures  having  different  nitrogen 
to  oxygen  ratios  (3  (O  £  p  s  3)  are  in  good  agree¬ 
ment  with  experimental  results.  As  can  be  seen  in 
Fig.  3,  the  predicted  dc  for  C2H4-air  mixtures,  us¬ 
ing  the  same  value  for  the  constant  of  proportion¬ 
ality,  are  also  in  reasonable  agreement  with  ex¬ 
perimental  results  over  a  fairly  wide  range  of  C2H4 
concentrations.  However,  the  predicted  dc  based 
on  Westbrook’s  model  is  out  by  more  than  a  factor 
of  two  at  4.5%  C2H4  when  normalized  to  dc  =  0.36 
m  at  stoichiometric  composition  (6.54%  C2H4).  Based 
on  C2H4 — 02 — N2  laboratory  results11  Westbrook’s 
model  predicts  a  dc  larger  than  0.36  m  at  stoichi¬ 
ometric  composition,  so  that  agreement  with  all 
available  experimental  data  on  the  ethylene  system 
is  not  improved  by  re-normalizing  to  obtain  better 
agreement  with  the  present  lean  C2H4-air  results. 

There  are  two  reasons  for  the  difference  between 
these  two  model  predictions.  Westbrook’s  results 
are  obtained  using  a  detailed  chemical  reaction 
scheme  for  ethylene  oxidation  which  appears  to 
predict  an  effective  activation  energy  much  larger 
(~48  kcal)  than  the  27.5  kcal  obtained  by  Hidaka 
et  al. 14  from  shock  tube  studies  on  C2I14— 02  mix¬ 
tures  highly  diluted  by  argon.  Furthermore,  West¬ 
brook  uses  the  equilibrium  temperature  behind  the 
shock  wave  to  calculate  the  induction  time.  This 
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with  experiment  over  a  wider  range  of  mixture  sen¬ 
sitivity.  In  view  of  the  rather  arbitrary  choice  of 
C-J  conditions  for  determining  the  cell  size,  this 
is  not  surprising.  It  is  well  established  that  the 
shock  velocity  within  a  detonation  cell  varies  from 
an  initial  value  of  about  1.6  l'Cj  to  a  minimum  of 
about  0.8  U(  i  near  the  end  of  the  cell.  Although 
C-J  conditions  represent  the  average  conditions  in 
the  cell,  the  re-initiation  process  within  the  cell 
observed  by  Libouton  et  al.2  does  not  occur  at 
these  average  conditions.  A  better  understanding 
of  the  coupled  chemical  gasdynamic  processes  within 
the  cell  is  clearly  required  in  order  to  properly 
model  the  link  between  cell  size  and  chemical  ki¬ 
netics. 

The  link  between  cell  size  and  critical  tube  di¬ 
ameter  is  also  not  understood.  However,  in  this 
case  the  13  S  relation  continues  to  hold  for  less 
sensitive  mixtures.  The  transverse  wave  spacings 
measured  in  the  tubes  are  plotted  as  a  function  of 
%C2H4  by  volume  in  Fig.  4.  The  results  are  in 
good  agreement  with  the  measurements  of  Bull  et 
al.,16  Borisov17  and  Knystautas  et  al.7  at  stoichio¬ 
metric  composition,  and  also  with  the  results  of 
Moen  et  al.11  obtained  by  monitoring  the  onset  of 
near-limit  phenomena  in  tubes.  The  curve  in  this 
figure  corresponds  to  the  solid  curve  in  Fig.  3  with 


(,'.‘Nr.lNr»>ATiON  Of  m4  IN  AIM  :oi 


Fit..  4.  Transverse  wave  spacing  versus  concen¬ 
tration  iTCTI,  in  C .'  1 1 4-air  by  volume).  Bull  et 
al  — Kef  !6.  Borisov — Ref.  17;  Moen  et  al. — Ref 
11.  Knvstautas  et  al.  —  Ret.  7. 


S  =  dc/ 13.  Although  there  is  considerable  scatter 
in  the  data,  the  empirical  13  S  formula  provides  a 
reasonably  good  correlation  over  a  wide  range  of 
C2H4-air  compositions.  The  scatter  in  the  data  is 
partly  a  reflection  of  the  complex  structure  ob¬ 
served  on  many  of  the  smoked-foil  records.  The 
structure  typically  consists  of  large  cells  superim¬ 
posed  on  a  smaller  much  more  irregular  structure, 
so  that  the  identification  and  measurement  of  the 
relevant  cell  structure  can  be  difficult. 

A  typical  smoked-foil  record,  together  with  a 
sketch  illustrating  our  interpretation  of  the  record, 
is  shown  in  Fig.  5.  Two  independent  measuring 
techniques  were  used  to  determine  the  most  rep¬ 
resentative  transverse  wave  spacings  from  these 
records.  The  first  method  was  simply  to  take  the 
average  value  of  the  cell  width  from  a  large  num¬ 
ber  of  cells.  In  this  approach  some  judgement  as 
to  what  constitutes  a  cell  is  required.  This  is  avoided 
to  some  extent  in  the  second  method  which  is 
based  on  identification  of  dominant  diagonal  bands 
of  similarly  shaped  cells.  For  low  mode  number 
detonations  in  the  tube  the  characteristic  cell  size 
can  also  be  deduced  from  the  frequency  of  the 
pressure  oscillations  behind  the  leading  detonation 
wave. 

An  interesting  observation  is  that  the  width  of 
the  longitudinal  plastic  strips  produced  by  a  deto¬ 
nation  in  the  bag  appears  to  be  equal  to  the  trans¬ 
verse  wave  spacing.  The  widths  of  the  plastic  strips 
observed  for  different  C2H4-air  compositions  are 
included  in  Fig.  4.  Many  of  the  plastic  strip  results 
were  obtained  in  a  subsequent  test  series  involving 
the  transmission  from  a  tube  into  a  plastic  bag  of 
the  same  diameter.18 

ii)  Transmission  Through  Orifices 

As  long  as  the  re-establishment  of  detonation 
upon  diffraction  is  localized  (i.e.,  confined  to  a  re¬ 
gion  in  the  immediate  neighbourhood  of  the  dif¬ 
fracted  detonation  core  where  the  complex  shock 
reflections  from  an  orifice  plate  would  not  influ¬ 
ence  the  gasdynamic  flow  field),  the  critical  con¬ 
ditions  for  transmission  through  a  circular  orifice 
should  be  the  same  as  those  for  a  circular  tube. 
The  equality  of  the  critical  tube  diameter  and  crit¬ 
ical  circular  orifice  diameter  was  first  demonstrated 
in  a  series  of  laboratory  tests  with  oxygen-enriched 
ethylene-air  and  acetylene-air  mixtures  for  orifices 
of  diameters  89,  51,  25  and  13  mm  at  the  end  of 
a  200  mm  diameter  tube  using  the  apparatus  de¬ 
scribed  by  Knystautas  et  al. 1  The  large-scale  tests 
with  orifices  of  diameters  0.96  m,  0.44  in  and  0.20 
m  at  the  end  of  the  1.36  m  diameter  tube  provide 
evidence  that  this  equality  also  holds  for  less  sen¬ 
sitive  fuel-air  mixtures.  The  Raufoss  1981  experi¬ 
mental  points  at  dr  =  0.96  in  and  0.44  m  in  Fig. 
3  were  obtained  from  these  transmission  tests.  In 
view  of  the  significant  disagreement  between  the 
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present  results  and  the  1980  Raufoss  results10  at  dc 
=  0.45  m  and  the  fact  that  only  a  Go  result  is  avail¬ 
able  for  the  0.96  m  diameter  orifice,  these  results 
are  not  very  convincing  evidence  of  the  equality. 
More  convincing  evidence  was  obtained  for  C2H2- 
air  mixtures  with  orifice  holes  of  diameters  0.44  m 
and  0.20  m.  As  can  be  seen  in  Fig.  6,  these  results 
agree  with  the  1980  critical  tube  results10  within 
the  experimental  uncertainty  in  fuel  concentration 
(±0.1%  C2H2).  Also  shown  in  Fig.  6  are  critical 
tube  results  at  stoichiometric  composition  (7.75% 
C2H2s6)  obtained  by  Knystautas  et  al.,7  and  by  Frei- 
wald  and  Koch,19  together  with  the  prediction  based 
on  13  S  from  Bull  et  al.16 

The  theoretical  curves  shown  in  Fig.  6  were  de¬ 
rived  in  the  same  manner  as  those  for  C2H4-air 
mixtures  in  Fig.  3.  The  solid  curve  was  arrived  at 
by  the  procedure  of  Moen  et  al.11  using  the  in¬ 
duction-time  formula  of  White20  with  an  activation 
energy  of  17.3  kcal,  and  the  dashed  curve  was  ob¬ 
tained  using  Westbrook’s  model  for  C2H2  oxida¬ 
tion.1213  As  was  the  case  for  ethylene-air  mixtures, 
Westbrook’s  model  appears  to  predict  too  strong  a 
variation  of  dc  with  fuel  concentration. 


Fig.  6.  Critical  tube  diameter  versus  concentra¬ 
tion  i'/cC.H.  in  C2H,-air  by  volume).  Raufoss  1980 — 
Ref.  10,  Knystautas  et  al. — Ref.  7;  Freiwald  and 
Koch — Ref.  19;  Bull  et  al. — Ref.  16;  Moen  et  al. — 
Ref.  11,  Westbrook— Refs.  12,  13. 


4.  Conclusion 

Large-scale  field  tests  on  the  diffraction  of  det¬ 
onations  in  ethylene-air  mixtures  have  provided  key 
data  on  the  detonability  of  ethylene-air  for  a  range 
of  fuel  concentrations  from  3.6%  to  10.5%  C2H4  by 
volume.  These  tests  show  that  the  critical  tube  di¬ 
ameter  for  any  ethylene-air  mixture  is  larger  than 
0.31  m.  The  critical  tube  diameter  near  stoichio¬ 
metric  composition  (6.54%  C2H4)  is  0.45  m,  which 
is  slightly  larger  than  the  values  of  0.36  m  to  0.38 
m  obtained  by  extrapolating  from  laboratory  results 
in  oxygen-enriched  ethylene-air  mixutres.7,11  This 
difference  may  be  partly  due  to  the  low  atmo¬ 
spheric  pressure  (92.5  kPa)  on  the  DRES  range 
where  most  of  the  large-scale  tests  were  per¬ 
formed.  In  the  largest  1.36  m  diameter  tube  the 
critical  fuel  concentration  is  found  to  be  4.5%  C2H4. 
According  to  the  work-done  model  of  Matsui  and 
Lee,1  the  critical  initiation  energy  Ec  varies  as  dc3. 
Based  on  this  assumption,  the  present  results  pre¬ 
dict  that  Ec  at  4.5%  C2H4  is  about  30  times  that 
at  stoichiometric  composition.  This  prediction  is  in 
good  agreement  with  large-scale  measurements  of 
the  critical  energies15  and  thus  provides  further 
evidence  that  either  Ec  or  dc  can  be  used  to  de¬ 
termine  the  relative  detonation  hazards  of  explosive 
mixtures,  as  originally  proposed  by  Matsui  and  Lee.1 

Measurements  of  transverse  wave  spacings  or 
detonation  cell  sizes  show  that  the  empirical  rela¬ 
tion  dc  =  13  S,  first  discovered  for  oxy-acetylene4 
and  found  to  be  valid  for  a  variety  of  oxygen-en¬ 
riched  fuel-air  mixtures,7  is  also  valid  for  fuel-air 
mixtures  with  characteristic  transverse  wave  spac¬ 
ings  up  to  about  100  mm.  However,  the  correlation 
between  cell  size  and  chemical  kinetics,  based  on 
a  linear  relationship  between  the  transverse  wave 
spacing  and  the  induction-zone  length  calculated  at 
C-]  conditions  using  the  C2H4  oxidation  scheme  of 
Westbrook,12 13  is  not  in  good  agreement  with  the 
present  data.  The  data  are  in  better  agreement 
with  a  similar  empirical  correlation  proposed  by 
Moen  ct  al.11  which  uses  a  smaller  effective  acti¬ 
vation  energy  and  a  larger  temperature  to  calculate 
the  induction-zone  length. 

High-speed  cinematographic  records  of  the  det¬ 
onation  diffraction  process  show  that  the  re-estab¬ 
lishment  of  detonation,  after  the  abrupt  area  ex¬ 
pansion  from  a  tube  to  a  large  bag,  is  via  the 
formation  of  localized  re-ignition  centres  at  sites 
along  the  head  of  the  expansion  which  originates 
at  the  area  change.  At  criticality,  re-ignition  occurs 
at  a  localized  ignition  centre  near  the  tube  axis, 
approximately  one  tube  diameter  from  the  tube 
exit.  The  localized  nature  of  this  re-ignition,  im¬ 
mediately  behind  the  shrinking  detonation  core, 
suggests  that  the  critical  conditions  for  re-cstablish- 
ment  of  detonation  depend  only  on  the  frontal  area 
of  the  emerging  planar  detonation  wave.  The 
equality  of  the  critical  tube  diameter  and  the  crit- 
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ieal  circular  orifice  hole  diameter  is  further  evi¬ 
dence  for  this  conclusion. 

Further  refinements  in  both  theory  and  experi¬ 
ment  are  necessary  in  order  to  obtain  a  more  fun¬ 
damental  understanding  of  the  detonability  prop¬ 
erties  of  fuel-air  mixtures.  The  present  results, 
together  with  the  critical  energy  results  from  pre¬ 
vious  large-scale  tests'3  provide  key  data  for  fuel- 
air  mixtures  with  characteristic  detonation  length 
scales  (i.e.,  S  and  dc)  that  cannot  be  obtained  in 
the  laboratory.  These  data  are  not  only  of  practical 
interest  for  accidental  explosions  and  military  ap¬ 
plications,  but  also  provide  critical  tests  of  theo¬ 
retical  models  and  correlations  linking  the  detona¬ 
bility  properties  to  the  chemical  kinetics  of  the 
explosive  mixture.  Based  on  detailed  observations 
of  the  detonation  diffraction  process,  the  monitor¬ 
ing  of  detonation  structure  and  comparison  of  re¬ 
sults  with  theoretical  models,  we  can  conclude  that 
a  more  fundamental  understanding  of  the  detona¬ 
bility  properties  can  only  be  achieved  through  fur¬ 
ther  elucidation  on  the  nature  of  the  re-ignition 
process,  both  in  the  diffracted  detonation  wave  and 
at  the  end  of  a  detonation  cell. 
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COMMENTS 


C.  Westbrook,  Lawrence  Livermore  Lab., 
USA.  This  paper  is  a  very  fine  collection  and  pre¬ 
sentation  of  experimental  data  which  are  of  unusual 
value  in  developing  new  models  for  detonation 
properties  and  understanding  such  phenomena. 

However,  the  use  of  simplified,  particularly  one- 
step,  kinetics  expressions  for  induction  times  in 
detonations  cannot  really  be  justified  any  longer  for 
fuels  where  detailed  models  exist,  such  as  C,H4  and 
CTL  The  Hidaka  expression  used  here  was  de¬ 
rived  for  highly  dilute  mixtures  in  shock  tubes  and 
covered  only  a  very  limited  range  of  fuel-oxygen 
ratio.  Its  agreeement  with  these  experimental  det¬ 
onation  data  can  only  be  fortuitous. 

Both  cell  size  and  critical  tube  diameter  for  C2ll(- 
air  must  increase  for  fuel  concentrations  in  excess 
of  7-8%  C,H4  as  predicted  by  my  own  kinetics 
model  as  referenced  in  this  paper.  However,  using 
the  simplified  expression,  the  curve  of  dr  vs  fuel 
concentration  (Fig.  3)  continues  to  fall  even  for  rich 
mixtures.  The  same  type  of  behavior  should  also 
occur  for  C2H2-air,  also  incorrectly  predicted  by  the 
simplified  mechanism. 

Author’s  Reply.  It  is  certainly  true  that  the  jus¬ 
tification  for  using  a  one-step  chemical  kinetic  model 
is  very  weak.  However,  there  is  also  little  justifi¬ 
cation  for  using  a  more  detailed  chemical  model 
without  taking  into  account  the  complex  3-D  struc¬ 
ture  of  the  detonation  wave.  This  is  not  to  say  that 
such  modefling  is  not  valuable,  but  such  parame¬ 
ters  as  the  cell  size,  critical  tube  diameter  and  crit¬ 
ical  energy  can  only  be  calculated  when  a  com¬ 
bined  gas  dynamic-chemical  kinetics  model,  which 
takes  into  account  the  cellular  detonation  structure, 


has  been  develped  Until  such  time,  the  correlation 
of  these  parameters  using  chemical  kinetic  models 
is  purely  empirical,  so  that  the  justification  for  us¬ 
ing  one  mode)  over  another  can  only  be  deter¬ 
mined  from  experiment.  In  the  particular  case  of 
CjH(-air  mixtures,  the  one-step  kinetic  model  ap¬ 
pears  to  correlate  the  results  for  lean  mixtures. 
However,  as  Dr.  Westbrook  states,  it  fails  to  ac¬ 
count  for  the  increase  in  cell  sizes  for  rich  mixtures 
recently  observed  by  Knystautas  et  al.  at  McGill 
University.  For  C2II2-air,  the  available  data  is  not 
extensive  enough  to  determine  which  correlation  is 
better. 

Although,  the  success  of  these  simplified  models, 
including  Dr.  Westbrook's  model,  may  be  fortu¬ 
itous  they  serve  as  useful  guides  for  experimental 
investigations. 

• 

A.  Lifshitz,  Hebrew  University,  Israel.  The  in¬ 
duction  time  for  the  COS  +  02  — »  CO  +  S02  re¬ 
action  system  has  a  very  weak  temperature  depen¬ 
dence;  —  17  Kcal/mole.  Since  the  kinetics  of  these 
systems  is  very  well  known  it  is  highly  recom¬ 
mended  for  all  cell  size  calculation  by  kinetic  and 
hydrodynamic  interaction.  The  system  will  only  be 
weakly  dependent  on  the  exact  temperature  profile 
across  the  cell. 

Authors  Reply,  litis  is  an  excellent  proposal,  since 
presumably  for  this  sytem  the  hydrodynamic  part 
of  the  interaction  can  be  identified  and  validated 
before  attempting  to  include  more  temperature 
sensitive  chemical  systems. 
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mgs,  S.  associated  with  the  detonations  were  measured'  from  smoked  foils  mounted  in  the 
tubes.  These  measurements  demonstrate  that  the  empirical  relation  d,  ■  13  S  provides  a 
good  correlation  between  the  critical  tube  diameter  and  the  cell  size  over  a  wide  range  of 
ethylene-air  compositions.  However,  a  better  understanding  of  the  coupled  chemical-gas- 
dynamic  processes  within  the  cell  is  required  in  order  to  clarify  the  link  between  dell  size 
and  chemical  kinetics. 

In  addition  to  the  open  tube  tests,  transmission  of  detonations  through  circular  orifice 
holes  was  investigated  in  both  ethylene-air  and  acetylene-air  mixtures.  Based  on  these  tests 
it  as  concluded  that  the  critical  tube  and  critical  orifice  diameters  are  equal. 
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DYNAMIC  PARAMETERS  OF 
GASEOUS  DETONATIONS 


John  H.  S.  Lee 

Department  of  Mechanical  Engineering,  McGill  University,  Montreal, 
Quebec,  Canada  H3A  2K6 


INTRODUCTION 

In  addition  to  gases,  flammable  liquids  and  solids  in  the  form  of  fine 
droplets  and  dust  particles  also  form  explosive  mixtures  with  air.  An 
explosive  mixture  can,  in  general,  support  two  modes  of  combustion.  The 
slow  laminar  deflagration  mode  is  at  one  extreme;  here  the  flame 
propagates  at  typical  velocities  of  the  order  1ms-1  relative  to  the  unburned 
gases  and  there  is  negligible  overpressure  development  when  the  explosion 
is  unconfined.  At  the  other  extreme  is  the  detonation  mode,  in  which  the 
detonation  wave  propagates  at  about  2000  ms-1  accompanied  by  an 
overpressure  rise  of  about  20  bars  across  the  wave.  The  propagation  of 
laminar  deflagrations  is  governed  by  the  molecular  diffusion  of  heat  and 
mass  from  the  reaction  zone  to  the  unburned  mixture.  The  propagation  of 
detonations  depends  on  the  adiabatic  shock  compression  of  the  unburned 
mixtures  to  elevated  temperatures  to  bring  about  autoignition.  The  very 
strong  exponential  temperature  dependence  of  chemical  reaction  rates  in 
general  makes  possible  the  rapid  combustion  in  the  detonation  mode.  Two- 
phase  liquid  droplets  or  dust-air  mixtures  are  similar,  but  they  require  more 
physical  processes  (e.g.  droplet  break-up,  phase  change,  mixing,  etc.)  prior 
to  combustion.  Thus,  characteristic  time  or  length  scales  associated  with 
the  combustion  front  are  usually  much  larger  than  those  of  homogeneous 
gaseous  fuel-air  mixtures.  The  essential  mechanisms  of  propagation  of  the 
combustion  waves,  however,  are  similar.  In  between  the  two  extremes  of 
laminar  deflagration  and  detonation,  there  is  an  almost  continuous 
spectrum  of  burning  rates  where  turbulence  plays  the  dominant  role  in  the 
combustion  process.  Due  to  space  limitations,  only  homogeneous  gaseous 
fuel-air  detonations  are  considered  in  this  article. 
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The  problem  of  fuel-air  explosions  has  received  considerable  attention  in 
recent  years  in  connection  with  the  safety  aspects  of  large-scale  transport 
and  storage  of  liquified  petroleum  or  natural  gases.  The  Three  Mile  Island 
incident  also  raises  the  question  of  the  possibility  and  consequences  of  a 
hydrogen-air  explosion  inside  a  nuclear-reactor  containment.  Two-phase 
liquid  sprays  and  dust  explosions  are  also  of  concern  in  accidental  fuel 
release  in  petrochemical  plants  and  in  grain  elevators  and  coal  mines.  As  a 
result  of  the  current  worldwide  interest  in  fuel-air  explosions,  active 
research  involving  large-scale  field  tests  has  been  carried  out  in  recent  years. 
In  this  paper,  we  single  out  the  gaseous  fuel-air  detonation  part  of  fuel-air 
explosions,  reporting  on  the  progress  made  and  the  current  state  of  the  art. 

Accidental  fuel-air  explosions  are  discussed  in  the  recent  monographs  of 
Gugan  (19781.  Bartknecht  (1978).  and  Baker  et  al.  (1983),  as  well  as  in  an 
excellent  general  review  article  by  Strehlow  (1980).  In  addition,  the 
proceedings  of  a  recent  specialist  meeting  on  fuel-air  explosions  contains 
progress  reports  of  current  activities  at  various  laboratories  around  the 
world  (Lee  &  Guirao  1982).  The  present  paper,  therefore,  avoids  dupli¬ 
cation  of  material  that  can  be  found  elsewhere.  Efforts  are  made,  however, 
to  keep  the  paper  self-contained. 

THE  DYNAMIC  PARAMETERS  OF  DETONATIONS 

The  classical  Chapman-Jouguet  theory,  in  essence,  seeks  the  unique 
solution  of  the  one-dimensional  conservation  equations  across  the  det¬ 
onation  front  in  which  the  flow  behind  the  wave  is  sonic.  It  does  not  require 
a  knowledge  of  the  structure  of  the  wave  itself,  and  involves  only  an 
equilibrium  thermodynamic  calculation  for  the  detonation  states  (i.e.  the 
detonation  velocity,  pressure,  temperature,  am.  density  ratios  across  the 
wave,  and  the  equilibrium  composition  of  the  product  gases).  The 
equilibrium,  or  “static."  detonation  states  obtained  from  the  classical 
C  hapman-Jouguet  theory  agree  surprisingly  well  with  experimental  obser¬ 
vations,  even  in  near-limit  conditions  when  the  flow  structure  near  the  front 
is  highly  three-dimensional. 

The  parameters  requiring  a  knowledge  of  the  structure,  and  hence  the 
chemical  reaction  rates,  are  the  detonability  limits,  the  initiation  energy,  the 
critical  tube  diameter,  and  the  thickness  of  the  reaction  zone.  We  henceforth 
refer  to  these  parameters  as  the  dynamic  detonation  parameters  to 
distinguish  them  from  the  equilibrium  “static”  detonation  states  obtained 
from  the  Chapman-Jouguet  theory.  A  one-dimensional  model  for  the 
detonation  structure  was  proposed  in  the  early  1940s  by  Ya.  B.  Zel’dovich, 
W.  Doring,  and  J.  von  Neumann  (i.e.  the  ZDN  model)  It  assumes  a  planar 
shock  followed  by  a  reaction  zone  after  an  induction  delay  period.  This 
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model  permits  the  computation  of  the  dynamic  parameters  when  a  model 
for  the  physical  processes  involved  is  given.  However,  results  from  such 
computations  using  the  ZDN  model  for  the  detonation  structure  are  quite 
far  from  those  obtained  from  experiments,  mainly  because  the  ZDN 
structure  is  unstable  and  is  never  observed  experimentally  except  under 
transient  conditions.  This  is  in  accord  with  experimental  observations  that 
all  self-sustained  detonations  have  a  three-dimensional  cellular  structure.  A 
quantitative  theory  for  the  real  cellular  detonation  structure  has  yet  to  be 
developed.  Thus,  it  seems  rather  surprising  that  a  century  after  the 
detonation  phenomenon  was  first  identified,  and  over  eighty  years  since  the 
successful  Chapman-Jouguet  theory  for  the  prediction  of  the  equilibrium 
properties  was  proposed,  there  is  still  no  quantitative  theory  for  estimating 
these  dynamic  detonation  parameters.  They  remain  as  experimental 
parameters,  and  their  measurements  are  not  without  some  fundamental 
difficulties  either. 

Detonation  research  up  until  the  early  1970s  was  mostly  aimed  at 
furthering  our  understanding  of  the  physical  and  chemical  processes 
involved  in  the  various  aspects  of  the  detonation  phenomena  (i.e.  transition 
from  deflagration  to  detonation,  onset  of  detonation,  blast  initiation,  shock 
interaction  processes  leading  to  the  so  called  multiheaded  or  cellular 
structure,  etc.).  Thus,  sensitive  mixtures  of  fuels  (e.g.  H2,  C2H2)  with  pure 
oxygen  were  used  in  general,  and  the  experiments  were  usually  carried  out 
at  subatmospheric  initial  pressures  for  convenience.  Measurements  of  the 
dynamic  parameters  were  also  mostly  made  on  fuel-oxygen  mixtures  at 
subatmospheric  initial  pressures,  since  the  objective  of  these  measurements 
was  to  study  the  dependence  of  these  dynamic  detonation  parameters  on 
the  other  experimental  variables.  As  a  result,  there  are  very  few  exper¬ 
imental  data  on  these  dynamic  parameters,  particularly  on  practical 
systems,  to  permit  even  empirical  correlations  to  be  made. 

In  Volume  5  of  this  series,  Oppenheim  &  Soloukhin  (1973)  gave  a  very 
comprehensive  historical  summary  of  detonation  research  up  to  that  time, 
as  well  as  an  extensive  bibliography  of  the  early  works.  Of  particular 
interest  in  their  article  are  the  beautiful  laser-schlieren  cinematographic 
records  of  the  various  physical  processes  associated  with  the  phenomena  of 
transition  and  initiation,  as  well  as  the  detailed  cellular  structure  of  the 
detonation  front  itself.  These  exceptional  experimental  records  provide  a 
good  qualitative  understanding  of  the  detonation  phenomena,  and  readers 
are  encouraged  to  refer  to  this  article  for  a  physical  background  on  the 
subject.  Other  noteworthy  review  articles  summarizing  the  detailed 
structure  of  cellular  detonations  are  those  of  Edwards  (1969)  and  Strehlow 
(1968).  The  recent  monograph  by  Fickett  &  Davis  (1979)  also  gives  a  good 
comprehensive  review  of  the  stability  and  structure  of  cellular  detonations. 


More  recent  reviews  emphasizing  the  direct  initiation  and  the  transition 
from  deflagration  to  detonation  have  been  given  by  Lee  (1977)  and  Lee  & 
Moen  (1980). 

It  may  be  concluded  that  by  the  early  1970s  all  aspects  of  the  detonation 
phenomena  were  quite  well  understood,  at  least  on  i  qualitative  basis. 
There  is,  however,  a  strong  deficiency  in  reliable  experimental  data  on  the 
dynamic  detonation  parameters,  particularly  in  fuel-air  mixtures  at 
atmospheric  initial  pressures  that  are  of  practical  interest.  There  are  also 
very  few  results  on  fuel-air  detonations  due  to  the  fact  that  very-high-energy 
sources  such  as  solid  explosive  charges  must  be  used  for  the  initiation  as  a 
result  of  their  relative  insensitivity.  Also,  the  size  of  the  apparatus  required 
necessitates  that  the  experiments  be  carried  out  in  the  Held.  However,  as  a 
result  of  the  current  interest  in  accidental  fuel-air  explosions,  a  number  of 
large-scale  experiments  have  been  carried  out  in  the  past  decade  on  fuel-air 
detonations.  We  emphasize  here  the  discussion  of  these  new  results  and 
their  correlations,  and  point  out  the  progress  made  toward  the  prediction  of 
the  dynamic  detonation  parameters. 

THE  DETONATION  CELL  SIZE  X 

Perhaps  the  most  important  dynamic  parameter  is  the  average  size  A  of  a 
cell  of  the  detonation  front.  Because  of  its  fundamental  significance,  it  is 
worthwhile  to  review  briefly  the  essential  features  of  a  cellular  detonation 
front.  Figure  1  shows  the  pattern  made  by  the  normal  reflection  of  a 
detonation  on  a  glass  plate  coated  lightly  with  carbon  soot,  which  may  be 
from  either  a  wooden  match  or  a  kerosene  lamp.  The  cellular  structure  of 
the  detonation  front  is  quite  evident.  If  a  similarly  soot-coated  polished 
metal  (or  mylar)  foil  is  inserted  into  a  detonation  tube,  the  passage  of  the 
detonation  wave  will  leave  a  characteristic  “fish-scale”  pattern  on  the 
smoked  foil.  Figure  2  is  a  sequence  of  Iaser-schlieren  records  of  a  detonation 
wave  propagating  in  a  rectangular  tube.  One  of  the  side  windows  has  been 
coated  with  smoke,  and  the  fish  scale  pattern  formed  by  the  propagating 
detonation  front  is  very  well  illustrated.  The  detailed  shock  configuration  of 
the  cellular  detonation  front  itself  is  illustrated  by  the  interferogram  shown 
in  Figure  3.  The  direction  of  propagation  of  the  detonation  is  toward  the 
right.  As  can  be  seen  in  the  sketch  at  the  top  left  corner,  there  are  two  triple 
points.  At  the  first  triple  point  A,  AI  and  AM  represent  the  incident  shock 
and  Mach  stem  of  the  leading  front,  while  AB  is  the  reflected  shock.  Point  B 
is  the  second  triple  point  of  another  three-shock  Mach  configuration  on  the 
reflected  shock  AB;  the  entire  shock  pattern  represents  what  is  generally 
referred  to  as  a  double  Mach  reflection.  The  hatched  lines  denote  the 
reaction  front,  while  the  dash-dot  lines  represent  the  shear  discontinuities 
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or  slip  lines  associated  with  the  triple-shock  Mach  configurations.  The 
entire  front  ABCDE  is  generally  referred  to  as  the  transverse  wave,  and  it 
propagates  normal  to  the  direction  of  the  detonation  motion  (down  in  the 
present  case)  at  about  the  sound  speed  Cx  of  the  hot  product  gases.  It  has 
been  shown  conclusively  that  it  is  the  triple-point  regions  at  A  and  B  that 
"write”  on  the  smoked  foil.  The  exact  mechanics  of  how  the  triple-point 
region  does  the  writing  is  not  clear.  It  has  been  postulated  that  the  high 
shear  at  the  slip  discontinuity  causes  the  soot  particles  to  be  erased.  Figure  4 
shows  a  schematic  of  the  motion  of  the  detonation  front.  The  fish-scale 
pattern  is  a  record  of  the  trajectories  of  the  triple  points.  It  is  important  to 
note  the  cyclic  motion  of  the  detonation  front.  Starting  at  the  apex  of  the  cell 
at  A.  the  detonation  shock  front  is  highly  overdriven,  propagating  at  about 


Fit/urc  I  End-on  pattern  from  the  normal  reflection  of  a  cellular  detonation  on  a  smoked 
glass  plate. 
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Figure  2  Laser-schlieren  cinematography  of  a  propagating  detonation  in  low-pressure 
H;-02  mixtures,  with  fish-scale  pattern  shown  on  a  soot-coatui  window  (courtesy  of 
A.  K.  Oppenheim). 

1.6  times  the  equilibrium  Chapman-Jouguet  detonation  velocity.  Toward 
the  end  of  the  cell  at  D,  the  shock  has  decayed  to  about  0.6  times  the 
Chapman-Jouguet  velocity  before  it  is  impulsively  accelerated  back  to  its 
highly  overdriven  state  when  the  transverse  waves  collide  to  start  the  next 
cycle  again.  For  the  first  half  of  the  propagation  from  A  to  BC,  the  wave 
serves  as  the  Mach  stem  to  the  incident  shocks  of  the  adjacent  cells.  During 
the  second  half  from  BC  to  D,  the  wave  then  becomes  the  incident  shock  to 
the  Mach  stems  of  the  neighboring  cells.  Details  of  .he  variation  of  the 
shock  strength  and  chemical  reactions  inside  a  cell  can  be  found  in  a  recent 
paper  by  Libouton  et  al.  ( 1981).  AD  is  usually  defined  a  >  the  length  Lc  of  the 
cell,  and  BC  denotes  the  cell  diameter  (also  referred  to  as  the  cell  width  or 
the  transverse-wave  spacing).  The  average  velocity  of  the  wave  is  close  to 
the  equilibrium  Chapman-Jouguet  velocity. 

We  thus  see  that  the  motion  of  a  real  detonation  front  is  far  from  the 
steady  and  one-dimensional  motion  given  by  the  ZDN  model.  Instead,  it 
proceeds  in  a  cvelic  manner  in  which  the  shock  vclocit  /  fluctuates  within  a 
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cell  about  the  equilibrium  Chapman-Jouguet  value.  Chemical  reactions  are 
essentially  complete  within  a  cycle  or  a  cell  length.  However,  the 
gasdynamic  flow  structure  is  highly  three-dimensional,  and  full  equili¬ 
bration  of  the  transverse  shocks,  so  that  the  flow  becomes  essentially  one¬ 
dimensional,  will  probably  take  an  additional  distance  of  the  order  of  a  few 
more  cell  lengths. 


figure  3  Interfcrogram  of  the  detailed  double  Mach-refiection  configuration  of  the  structure 
of  a  cellular  detonation  front  (courtesy  of  D.  H.  Hdwards). 


Figure  4  Schematic  of  the  propagation  of  a  cellular  detonation  front,  showing  the  trajectories 
of  thje  triple  points. 

From  both  the  cellular  end-on  or  the  axial  fish-scale  smoke  foil,  the 
average  cell  size  X  can  be  measured.  The  end-on  record  gives  the  cellular 
pattern  at  one  precise  instant.  The  axial  record,  however,  permits  the 
detonation  to  be  observed  as  it  travels  along  the  length  of  the  foil.  It  is  much 
easier  by  far  to  pick  out  the  characteristic  cell  size  X  from  the  axial  record ; 
thus,  the  end-on  pattern  is  not  used,  in  general,  for  cell-size  measurements. 

Early  measurements  of  the  cell  size  have  been  carried  out  mostly  in  low- 
pressure  fuel-oxygen  mixtures  diluted  with  inert  gases  such  as  He,  Ar,  and 
N2  (Strehlow  &  Engel  1969).  The  purpose  of  these  investigations  is  to 
explore  the  details  of  the  detonation  structure  and  to  find  out  the  factors 
that  control  it.  It  was  not  until  very  recently  that  Bull  et  al.  (1982)  made 
some  cell-size  measurements  in  stoichiometric  fuel-air  mixtures  at  atmos¬ 
pheric  pressure.  Due  to  the  fundamental  importance  of  the  cell  size  in  the 
correlation  with  other  dynamic  parameters,  a  systematic  program  has 
been  carried  out  by  Knystautas  to  measure  the  cell  size  of  atmospheric  fuel- 
air  detonations  in  all  the  common  fuels  (e.g,  H2,  C2H2,  C2H4,  C3H6,  C2H6, 
C3H8,  C4H10,  and  the  welding  fuel  MAPP)  over  the  entire  range  of  fuel 
composition  between  the  limits  (Knystautas  et  al.  1983).  Stoichiometric 
mixtures  of  these  fuels  with  pure  oxygen,  and  with  varying  degrees  of  N2 
dilution  at  atmospheric  pressures,  were  also  studied  (Knystautas  et  al. 
1982).  To  investigate  the  pressure  dependence,  Knystautas  et  al.  (1982)  have 
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also  measured  the  cell  size  in  a  variety  of  stoichiometric  fuel-oxygen 
mixtures  at  initial  pressures  10  <  p0  <  200  torr.  Figure  5  gives  the  typical 
U-shaped  curves  of  the  variation  of  the  cell  size  X  for  some  of  the  common 
hydrocarbon  fuel-air  mixtures  with  the  equivalence  ratio  4>  (i.e.  <p  =  1  for 
stoichiometric  composition)  of  the  mixture.  The  minimum  cell  size  usually 
occurs  at  about  the  most  detonable  composition  {(p  =  1).  The  cell  size  X  is 
representative  of  the  sensitivity  of  the  mixture.  Thus,  in  descending  order  of 
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Figure  5  Cell  size  of  fuel-air  mixtures  at  atmospheric  piessure. 


sensitivity,  we  have  C2H2,  H2,  C2H4,  and  the  alkanes  C3H8,  C2H6,  and 
C4H10.  Methane  (CH4),  although  belonging  to  the  same  alkane  family,  is 
exceptionally  insensitive  to  detonation,  with  an  estimated  cell  size  X  ~  33 
cm  for  stoichiometric  composition  as  compared  with  the  corresponding 
value  of  X  ^  5.35  cm  for  the  other  alkanes.  That  the  cell  size  X  is 
proportional  to  the  induction  time  of  the  mixture  had  been  suggested  by 
Shchelkin  &  Troshin  (1965)  long  ago.  However,  to  compute  an  induction 
time  requires  that  the  model  for  the  detonation  structure  be  known,  and  no 
theory  exists  as  yet  for  the  real  three-dimensional  structure.  Nevertheless, 
one  can  use  the  classical  ZDN  model  for  the  structuie  and  compute  an 
induction  time  or,  equivalently,  an  induction-zone  length  l.  While  this  is  not 
expected  to  correspond  to  the  cell  size  x  (or  cell  length  Lc),  it  may  elucidate 
the  dependence  of  X  on  /  itself  (e.g.  a  linear  dependence  /.  =  Al  as  suggested 
by  Shchelkin  &  Troshin  1965).  Westbrook  (1982,  Westbrook  &  Urtiew 
1 982)  has  made  computations  of  the  induction-zone  length  /  using  the  ZDN 
model,  but  his  calculations  are  based  on  a  constant-volume  process  after 
the  shock,  rather  than  integration  along  the  Rayleigh  line.  Very  detailed 
kinetics  of  the  oxidation  processes  are  employed.  By  matching  with  one 
experimental  point,  the  proportionality  constant  A  can  be  obtained.  This  is 
done  in  Figure  5  by  using  Westbrook’s  data  for  l  and  matching  at  </>  =  1. 
The  solid  lines  in  the  figure  represent  the  linear  relationship  X  =  Al,  where  / 
is  computed  from  detailed  kinetics.  The  general  dependence  of  X  on  0  is 
reproduced  qualitatively,  but  quantitatively  the  agreement  is  not  as  good  as 
it  appears  because  errors  are  hidden  by  the  very  steep  nature  of  the  U- 
shaped  curve.  The  agreement,  however,  is  whhin  one  order  of  magnitude. 
The  constant  A  differs  for  different  gas  mixtures  (e.g.  A  —  1 0.14  for  C2H4,  A 
=  52.23  for  H2);  thus,  the  three-dimensional  gasdynamic  processes  cannot 
be  represented  by  a  single  constant  alone  over  a  range  of  fuel  composition 
for  all  the  mixtures.  The  chemical  reactions  in  a  detonation  wave  are 
strongly  coupled  to  the  details  of  the  transient  gasdynamic  processes,  with 
the  end  product  of  the  coupling  being  manifested  by  a  characteristic 
chemical  length  scale  X  (or  equivalently  Lc)  or  time  scale  tc  =  X/Cx  (where 
C,  denotes  the  sound  speed  in  the  product  gases,  which  is  approximately 
the  velocity  of  the  transverse  waves)  that  describes  the  global  rate  of  the 
chemical  reactions.  Since  X  ^  0.6 Lc  and  Cx  ~  0.5D,  where  D  is  the 
Chapman-Jouguet  detonation  velocity,  we  have  rc  ~  LJD,  which  cor¬ 
responds  to  the  fact  that  the  chemical  reactions  are  essentially  completed 
within  one  cell  length  (or  one  cycle). 

It  appears  that  a  correct  model  for  the  cellular  structure  must  be  a  time- 
dependent  one  in  which  the  nonlinear  coupling  mechanism  between 
gasdynamics  and  chemical  kinetics  can  be  properly  modeled.  Computer 
codes  for  transient  reacting  flows  with  shock  waves  are  currently  available. 
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The  feasibility  of  modeling  the  real  cellular  structure  of  a  two-dimensional 
detonation  front  has  also  been  demonstrated  recently  by  Taki  &  Fujiwara 
(1981)  and  Oran  et  al.  (1981).  However,  an  extension  to  three  dimensions 
with  the  appropriate  spatial  resolutions  to  account  for  the  details  of  the 
complex  double  Mach  configurations,  although  possible  in  principle, 
would  require  computer  time  and  storage  capacity  much  in  excess  of  what 
even  the  current  fast  machines  are  capable  of.  Thus,  the  use  of  numerical 
simulation  as  a  tool  for  obtaining  cell-size  data  would  not  be  practical,  and 
direct  experimental  measurement  remains  as  the  most  convenient  means  of 
determining  he  cell  size  /.  Although  other  methods  have  been  tried  to 
measure  the  cell  size  directly,  the  simple  smoke-foil  technique  has  not  been 
improved  upon  to  date  and  still  remains  as  the  only  successful  method.  For 
near-limit  fuel-air  detonations,  where  the  cells  are  very  large  and  foils  on  the 
order  of  meters  in  length  have  to  be  used,  the  actual  deposition  of  a  uniform 
coating  of  soot  from  numerous  burners  poses  a  rather  difficult  engineering 
problem.  A  more  serious  problem  in  the  measurement  of  cell  size  is  the 
actual  interpretation  of  the  fish-scale  pattern.  For  most  fuel-air  mixtures, 
this  pattern  is  highly  irregular,  and  thus  the  selection  of  the  “correct”  cell 
size  requires  a  certain  amount  of  experience.  This  introduces  a  subjective 
element  into  the  measurement  of  A,  which  has  to  be  resolved.  The  use  of  long 
foils  so  that  the  wave  is  recorded  over  a  long  travel  length  certainly 
facilitates  the  interpretation  of  the  foil.  Alternatively,  if  a  large  number  of 
experiments  are  made  in  order  to  accumulate  an  ensemble  of  records  under 
identical  conditions,  then  this  also  improves  the  accuracy  of  the  measure¬ 
ment.  However,  other  techniques  must  be  developed  to  facilitate  the 
measurement  of  this  important  fundamental  dynamic  parameter,  which 
characterizes  the  real  structure  of  the  detonation  wave. 

THE  CRITICAL  TUBE  DIAMETER 

Another  important  dynamic  parameter  that  has  received  considerable 
attention  in  recent  years  is  the  so-called  critical  tube  diameter. 
Experimentally,  it  is  found  that  if  a  planar  detonation  wave  propagating  in 
a  circular  tube  emerges  suddenly  into  an  unconfined  volume  containing  the 
same  mixture,  the  planar  wave  will  transform  into  a  spherical  wave  if  the 
tube  diameter  d  exceeds  a  certain  critical  value  dc  (i.e.  d  >  dc).  If  d  <  dc,  the 
expansion  waves  will  decouple  the  reaction  zone  from  the  shock,  and  a 
spherical  deflagration  wave  results.  Excellent  schlieren  records  of  these 
transmission  phenomena  for  d  >  dc  and  d  <  dc  can  be  found  in  the  review 
by  Oppenheim  &  Soloukhin  (1973). 

Perhaps  the  most  important  progress  in  recent  years  has  been  in  the 
linking  together  of  the  various  dynamic  parameters.  This  work  originates 


from  the  demonstration  of  a  universal  correlation  between  the  critical  tube 
diameter  and  the  cell  size,  an  observation  first  made  some  20  years  ago  by 
Mitrofanov  &  Soloukhin  (1965).  In  studying  the  diffraction  of  a  planar 
detonation  wave  as  it  emerges  from  a  circular  tube  into  unconfined  space 
containing  the  same  mixture,  they  noticed  that  the  critical  tube  diameter  is 
about  13  times  the  cell  size  of  the  mixture  (i.e.  dc  =  132).  For  a  square  tube, 
the  critical  width  We  is  found  to  be  of  the  order  of  10  times  the  cell  size  (i.e. 
Wc  =  10/.).  Since  they  based  their  observation  on  one  particular  mixture  of 
stoichiometric  acetylene  and  oxygen  over  a  narrow  range  of  initial 
pressures  (i.e.  p0  <  100  torr),  the  significance  of  this  observation  was  not 
recognized.  Over  a  decade  elapsed  before  Edwards  et  al.  (1979)  brought  up 
this  observation  again.  In  repeating  Mitrofanov  &  Soloukhin’s  experi¬ 
ments,  they  confirmed  the  earlier  result  that  dc  =  132.  Furthermore,  they 
also  pointed  out  that  there  is  no  reason  why  this  correlation  should  not  be 
applicable  to  all  other  detonative  systems.  Moen  et  al.  (198 1)  followed  up  on 
this  suggestion  and  demonstrated  the  validity  of  the  132.  correlation  for 
stoichiometric  ethylene-oxygen  mixtures  with  nitrogen  dilution.  A  sys¬ 
tematic  experimental  program  was  then  carried  out  by  Knystautas  et  al.,  in 
which  cell  sizes  and  critical  tube  diameters  were  simultaneously  measured 
for  all  the  common  hydrocarbons  (i.e.  H2,  C2H2,  C2H4,  C3H6,  C2H6,  C3H8, 
CH4,  and  MAPP).  Cell  sizes  were  measured  in  stoicniometric  mixtures  of 
these  fuels  with  pure  oxygen  at  low  initial  pressures  (p0  <  150  torr),  as  well 
as  in  stoichiometric  fuel-oxygen  mixtures  at  atmospheric  pressure,  but  with 
different  oegrees  of  nitrogen  dilution.  A  comparison  with  the  direct 
measurements  of  the  critical  tube  diameter  confirmed  the  universality  of  the 
13/.  correlation.  A  comprehensive  report  of  the  detailed  results  of  this 
important  study  car  be  found  in  Knystautas  et  al.  (1981).  The  cell-size 
measurement  has  recently  been  further  extended  to  fuel-air  mixtures  over 
the  range  of  fuel  concentration  between  the  limits.  We  have  seen  some  of 
these  cell-size  data  already  in  Figure  5.  There  have  also  been  a  number  of 
large-scab  field  experiments  on  the  direct  measurements  of  the  critical  tube 
diameter  in  fuel-air  mixtures  (Rinnan  1982,  Guirao  et  al.  1982,  Moen  et  al. 
1 983).  The  results  are  summarized  in  Figure  6,  which  shows  the  dependence 
of  the  critical  tube  diameter  dc  on  the  equivalence  ratio  (j)  (i.e.  mixture 
composition)  for  some  of  these  common  fuel-air  mixtures.  The  solid  curves 
are  the  cell-size  data  from  Figure  5  multiplied  by  1 3.  As  can  be  observed,  the 
agreement  is  sufficiently  good  to  confirm  the  universal  nature  of  the  dc 
-  13/.  correlation.  Some  doubt  has  been  expressed  by  Moen  et  al.  (1983), 
who  simultaneously  measured  both  cell  size  and  critical  tube  diameter. 
They  reported  that,  based  on  their  own  cell-size  data,  142  <  dc  ;$  242. 
However,  their  critical  tube  diameter  data  agree  perfectly  with  the  132 
correlation  when  the  cell-size  data  obtained  independently  by  Knystautas 
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et  al.  (1983)  are  used.  Knystautas  et  al’s  cell-size  data  are  judged  more 
reliable  in  that  very  long  smoke  foils  were  used  and  usually  a  number  of 
records  for  the  same  mixture  were  taken  to  ensure  a  correct  interpretation 
of  the  true  cell  size.  Furthermore,  independent  estimates  of  the  cell  size  have 
also  been  obtained  from  the  frequencies  of  the  pressure  oscillations  from  the 
pressure  profile  recorded  by  the  transducers  mounted  on  the  detonation 
tube.  An  assessment  of  the  data  available  thus  far  indicates  that  dc  =  \3X  is 
indeed  a  general  relationship,  at  least  for  the  range  of  fuels  tested.  It  should 
be  pointed  out  that  the  dc  =  13X  correlation  is  extremely  useful  in  practice. 
The  cell  size  /.  can  be  determined  relatively  easily  from  laboratory-scale 
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Figure  6  Comparison  of  experimental  data  for  the  critical  tube  diameter  in  fuel-air  mixtures 
with  the  13/.  correlation. 


detonation  tubes  under  easily  controlled  conditions,  thus  enabling  the 
critical  tube  diameter  to  be  obtained  immediately.  Direct  measurements  of 
the  critical  tube  diameter  require  rather  large-scale  field  experiments,  in 
general,  for  the  relatively  insensitive  fuel-air  mixtures. 

There  are  no  quantitative  theories  for  the  prediction  of  the  critical  tube 
diameter  as  yet.  From  physical  considerations,  a  qualitative  criterion  for 
the  transmission  can  be  developed.  When  the  detonation  emerges  from  the 
tube,  rarefaction  waves  are  generated  at  the  circumference  that  propagate 
toward  the  tube  axis  as  the  high-pressure  detonation  products  expand 
radially  outwards.  These  waves  cool  the  shocked  gases  and  thus  increase 
the  induction  time,  causing  the  reaction  zone  to  decouple  from  the  leading 
shock  front.  The  characteristic  time  for  this  gasdynamic  quenching  process 
will  be  of  the  order  of  the  tube  radius  divided  by  the  sound  speed  of  the 
product  gases  (i.e.  tc  =  RJCJ.  It  seems  reasonable  to  assume  that  if  the 
effective  thickness  of  the  detonation  wave  (i.e.  from  the  leading  shock  to  the 
equilibrium  Chapman-Jouguet  plane)  is  AH,  then  the  detonation  must 
propagate  a  distance  of  at  least  2AH  prior  to  the  rarefaction  waves  reaching 
the  tube  axis  for  it  not  to  be  quenched.  In  this  way,  at  least  a  small 
detonation  core  near  the  tube  axis  is  not  influenced  by  the  gasdynamic 
quenching,  and  this  core  serves  as  a  kernel  that  subsequently  develops  into 
a  spherical  wave.  The  time  fD  for  the  detonation  to  propagate  an  axial 
distance  of  2AH  will  simply  be  2Ah jD,  where  D  is  the  Chapman-Jouguet 
velocity.  Equating  the  two  characteristic  times  (i.e.  tD  =  tc)  gives  Rc 
=  2C,  Ah/D.  Since,  for  most  detonable  mixtures,  the  detonation  velocity  D 
^  2C,,  we  see  that  the  above  criterion  gives  Rc  ~  AH  or  the  critical  tube 
diameter  dc  ~  2AH.  The  effective  thickness,  or  the  so-called  hydrodynamic 
thickness,  of  a  detonation  wave  is  reported  by  Edwards  et  al.  (1976)  to  be 
between  2.5  to  4  cell  lengths,  or  5  to  8  cell  diameters.  They  base  this  result  on 
measurements  of  the  decay  of  the  transverse  pressure  vibrations  (i.e. 
equilibration  of  the  transverse  shocks)  in  oxyhydrogen  and  oxyacetylene 
mixtures.  Their  results  are  in  accord  with  earlier  independent  measure¬ 
ments  by  Vasiliev  et  al.  (1972).  Taking  an  averaged  value  of  AH  from 
Edwards  et  al’s  observations  (i.e.  AH  =  6.5A),  we  see  that  the  proposed 
criterion  (tJtD  =  1)  leads  to  the  result  dc  ~  2AH  =  13A,  which  is  in  perfect 
agreement  with  experiment.  From  the  above  discussions,  the  critical-tube- 
diameter  problem  is  thus  reduced  to  the  resolution  of  the  hydrodynamic  or 
effective  thickness  of  the  cellular  detonation  front.  Knowledge  of  AH  permits 
dQ  to  be  determined. 

The  recent  experimental  work  of  Liu  et  al.  (1983)  has  also  produced  some 
interesting  results  that  shed  further  light  on  the  critical-tube-diameter 
problem.  In  studying  the  transmission  through  an  orifice  plate,  instead  of 
from  a  straight  tube,  it  is  found  that  the  critical  orifice  diameter  is  identical 
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to  the  critical  tube  diameter.  Thus,  the  upstream  pressure  enhancement 
from  the  reflected  detonation  off  the  orifice  plate  does  not  influence  the 
transition  process  from  a  planar  to  a  spherical  wave.  The  reflected 
detonation  generates  a  radially  converging  shock  wave,  which,  like  the 
rarefaction  waves,  propagates  toward  the  tube  axis  at  about  the  sound 
speed  of  the  product  gases  as  a  result  of  the  high  temperatures  (and,  hence, 
weak  shock).  The  fact  that  there  is  no  difference  between  the  corresponding 
critical  orifice  and  the  critical  tube  diameters  suggests  that  the  unat¬ 
tenuated  core  near  the  tube  axis,  rather  than  the  details  of  the  exact 
gasdynamic  flow  structure  behind  the  diffracted  wave,  dominates  the 
development  of  the  spherical  detonation. 

Experiments  on  the  transmission  through  different-shaped  orifice  plates 
(i.e.  square,  rectangular,  triangular,  and  elliptical)  indicate  that  the  13A 
correlation  still  holds  if  an  effective  diameter  de{{,  which  is  the  mean  value  of 
the  largest  and  smallest  dimensions  that  characterize  the  orifice  opening,  is 
used  (i.e.  de{{  =  13A).  For  a  square  orifice  of  side  length  W,  the  smallest 
characteristic  dimension  would  be  W,  while  the  largest  would  be  the 
diagonal  y/lW.  Thus,  the  effective  diameter  of  the  square  orifice  would  be 
d'„=\/2(W+y/2W)^\2W.  Hence,  de{{  =  13A  gives  Wc  =  (13/1.2)A 
~  10.8/,  in  close  agreement  with  the  earlier  observations  of  Wc  =  1(U  by 
Mitrofanov  &  Soloukhin  (1965)  and  Edwards  et  al.  (1979).  The  effective 
diameter  is  also  the  mean  value  of  the  diameter  of  the  circumscribed  and 
inscribed  circles  of  the  orifice  opening.  It  would  appear  that  the  smallest 
dimension  should  be  more  important,  since  it  controls  the  time  for  the 
rarefaction  to  penetrate  the  wave.  However,  when  the  two  characteristic 
dimensions  are  not  too  different,  the  expansion  is  three-dimensional,  and 
thus  the  mean  value  would  be  more  appropriate  to  characterize  the 
gasdynamic  process. 

In  generalizing  the  results,  one  could  also  introduce  the  concept  of  wave 
curvature.  The  critical  condition  for  transmission  could  then  be  stated  in 
terms  of  a  minimum  radius  of  curvature  of  the  diffracted  wave.  If  the 
rarefaction  waves  give  rise  to  a  curvature  of  the  diffracted  detonation 
exceeding  a  certain  critical  value,  then  failure  of  the  wave  results.  The 
maximum  curvature  (or  minimum  radius  of  curvature  R)  would  be  of  the 
order  of  about  the  hydrodynamic  thickness  AH  (i.e.  R  =  AH  =  6.5A)  for  the 
three-dimensional  case  of  the  transformation  from  arbitrarily  shaped 
planar  to  spherical  geometry. 

For  rectangular  orifices  (or  slots)  where  the  aspect  ratio  L/W  becomes 
large  (e.g.  L/W  >  7)  and  the  planar  wave  transforms  to  a  cylindrical  wave, 
Liu  et  al.  (1983)  found  that  the  smaller  length  W  now  controls  the 
transmission  process,  rather  than  the  mean  value  of  the  two  characteristic 
dimensions  L  and  W.  Based  on  the  wave-curvature  concept,  the  cor- 


responding  radius  of  curvature  of  a  two-dimensional  cylindrical  wave 
would  be  one  half  that  of  the  spherical  wave  for  the  same  curvature.  Thus,  if 
R  —  Ah  for  the  three-dimensional  case,  R  =  AH/2  =  3.252  for  the  two 
dimensional  case.  This  is  in  agreement  with  the  experimental  result  of  Liu 
that  WQ  ~  3.252. 

Since  Liu’s  experiments  are  based  on  rectangular  slot  orifices  of 
L/W  <  7,  some  large-scale  experiments  on  two-dimensional  rectangular 
channels  with  L/W  sj  35  have  been  carried  out  recently  by  Benedick  et  al. 
( 1 983)  to  verify  the  Wc  ~  3/  limit.  Their  experiments  used  H2-air  and  C2H4- 
air  mixtures  at  atmospheric  pressures,  instead  of  the  H2-02-N2  and  C2H4- 
02-N2  mixtures  used  by  Liu.  The  results,  however,  confirmed  the  two- 
dimensional  limit  of  Wc  =  32  found  by  Liu.  Furthermore,  a  half-channel 
experiment  was  also  performed  that  demonstrated  symmetry  about  the 
centerline,  in  that  a  value  of  14^/2  =  1.52  was  obtained. 

These  very  recent  experimental  findings  are  extremely  interesting.  It  is 
clear  that  the  critical-tube-diameter  problem  provides  much  insight  into 
the  coupling  mechanisms  between  gasdynamics  and  chemical  kinetics.  The 
sudden  emergence  from  confined  to  unconfined  space,  in  essence,  subjects 
the  detonation  to  a  finite  perturbation.  Whether  it  survives  and  readjusts  to 
the  ne  geometry  or  fails  will  define  the  minimum  requirements  for  the  self- 
sustained  propagation  of  the  wave  under  its  new  environment.  The  recent 
study  of  Murray  &  Lee  (1983)  on  detonation  failure  when  the  detonation 
transmits  from  a  rigid  steel  tube  to  a  thin  plastic  tube  is  of  particular  interest 
in  this  respect.  By  using  plastic  tubes  of  different  wall  thickness,  the  severity 
of  the  gasdynamic  expansion  can  be  controlled,  and  the  failure  criterion  in 
this  case  can  then  be  related  to  the  wave  curvature.  It  should  be  noted  that 
even  in  rigid  tubes,  the  negative  displacement  thickness  of  the  boundary 
layer  causes  flow  divergence  much  like  the  yielding  plastic  wall  experiment 
of  M  urray  &  Lee  ( 1 983).  Thus,  rarefaction  waves  are  generated  at  the  walls, 
which  then  propagate  toward  the  tube  axis  and  cause  the  wave  to  be  curved. 
Direct  measurements  of  the  radius  of  curvature  of  a  detonation  wave  have 
been  carried  out  recently  by  Desbordes  et  al.  ( 1 98  i ).  They  reported  a  radius 
of  curvature  of  about  1  m  fora  planar  detonation  in  a  10-em-diameter  tube. 
The  corresponding  curvature  of  the  wave  is  much  smaller  than  the  critical 
value  for  the  failure  of  a  spherical  wave  (i.e.  l/R  ~  1/6.52).  Thus,  it  appears 
that  the  minimum  radius  of  curvature  of  a  spherical  wave  would  be  about 
the  hydrodynamic  thickness  A„  ~  3.5 Lc  ar  6.52.  For  waves  with  radii  of 
curvature  much  larger  than  the  critical  value,  a  velocity  deficit  results.  Thus, 
failure  in  rigid  and  plastic  tubes,  as  well  as  sudden  expansion  into 
unconfined  space,  are  all  similar  processes  and  differ  only  by  the  degree  of 
the  transverse  expansion  itself.  Since  the  lateral  expansion  results  in  a 
curved  front,  the  propagation  criterion  under  different  boundary  con¬ 
ditions  can  be  unified  under  the  wave-curvature  concept. 
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CRITICAL  ENERGY  FOR  DIRECT  INITIATION 

Experimentally,  it  is  found  that  for  a  given  mixture  at  given  initial 
conditions,  a  definite  quantity  of  energy  must  be  used  to  initiate  a 
detonation  “instantaneously.”  By  “instantaneously”  is  meant  that  the 
initial  strong  blast  wave  generated  by  the  powerful  igniter  upon  the  rapid 
deposition  of  its  energy  decays  asymptotically  to  a  Chapman-Jouguet 
detonation.  If  the  igniter  energy  is  less  than  a  certain  critical  value,  the 
reaction  zone  progressively  decouples  from  the  blast  as  it  decays  and  a 
deflagration  results.  Since  the  transition  from  deflagration  to  detonation 
under  unconfined  conditions  is  extremely  difficult,  if  even  possible, 
spherical  detonations  are  almost  always  initiated  in  practice  via  the  direct 
or  blast-initiation  mode.  Since  the  subject  of  direct  initiation  has  been 
thoroughly  reviewed  previously  (Lee  1977),  it  suffices  to  consider  only  the 
recent  developments  and  results. 

The  interest  in  unconfined  fuel-air  detonations  in  recent  years  has  led  to  a 
number  of  experimental  studies  on  the  measurement  of  the  critical 
initiation  energy.  The  most  extensive  of  these  measurements  for  fuel-air 
mixtures  are  those  by  Bull  et  al.  (1978).  The  dependence  of  the  critical 
energy  (or,  equivalently,  the  weight  of  the  solid  explosive  charge  used)  on 
the  equivalence  ratio  cf)  is,  like  all  dynamic  parameters,  in  the  form  of  U- 
shaped  curves  with  the  minimum  around  stoichiometric  composition  (i.e.  $ 
=  1).  The  minimum  energy  at  4>  —  1  has  been  used  as  a  relative  measure  of 
the  sensitivity  of  the  various  fuels  to  detonation  (Matsui  &  Lee  1978).  Also, 
by  specifying  somewhat  arbitrarily  an  upper  limit  of  the  critical  energy 
above  which  the  mixture  is  rendered  nondetonable,  the  L -shaped  curve  can 
be  used  to  define  the  detonability  limits  for  unconlined  detonations. 
Although  there  are  a  number  of  experimental  parameters  that  may 
influence  the  initiation  energy  (e.g.  the  type  of  igniter,  its  geometry,  energy¬ 
time  characteristics,  etc.),  it  can  be  concluded  that  as  long  as  the  time  for  the 
energy  deposition  is  short  compared  with  the  characteristic  time  of  the  blast 
wave  [i.e.  R0/C0 ,  where  R0  =  (£0/p0)1/3  is  the  explosion  length  and  C0  is  the 
sound  speed],  the  blast  wave  can  be  considered  as  an  ideal  point  spherical 
blast  characterized  by  the  total  energy  £0  only.  Thus,  experimental 
measurements  of  the  critical  energy  in  fuel-air  mixtures  using  concentrated 
solid  explosive  charges  satisfy  this  condition,  and  we  need  not  consider 
other  details  of  the  source  itself. 

The  various  theories  for  the  critical  energy  have  already  been  reviewed 
(Lee  1 977).  They  all  show  a  cubic  dependence  of  the  critical  energy  E0  on  the 
induction  time  x  for  spherical  detonations  (i.e.  £0  ~  r3).  However,  the  lack 
of  an  overall  chemical  length  or  time  that  characterizes  the  real  three- 
dimensional  cellular  structure  of  the  detonation  wave  prevents  all  these 
theoretical  efforts  from  being  completely  quantitative.  The  use  of  an 


experimental  data  point  is  usually  necessary  to  evaluate  a  certain  constant 
of  proportionality  between  £0  and  t3.  As  an  example,  consider  the  simplest 
theory,  based  on  Zel’dovich’s  criterion  (see  Lee  1977)  that  when  the 
initiating  blast  has  decayed  to  the  Chapman-Jouguet  strength,  the  decay 
time  must  be  of  the  order  of  the  induction  time.  Using  strong-blast  theory, 
Lee  (1977)  derived  the  following  expression  for  the  critical  energy  £0: 

E0  =  kjPoID^3l(j  +  3)/2V+l^+l,  (1) 

where  kj  =  1, 2n,  4n  for  j  —  0, 1, 2,  corresponding  to  the  planar,  cylindrical, 
and  spherical  geometry,  respectively;  1  is  a  numerical  constant  that  is  a 
function  of  the  specific  heat  ratio  y  (for  y  =  1.4, 1  =  0.423, 0.626, 1.212  for  j 
=  2,  1, 0,  respectively);  p0  is  the  initial  density;  D  is  the  Chapman-Jouguet 
velocity,  and  r  is  a  characteristic  chemical  time.  If  the  induction  time 
evaluated  at  the  shock  temperature  corresponding  to  a  Chapman-Jouguet 
detonation  is  used  for  t,  Equation  ( 1 )  gives  values  for  £0  about  three  orders 
of  magnitude  smaller  than  the  experimental  values.  A  more  appropriate 
overall  chemical  time  can  be  derived  from  the  cell  length,  which  is 
approximately  twice  the  cell  diameter  A,  because  chemical  reactions  are 
essentially  complete  within  one  cell  length.  Thus,  we  may  define 

tc  =  2  A/D, 

and  Equation  1  then  becomes 

£o  —  kjp0I(j  +  3Y+  lD2Aj+  \  (2a) 

£0  =  500npoID2P  for  j  —  2.  (2b) 

A  knowledge  of  the  cell  diameter  A  from  experiment  permits  the  evaluation 
of  the  critical  energy  £0.  The  Chapman-Jouguet  velocity  D  and  the  value  for 
the  constant  /  can  be  determined  easily  when  the  initial  state  and  the 
mixture  composition  are  specified. 

Another  simple  relationship  linking  the  cell  diameter  to  the  critical 
energy  can  be  obtained  via  the  critical  tube  diameter.  Based  on  the  work- 
done  concept,  Lee  &  Matsui  (1977)  derived  a  simple  expression  for  £0  and 
dc  as 


(3) 


where  pu  u,,  and  Cj  are  the  pressure,  particle  velocity,  and  sound  speed  at 
the  Chapman-Jouguet  plane.  Using  the  Rankine-Hugoniot  conditions 
across  a  Chapman-Jouguet  detonation,  Equation  (3)  may  also  be  written  as 


np0D2  _  21977rp0Z)2A3 
24-/(y+l)  c  _  24y(y  +  1 ) 


(4) 


when  the  13A  correlation  is  used  for  dc. 
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Recently,  Lee  et  al.  (1982)  have  proposed  a  more  direct  link  between  the 
initiation  energy  and  the  cell  size  via  the  critical  tube  diameter.  The  idea  is 
based  on  the  requirement  of  a  minimum  surface  energy  before  a  planar 
wave  can  evolve  into  a  spherical  wave  without  failure.  This  minimum 
surface  energy  should  correspond  to  the  area  of  the  critical  tube,  ndl/4. 
Thus,  in  blast  initiation,  the  blast  energy  E0  must  be  such  that  when  the 
wave  has  decayed  to  the  Chapman-Jouguet  strength,  the  surface  energy  of 
the  blast  sphere  must  be  at  least  proportional  to  71^/4.  Equating  directly 
the  surface  area  of  the  blast  sphere  to  the  area  of  the  critical  tube,  we  write 

4  tcR2  =  nd2/4, 
and  hence, 

Rc  —  dJ4, 

where  Rc  is  the  blast  radius  when  its  strength  has  decayed  to  the  Chapman- 
Jouguet  value  (i.e.  Rs  -*■  Rc,  Ms  -+  AfCJ,  where  Ms  is  the  shock  Mach 
number).  Using  strong-blast  theory,  one  can  easily  determine  the  blast 
energy  as  a  function  of  the  blast  radius,  and  the  resultant  expression  is  given 
by 

E0  =  4nyp0M2aI(m/4)\ 

=  (2\91/\6)np0ID2P,  (5) 

which  is  quite  similar  to  Equations  (2b)  and  (4).  Using  the  cell-size  data  of 
Knystautas  et  al.  (1982,  1983),  the  critical  energies  have  been  computed  for 
the  common  fuel-air  mixtures  where  experimental  data  are  available  (J.  E. 
Elsworth,  private  communication).  The  results  indicate  that  Equation  (5), 
based  on  the  surface-energy  concept,  gives  the  best  correspondence.  Figure 
7  shows  a  comparison  of  the  predicted  results  using  Equation  (5)  with  those 
obtained  by  Elsworth  from  direct  experimental  measurements  of  the 
critical  weight.  The  critical  energy  in  Figure  7  has  been  expressed  as  an 
equivalent  charge  weight  of  tetryl  ( 1  g  tetryl  is  equivalent  to  4270  J)  for  easy 
comparison  with  Elsworth’s  experimental  data.  The  agreement,  in  general, 
is  reasonably  good  in  view  of  the  simplicity  of  the  model.  No  experimental 
data  of  critical  energies  have  been  reported  for  C2H2-air  mixtures  to  permit 
a  comparison  with  the  present  prediction  from  Equation  (5).  For  most 
practical  situations.  Equation  (5)  can  be  used  to  predict  critical  initiation 
energies  (or  charge  weights)  with  quite  acceptable  accuracy. 


DETONABILITY  LIMITS 

Composition  limits  refer  to  the  minimum  and  maximum  fuel  concentration 
in  which  a  self-sustained  detonation  can  propagate.  From  the  previous 
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Figure  7  Comparison  of  critical-energy  data  for  fuel-air  mixtures  with  the  surface-energy 
theory. 

discussion,  it  is  clearly  meaningless  to  specify  the  composition  limit  without 
giving  the  details  of  the  boundary  conditions  or  the  environment  in  which 
the  wave  propagates.  Thus,  for  a  planar  detonation  propagating  in  a  rigid 
tube  of  circular  cross  section,  the  tube  diameter  must  be  specified 
simultaneously  with  the  composition  limits.  Strictly  speaking,  the  nature  of 
the  wall  also  influences  the  detonation  propagation.  Thus,  the  wall  material 
and  its  surface  roughness  must  also  be  specified.  In  general,  different  tubes 
of  different  diameters,  or  of  different  geometrical  cross-sectional  area,  wall 
roughness,  etc.,  will  all  have  their  respective  lean  and  rich  limits.  Cylindrical 
and  spherical  detonations  accordingly  will  also  have  their  particular 
limiting  concentrations.  It  appears  that  to  be  rid  of  the  boundary  problem 
completely,  only  cylindrical  or  spherical  detonation  limits  can  be  con¬ 
sidered  as  “true  limits”  characteristic  of  the  mixture.  Wt  shall  return  to  our 
discussion  of  the  limits  for  unconfined  detonations  later. 
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Consider  the  case  of  “smooth”  circular  tubes  where  only  the  tube 
diameter  need  be  specified  in  relation  to  the  limiting  composition.  There 
exists  the  problem  of  an  experimentally  derived  operational  definition  of 
the  limits.  Since  very  strong  initiators  must  be  used  for  the  near-limit 
mixtures,  the  first  question  to  be  resolved  is  how  long  must  the  tube  be 
before  the  influence  of  the  initiator  can  be  considered  negligible  and  the 
wave  is  truly  self-sustained.  Experiments  indicate  that  the  effect  of  the 
initiator  takes  an  extremely  long  time  to  decay  (e.g.  hundreds  of  tube 
diameters).  There  are  strong  indications  that  in  confined  tubes  the  single¬ 
head  spinning  mode  can  be  initiated  in  mixtures  outside  the  limits  for  a 
given  tube  size  if  a  powerful  igniter  is  used.  Once  initiated,  the  single-head 
spinning  detonation  does  not  appear  to  decay  at  all  as  it  propagates  down 
the  tube.  It  is  analogous  to  an  unstable  rocket  motor  resonating  at  the 
fundamental  transverse  mode  of  instability  as  it  travels  down  the  tube  at  the 
Chapman-Jouguet  velocity.  The  experiments  of  Wolanski  et  a).  (1981)  on 
methane-air  detonation  are  an  example  of  this,  since  the  5-cm-diameter 
tube  used  is  much  smaller  than  that  required  for  the  propagation  of 
methane-air  detonations.  Yet,  no  noticeable  attenuation  has  been  observed, 
and  over  a  wide  range  of  methane  concentration,  the  detonation  propa¬ 
gates  at  the  fundamental  spinning  mode.  The  existence  of  such  phenomena 
makes  it  extremely  difficult  to  establish  the  limits  experimentally.  Perhaps 
there  correspond  certain  limits  where,  irrespective  of  the  strength  of  the 
igniter,  a  supersonic  wave  at  the  fundamental  spinning  mode  propagating 
at  about  the  Chapman-Jouguet  velocity  is  not  possible.  However,  there 
exists  a  range  of  composition  between  the  first  onset  of  the  fundamental 
spinning  mode  in  a  given  tube  and  the  composition  when  it  is  impossible  to 
initiate  a  supersonic  combustion  wave.  The  question  then  arises  as  to 
whether  the  supersonic  combustion  waves  in  this  range  are  truly  self- 
sustained  detonations.  That  the  single-head  spinning  detonation  is  the 
lowest  possible  stable  mode  in  a  given  tube  (thus  the  limit)  was  proposed 
long  ago  by  Dove  &  Wagner  (1960),  among  others.  However,  the  recent 
study  by  Donato  ( 1 982)  perhaps  sheds  some  light  on  this  question  of  what  is 
considered  as  truly  self-sustained  single-headed  spin.  In  studying  the 
stability  of  the  near-limit  detonations,  Donato  found  thai  by  adding  a  finite 
perturbation  (such  as  a  few  turns  of  a  wire  spiral  inside  the  tube)  to  destroy 
the  spinning  wave,  the  detonation  will  recover  and  continue  to  propagate  as 
a  spinning  wave  only  when  the  composition  of  the  mixture  corresponds  to 
the  first  onset  of  single-headed  spinning  detonation  in  that  particular  tube. 
In  other  words,  although  the  use  of  a  stronger  igniter  can  initiate  a  single¬ 
head  spinning  detonation  wave  below  the  composition  when  spinning 
detonation  is  first  observed,  the  wave  is  not  stable  to  the  finite  perturbation 
by  the  wire  spirals.  Once  destroyed,  the  wave  will  not  transit  back  to  the 


detonation  mode.  The  more  recent  studies  by  Lee  et  al.  ( 1983)  of  turbulent- 
flame  accelerations  in  tubes  and  transition  to  detonation  confirm  the  results 
of  Donato.  It  was  found  that  transition  to  detonation  is  only  possible  in 
mixtures  bounded  by  the  composition  limits  that  correspond  to  the  first 
onset  of  single-head  spin.  In  other  words,  transition  to  detonation  is  not 
possible  in  mixtures  where  the  composition  is  poorer  than  that  correspond¬ 
ing  to  the  onset  of  single-headed  structure.  For  the  onset  of  a  single-head 
spin,  the  tube  circumference  corresponds  to  one  complete  cell  X.  Thus,  we 
have  nD  =  /,  which  for  a  given  tube  diameter  is  the  critical  value  of  X. 
Hence,  the  corresponding  mixture  composition  can  be  determined.  For 
example,  for  stoichiometric  H2-air  mixtures,  X  =  1.5  cm.  Thus,  the 
minimum  tube  diameter  that  could  sustain  a  stable  spinning  detonation 
would  be  D  ~  0.5  cm.  For  the  alkanes,  the  minimum  diameter  would  be  of 
the  order  of  D  ~  1.7  cm  for  stoichiometric  composition.  It  is  of  interest  to 
note  that  according  to  the  cell  size  reported  by  Moen  et  al.  (1983)  for 
methane  (i.e.  X  =s  33  cm),  a  minimum  tube  diameter  for  stable  detonation 
propagation  in  stoichiometric  methane-air  mixtures  would  be  of  the  order 
of  at  least  10  cm.  Thus,  the  results  reported  by  Wolanski  et  al.  (1981)  in  a  5- 
cm  tube  should  correspond  to  overdriven  transient  waves  only.  This  is 
supported  by  the  fact  that  they  observed  single-head  spinning  waves  over 
the  whole  range  of  fuel  concentrations  studied.  Based  on  stability 
considerations,  the  composition  limit  can  now  at  least  be  defined 
experimentally. 

With  the  cell  size  known  as  a  function  of  composition,  the  limits  for  any 
given  circular  tube  can  now  be  specified  by  the  criterion  X  —  nD.  For  two- 
dimensional  planar  channels,  Vasiliev  (1982)  recently  found  that  the  critical 
condition  for  a  stable  wave  corresponds  to  a  channel  width  w  equal  to  X. 
Thus,  w  =  X  represents  the  criterion  for  two-dimensional  rectangular 
channels  of  large  aspect  ratios.  What  remains  to  be  done  is  an  experimental 
study  of  the  conditions  for  stable  propagation  in  tubes  of  arbitrary 
geometries.  It  appears  that  a  unified  criterion  for  limits  in  tubes  can  be 
derived  on  the  basis  of  the  minimum  wave-curvature  concept  discussed 
earlier. 

The  question  of  the  limits  for  unconfined  cylindrical  and  spherical 
detonations  is  much  more  difficult.  Here  there  are  no  boundaries  to  be 
considered,  and  the  answer  must  lie  within  the  detailed  processes  of  the 
cellular  structure.  In  particular,  the  limits  must  be  determined  from  the 
conditions  whereby  the  cells  can  multiply  in  an  unconfined  wave.  Stable 
propagation  requires  that  the  averaged  cell  size  be  constant  for  a  given 
mixture.  Hence,  the  detonation  cells  must  multiply  in  a  diverging  wave 
front  in  accordance  with  the  rate  at  which  the  surface  area  increases.  The 
mechanism  for  the  multiplication  of  the  cells  is  in  the  transient  processes  of 
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double  Mach  reflections  (in  particular,  how  the  second  Mach  stem  is 
formed  and  develops  to  form  a  new  cell).  Even  in  nonreacting  media,  the 
problem  of  describing  various  types  of  triple-shock  interactions  has  not 
been  fully  resolved.  For  the  cellular-detonation  case,  further  complexities 
are  introduced  by  the  addition  of  flame  fronts  into  the  already  complex 
shock  pattern.  The  interaction  process  is  also  highly  transient,  and  similar 
types  of  approximations  may  not  apply.  Apart  from  the  fundamental  limit 
mechanism  of  cell  reproduction  for  unconfined  detonations,  it  is  also 
difficult  to  formulate  the  necessary  operational  definition  for  the  experi¬ 
mental  measurements  of  limits  for  diverging  waves.  It  appears  that  the 
current  method  of  arbitrarily  specifying  an  upper  initation  energy  above 
which  the  mixture  is  rendered  nondetonable  is  reasonable,  at  least  from  a 
practical  point  of  view.  As  such,  the  cell  size  also  enters  through  the 
initiation  energy.  Specifying  an  upper  value  for  the  initiation  energy  is  then 
equivalent  to  specifying  a  maximum  value  of  the  cell  size  beyond  which  the 
mixture  is  too  insensitive  to  be  detonable.  Although  there  are  other  means 
of  specifying  reasonable  limits  for  unconfined  detonation  from  the  practical 
point  of  view,  such  as  the  use  of  the  size  of  the  detonation  kernel  as  the 
minimum  detonable  volume,  the  limit  problem  remains  fundamentally 
linked  to  the  nature  of  the  cellular  structure. 

CONCLUSIONS 

The  cell  size  has  been  demonstrated  to  be  the  most  fundamental  parameter 
characterizing  the  dynamic  detonation  properties.  It  has  also  been 
demonstrated  that  a  knowledge  of  the  detonation  cell  size  permits  the 
dynamic  parameters  to  be  determined  reasonably  well.  However,  there  is 
much  room  for  refining  the  existing  relationships  that  link  A  to  these 
dynamic  parameters.  The  13A  correlation  between  the  cell  size  and  the 
critical  tube  diameter  dc  is  presently  a  purely  empirical  relationship.  Its 
validity  is  based  on  experiments  on  the  common  hydrocarbon  fuels. 
Although  the  ranges  of  initial  conditions  and  mixture  compositions  studied 
are  fairly  wide,  all  the  fuels  studied  have  similar  kinetics  of  hydrocarbon 
oxidation.  Thus  it  seems  worthwhile  to  test  the  de  =  1 3A  correlation  in 
other  detonable  systems,  such  as  H2-C12  or  CS2-02  mixtures,  where  other 
fuel  molecules  are  involved.  Also,  if  the  13A  correlation  is  related  to  the 
concept  of  an  effective  or  hydrodynamic  thickness  of  the  detonation  wave, 
then  the  use  of  various  inert  diluents  to  modify  the  sound  speed  in  the 
product  gases  should  influence  the  equilibrations  of  the  transverse  wave 
and,  hence,  the  hydrodynamic  thickness.  This  may  have  an  effect  on  the 
correlation  dc  =  13/..  Strehlow  (1968)  has  already  demonstrated  that  the 
use  of  heavy  argon  dilution  has  significant  influence  on  the  “regularity”  of 


the  cellular  structure  of  H2-02  detonations.  Further  study  of  the  transmis¬ 
sion  from  different  geometries,  as  well  as  the  transmission  from  rigid  to 
flexible  wall  tubes,  could  provide  additional  insight  into  the  physics  behind 
the  dc  =  13A  correlation.  A  correct  theoretical  model  to  explain  this  simple 
correlation  will  pave  the  way  to  the  development  of  quantitative  theories 
for  the  other  dynamic  parameters. 

Regarding  the  prediction  of  k  from  the  given  kinetics  and  physical 
properties  of  the  mixture,  it  is  the  opinion  of  the  author  that  “brute  force” 
numerical  computations  will  not  lead  to  the  solution  of  the  problem.  The 
feasibility  of  numerical  simulations  and  the  propagation  of  two- 
dimensional  cellular  detonations  have  already  been  demonstrated.  The 
extension  to  the  simulation  of  three-dimensional  cellular  detonations  is,  in 
principle,  possible.  It  is,  however,  not  clear  what  additional  physics  can  be 
learned  from  these  complex  computations.  Certainly,  the  use  of  these  codes 
to  compute  the  cell-size  data  from  the  basic  kinetics  of  the  system  is  not 
feasible  from  the  economics  standpoint.  It  would  be  much  simpler  to  obtain 
cell-size  data  experimentally.  There  is,  however,  an  urgent  need  to  improve 
both  the  technique  and  accuracy  for  cell-size  measurements  and  to 
eliminate  the  “experience  factor”  in  the  identification  of  the  dominant  cell 
size  in  an  irregular  pattern.  If  a  relationship  between  the  cell  size  and 
another  more  easily  measured  dynamic  parameter  can  be  firmly  es¬ 
tablished,  then  perhaps  the  cell  size  can  be  determined  indirectly  from  the 
measurement  of  this  parameter. 

Thus,  it  seems  that  if  k  can  be  treated  as  a  fundamental  parameter 
characterizing  the  detonation  mode  of  combustion,  then  we  are  well  on  the 
way  to  developing  predictive  theories  for  the  important  dynamic  para¬ 
meters  of  practical  interest.  The  cell  size  k,  like  the  fundamental  kinetic-rate 
constants  of  the  elementary  reactions,  can  be  considered  as  a  fundamental 
property  of  the  mixture.  Adopting  such  a  point  of  view,  the  theory  of 
gaseous  detonation  can  be  said  to  be  very  close  to  completion  after  a 
century  since  the  initial  discoveries. 
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High  Speed  Turbulent  Deflagrations  and  Transition  to  Detonation 

in  H2-Air  Mixtures 

JOHN  H.  S.  LEE,  R.  KNYSTAUTAS,  and  A.  FREIMAN* 

Department  of  Mechanical  Engineering.  McGill  University,  Montreal.  Canada 

An  experimental  investigation  on  flame  deceleration  and  transition  to  detonation  in  H  ,-air  mixtures  has  been  earned  out 
in  a  tube  which  had  a  5  cm  cross-sectional  diameter  and  was  1 1  m  long  Obstacles  in  the  form  of  a  spiral  coil  (6  mm 
diameter  tubing,  pitch  5  cm.  blockage  ratio  BR  =  0.44)  and  repealed  onfice  plates  spaced  5  cm  apart  with  blockage 
ratios  of  BR  =  0.44  and  0.6  were  used.  The  obstacle  section  was  3  m  long.  The  compositional  range  of  H2  in  air 
extended  from  I0io45%,  the  initial  pressure  of  the  experiment  was  I  atm,  and  the  mixture  was  at  room  temperature  The 
results  indicate  that  steady-state  flame  (or  detonation)  speeds  are  attained  over  a  flame  travel  ol  10-40  lube  diameters. 
ForH,  <  1 3%  maximum  flame  speeds  are  subsonic,  typically  below  200  m/s  A  sharp  transition  occurs  at  about  I3%H2 
when  the  flame  speed  reaches  supersonic  values.  A  second  transition  to  the  so-called  quasi-detonation  regime  occurs 
near  the  stoichiometric  composition  when  the  flame  speed  reaches  a  critical  value  of  the  order  of  800  m/s.  The  maximum 
value  of  the  averaged  pressure  is  found  to  be  between  the  normal  C-J  detonation  pressure  and  the  constant  volume 
explosion  value.  Of  particular  interest  is  the  observation  that  at  a  critical  composition  of  about  17%  H2  transition  to 
normal  C-J  detonation  occurs  when  the  flame  exits  into  the  smooth  obstacle-free  portion  of  the  tube.  For  compositions 
below  17%  H2.  the  high  speed  turbulent  deflagration  is  observed  to  decay  in  this  portion  of  the  lube.  The  detonation  cell 
size  for  1 7%  H  2  is  about  1 50  mm  and  corresponds  closely  to  the  value  of  trD  that  has  been  proposed  to  designate  the 
onset  of  single-head  spinning  detonation,  in  this  case  for  the  5  cm  diameter  tube  used.  This  supports  the  limit  criterion, 
namely,  that  for  confined  detonations  in  tubes,  the  onset  of  single-head  spin  gives  the  limiting  composition  for  stable 
propagation  of  a  detonation  wave. 


1  INTRODUCTION 

The  study  of  flame  propagation  through  an  array 
of  orifice  plates  arranged  in  periodic  sequence 
along  the  length  of  a  flame  tube  was  pioneered 
by  Chapman  and  Wheeler  [1)  more  than  50 
years  ago.  In  this  early  work,  the  rapid  accelera¬ 
tion  of  the  flame  as  a  result  of  the  action  of  such 
restrictive  obstructions  to  velocities  in  excess  of 
400  m  s  in  premixed  methane-air  mixtures  was 
observed.  In  the  midseventies  Wagner  [2)  extended 
the  preliminary  nature  of  Wheeler's  early  study 
to  a  mure  exnaustive  parametric  investigation  in 
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which  the  effects  of  blockage  ratio  and  orifice 
spacing  were  considered.  Wagner’s  work  indicates 
that  under  optimal  conditions,  flame  velocities 
in  excess  of  700  m/s  for  stoichiometric  CH4-air 
mixtures  can  be  achieved  just  beyond  a  short 
array  of  21  equispaced  obstacles  in  a  4  cm  diameter 
tube.  Subsequently,  research  on  turbulent  flame 
acceleration  by  repeated  obstacles  in  one-dimen¬ 
sional  channels  and  in  cylindrical  diverging  con¬ 
figurations  has  been  carried  out  at  McGill  Uni¬ 
versity  [3,  4]  and  extended  to  fairly  large  scale 
in  collaboration  with  Eckhoff  at  theChr.Michelson 
Institute  (5 J  (2£  m  diam  X  10  m  long  tube) 
and  with  Wagner  at  the  Physical  Chemistry  In¬ 
stitute  (6 J  (2^  m  X  2£  m  channel).  However, 
in  all  these  previous  studies,  the  flame  tubes 
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were  too  short  on  a  relative  scale  (L/D  small)  interest  in  hydrogen  deflagrations  has  additional 

to  permit  steady-state  conditions  to  be  achieved  impetus  in  view  of  the  current  concern  of  explo- 

at  the  end  of  the  flame  travel.  Hence,  the  maxi-  sion  hazards  in  light  water  nuclear  reactors  due 

mum  turbulent  flame  speed  for  a  given  mixture  to  the  possibility  of  large  scale  hydrogen  release 

and  obstacle  configuration  was  not  reached,  and  following  loss  of  coolant  and  core  meltdown- 

important  information  on  whether  transition  to  type  accidents,  exemplified  by  the  Three  Mile 

detonation  is  possible  could  not  be  deduced  from  Island  incident  in  1979. 

these  experiments.  From  a  fundamental  point  of 

new.  it  is  of  great  interest  to  explore  the  media-  2  EXPERIMENTAL  DETAILS 
nisms  involved  in  curbing  the  powerful  flame 

acceleration  mechanism  to  yield  eventually  an  A  steel  combustion  tube  5  cm  in  diameter  and  1 1 

optimal  turbulent  burning  rate.  It  is  also  of  great  m  long  (L/D  =  220)  was  used  in  the  present  ex- 

interest  to  examine  the  propagation  mechanism  periments.  The  tube  was  equipped  with  three  types 

of  very  high  speed  deflagrations  which  may  lead  of  turbulence  inducing  obstacles:  (i)  a  spiral  coil 

to  the  understanding  of  the  critical  conditions  (Shcheikin  spiral)  made  from  6  mm  diameter 

underlying  the  transition  process  to  the  detonation  tubing  with  a  pitch  of  5  cm  and  effective  blockage 

regime.  From  a  practical  point  of  view,  informa-  ratio  BR  =  1  -  ( d/D )2  =  0.44;  (ii)  circular  orifice 

tion  on  maximum  burning  rates  which  would  per-  plates  with  d  =  37.5  mm,  BR  =  0.44,  and  5  cm 

rnit  the  estimation  of  pressure  development  in  plate  separation;  (iii)  circular  orifice  plates  with 

explosions,  run-up  distances,  and  critical  condi-  d  =  31.6  mm,  BR  =  0.6,  and  also  5  cm  plate 

tions  for  transition  to  detonation  is  of  importance  separation  (Fig.  1).  The  length  of  each  obstacle 

in  the  context  of  probabilistic  risk  assessment.  section  was  3  m  and  the  obstacles  were  located 

The  present  paper  reports  some  recent  experi-  to  start  at  the  ignition  end.  Thus  there  was  an  8  m 

mental  results  on  flame  acceleration  in  long  section  of  smooth  obstacle-free  tube  after  the 

tubes  (L/D  >  200)  in  the  presence  of  repeated  flame  left  the  obstacle  region.  Ionization  probes 

obstacles.  Hydrogen-air  mixtures  have  been  used  spaced  at  50  cm  intervals  along  the  tube  lengths 

because  they  are  sufficiently  sensitive  so  that  run-  were  used  for  flame  speed  measurement.  The 

up  distances  to  the  steady  state  are  short  and  probe  tips  were  extended  to  the  axis  of  the  tube, 

transition  to  detonation  can  be  studied  over  a  and  thus  the  flame  speeds  reported  correspond 

wider  range  of  H2  concentrations  than  less  sensi-  to  those  along  the  tube  axis.  The  ionization 

live  hydrocarbon-air  mixtures  (e.g.,  CH4).  The  gauges  could  be  replaced  by  piezoelectric  pres- 


ORIFICE  PLATE  OBSTACLES  —  pitch  ■  0  ,  blockage  rotio  •  I  -  0-60 

D 

Hg.  1.  Schematic  diagram  of  the  flame  tube  and  obstacle  arrays. 
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sure  transducers  (PCB  1 13A24)  for  recording  the 
pressure-time  development.  The  signals  were 
recorded  using  a  Biomation  digital  waveform 
recorder  (Model  805)  to  permit  the  control  of 
signal  amplitude  and  time  scales  for  data  analysis. 
For  most  experiments,  ignition  of  the  mixture 
was  achieved  by  a  glow  plug  igniter;  in  weaker 
mixtures,  a  more  powerful  electric  spark  discharge 
was  used  to  improve  the  reliability  of  initiation. 
The  explosive  mixture  was  prepared  in  a  con¬ 
tinuous  flow  system,  with  the  desired  composi¬ 
tion  achieved  by  regulating  the  respective  flow 
rates  of  hydrogen  and  air.  Standard  calibrated 
rotameters  (Gilmont  G-2 155)  were  used  to  monitor 
the  flow  rates.  To  ensure  a  homogeneous  mix¬ 
ture,  the  gases  were  (lowed  through  a  mixing 
chamber  equipped  with  "baffle  plates  prior  to 
entering  the  combustion  tube.  The  tube  was 


first  evacuated,  and  then  no  less  than  5  tube 
volumes  of  the  explosive  mixture  were  flushed 
through  the  tube  before  each  experiment.  The 
initial  pressure  of  the  experiments  was  atmos¬ 
pheric,  and  the  temperature  was  the  ambient 
value.  The  range  of  H2  concentrations  in  the 
present  study  was  10%  £  H2  ~  45%. 

3.  RESULTS  AND  DISCUSSION 

The  variation  of  the  flame  speed  (relative  to 
laboratory  coordinates)  with  distance  along  the 
combustion  tube  for  the  case  of  the  spiral  obstacle 
(pitch  equals  diameter;  blockage  ratio  BR  =  0.44) 
is  shown  in  Fig.  2  for  various  H2  concentrations. 
Although  the  total  tube  length  was  IT  in,  with 
the  First  3  m  constituting  the  obstacle  field,  only 
the  First  8  m  were  equipped  for  measurement. 


!  '-•  -•  Variation  ot  flame  velocity  along  the  flame  tube  when  a  Shchclkin  spiral  ob¬ 
stacle  array  <  BR  =  0.44)  is  used  in  the  first  3  m  of  the  tube. 
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Within  the  obstacle  field,  for  H2  <  17%,  the 
deflagration  was  observed  to  accelerate  to  a 
steady  state  after  about  1.5-2  m  (or  30-40  tube 
diameters)  of  travel.  However,  for  higher  concen¬ 
trations  of  H2.  the  steady  state  was  achieved  in 
less  than  1  m  for  the  present  case  of  the  spiral 
obstacle.  It  is  of  particular  interest  to  note  that 
there  appears  to  be  a  sudden  increase  in  the 
steady  velocity  from  about  200/ms  for  12% 
H2  to  almost  500  m/s  at  13%  H2.  This  sharp 
increase  in  the  flame  speed  occuring  at  around 
13%  H2  is  common  to  all  three  obstacle  con¬ 
figurations  investigated.  An  explanation  that 
suggests  itself  that  can  possibly  account  for  this 
rather  sudden  transition  is  that  there  is  a  change¬ 
over  in  the  dominant  kinetic  step  of  the  chemical 
reaction.  As  pointed  out  by  Atkinson  et  al.  [7], 
among  others  [8.  9] .  the  reaction  H  +  02  -»  OH  + 
O  is  in  competition  with  the  three-body  reaction 
H  +  02  +  M  -*  H02  +  M.  with  the  latter  dominat¬ 
ing  at  temperatures  below  about  1350°K.  At 
about  13%  H2,  the  flame  temperature  is  of  the 
order  of  1300°K.  Thus,  a  changeover  to  the  more 
rapid  OH  branching  reaction  may  take  place  and 
result  in  a  rapid  increase  in  the  maximum  turbulent 
flame  speed  observed  at  13%  H2.  A  second  sharp 
transition  is  also  noted  at  about  25%  H2,  when 
the  flame  speed  jumps  from  about  800  to  1500 
rms.  In  fact,  25%  H2  seems  to  he  the  critical 
composition  in  which  the  flame  speed  fluctuates 
between  two  regimes  randomly  from  experiment 
to  experiment.  The  high  velocity  regime  cor¬ 
responds  to  the  regular  detonation  mode  of  com¬ 
bustion.  but  the  pressure  losses  suffered  by  the 
detonation  in  propagating  through  the  obstacles 
account  for  the  large  velocity  deficit  observed 
(i.e..  1500  m/s,  as  compared  with  the  normal 
C-J  value  of  1861  m/s  for  25%  H2).  These  low 
velocity  detonations  resulting  from  severe  pres¬ 
sure  losses  inflicted  by  the  obstacles  are  hence¬ 
forth  referred  to  as  quasi-detonations.  Their 
propagation  mechanism  is  strongly  influenced 
by  shock  reflections  [IOj  and  violent  turbulent 
mixing  caused  by  the  obstacles.  This  differs  from 
the  usual  transverse  wave  collision  mechanism 
responsible  for  the  propagation  of  the  normal 
cellular  detonations  with  velocities  quite  close  to 
the  theoretical  C-J  values. 


The  sudden  transition  to  the  quasi-detonation 
regime  suggests  that  some  critical  value  for  the 
flame  speed  needs  to  be  reached  before  the  con¬ 
ditions  for  the  onset  of  detonation  can  occur. 
For  the  case  of  the  spiral  obstacle,  the  steady 
flame  speed  increases  with  H2  concentration, 
and  at  25%  H2  a  value  of  about  800  m/s  is  at¬ 
tained.  Such  a  flame  speed  drives  a  shock  wave 
of  about  4/s  =2.7,  which  when  reflected  from  the 
obstacles  yields  temperatures  of  the  order  of 
1200°K.  This  corresponds  closely  to  the  auto¬ 
ignition  temperatures  for  shock  heated  mixtures. 
These  results  suggest  that  the  necessary  condi¬ 
tion  for  the  sudden  transition  is  associated  with  a 
requirement  for  some  minimum  critical  value  of 
the  flame  speed.  This  is  compatible  with  previous 
studies  of  transition  from  deflagration  to  detona¬ 
tion  in  tubes  which  indicate  certain  critical  thresh¬ 
old  values  for  the  flame  speed  prior  to  the  onset 
of  detonation.  The  various  modes  of  transition 
to  detonation  revealed  in  the  classical  paper  of 
Urtiew  and  Oppenheint  [11]  all  indicate  the 
necessity  for  a  certain  minimum  flame  speed 
for  the  onset  of  detonation. 

To  verify  the  above  suggestion  that  a  minimum 
flame  speed  is  required  for  transition,  experiments 
were  carried  out  by  propagating  a  normal  C-J 
detonation  initiated  in  the  obstacle-free  section 
back  into  the  obstacle  region  (designated  as  back¬ 
ward  propagation).  A  hysteresis  effect  is  noted 
(e.g.,  at  20%  H2  shown  in  Fig.  3)  in  that  when  a 
fully  established  C-J  detonation  is  propagated 
back  from  the  obstacle-free  into  the  obstacle- 
filled  section,  it  settles  down  to  a  velocity  of  about 
1200  m/s,  as  compared  with  the  maximum  steady- 
state  turbulent  flame  speed  of  about  800  m/s 
that  is  attained  via  the  continuous  flame  accelera¬ 
tion  process  (designated  as  forward  propagation). 
Under  the  conditions  existing  in  the  present  ex¬ 
periment,  transition  to  detonation  via  continuous 
flame  acceleration  does  not  occur,  possibly  due 
to  the  fact  that  the  critical  value  of  flame  speed 
required  for  transition  is  not  attained.  However, 
when  started  out  as  a  fully  developed  C-J  detona¬ 
tion  prior  to  entering  the  obstacle  region,  the 
velocity  is  seen  to  decay  to  the  lower  quasi¬ 
detonation  value  corresponding  to  the  obstacle 
configuration  and  mixture  composition.  For  mix- 
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Fig.  3.  Variation  of  pressure  along  the  flame  tube  and  within  the  Shchelkin  spiral  ob¬ 
stacle  array  (BR  =  0.44). 

ture  compositions  where  transition  from  flame  speed  decays  rapidly  downstream  of  the  obstacle 

to  quasi-detonation  can  occur  in  the  obstacle  section.  This  is  due  to  the  absence  of  obstacles 

region  (i.e..  the  critical  value  of  the  flame  speed  to  randomize  the  mean  flow  to  produce  the  level 

is  attained),  no  hysteresis  effect  is  observed.  of  turbulence  intensity  required  to  sustain  the  high 

The  quasi-detonation  velocity  obtained  is  identi-  burning  rate.  Below  about  13%  H2,  the  flame soeed 

cal  for  both  cases,  whether  the  starting  point  is  decays  drastically  to  a  few  tens  of  meters*  per 

a  normal  C-J  detonation  or  a  flame  undergoing  second  within  about  20-40  tube  diameters.  For 

transition  to  a  quasi-detonation  later.  This  occurs  H2  >  13%,  the  flame  speed  first  suffers  a  sharp 

for  mixtures  ,n  the  range  25%  <  H2  <  45%  for  the  initial  drop,  then  decays  slowly  to  about  a  couple 

(Shchelkin  i  spiral  obstacle  as  shown  in  Fig,  3.  of  hundred  meiei.  per  second.  At  compositions 

Observations  downstream  of  the  obstacle  sec-  for  H2  >  17%  therr  clearly  is  a  dramatic  change 

tion  when  the  flame  or  quasi-detonation  propa-  in  the  nature  of  the  flame  decay  phenomenon, 

gates  into  the  obstacle-free  portion  of  the  tube  After  the  sharp  initial  drop  the  flame  is  seen  to 

also  yield  some  interesting  results.  From  Fig.  2  reaccelerate  and  undergo  transition  to  a  normal 

it  can  be  noted  that  below  17%  H2,  the  flame  C-J  detonation  after  a  few  meters  of  travel.  The 
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data  indicate  that  there  is  considerable  scatter  of 
results  in  this  region  typically  symptomatic  of 
the  abruptly  sporadic  nature  of  the  transition 
process.  Below  1 79?  H2,  the  flame  simply  decays 
monotonieally  to  low  velocity  with  increasing 
distance.  Thus,  17%  H2  appears  to  be  the  critical 
concentration  that  gives  rise  to  transition  to  nor¬ 
mal  detonation  in  the  obstacle-free  portion  of  the 
5  cm  diameter  tube.  However,  inside  the  obstacle 
filled  section  with  the  particular  spiral  obstacle 
configuration  used,  25%  H2  appears  to  be  the 
critical  concentration  for  which  transition  to  a 
quasi-detonation  wave  is  seen  to  occur  within 
the  obstacle  section.  In  examining  the  detonation 
cell  size  for  H2-air  mixtures  [12] ,  it  is  found  that 
for  17%  H2.  '-he  cell  size  is  about  15.0  cm.  This 
corresponds  very  closely  to  the  value  of  itD  for 
the  detonation  cell  size  for  mixtures  at  the  onset 
of  the  single-head  spinning  mode  in  a  circular 
tube  of  diameter  D.  This  confirms  the  criterion 
proposed  for  detonation  limits  in  circular  tubes 
that  the  critical  composition  should  be  based  on 
the  onset  of  single-head  spinning  detonations 
[13] -[15].  This  limit  criterion  was  proposed 
on  the  hasis  of  experimental  observations  [15] 
that  only  mixtures  whose  cell  sizes  are  less  than 
that  corresponding  to  the  single-head  spin  in  the 
given  tube  permit  the  propagation  of  a  stable 
detonation.  A  stable  detonation  in  this  case  means 
one  which,  having  been  destroyed  by  a  finite 
perturbation,  such  as  that  induced  by  blockage 
in  the  tube,  say,  via  an  orifice  plate,  will  rapidly 
undergo  transition  back  to  a  normal  detonation. 
In  the  context  of  the  above  definition  of  detona- 
bility  limits,  it  is  of  interest  to  refer  to  the  work 
of  Wolanski  et  al.  [16]  on  methane-air  detona¬ 
tions  in  a  5  cm  tube  which  yielded  single-head 
spin  waves  for  all  CH4  concentrations.  Clearly, 
their  experimental  conditions  are  beyond  the 
composition  limit  for  the  onset  of  single-head 
spin,  and  the  wave  is  not  stable  to  finite  perturba¬ 
tions.  Once  the  detonation  wave  is  destroyed, 
it  will  never  recover  to  propagate  in  the  detonation 
mode  again.  Such  waves  are  formed  in  the  first 
place  due  to  the  use  of  very  powerful  igniters. 
The  present  result  is  therefore  significant  in  that  it 
demonstrates  that  if  the  mixture  composition  is 
within  the  detonability  limits  for  the  given  tube, 


then  it  is  possible  for  transition  from  deflagration 
to  detonation  to  take  place.  It  is  evident  that  when 
and  if  transition  occurs  depends  on  the  effective¬ 
ness  of  the  flame  acceleration  mechanisms  availa¬ 
ble  to  bring  the  flame  speed  up  to  meet  certain 
critical  requirements. 

The  coi.  .sponding  pressure  variation  with 
distance  from  the  igniter  for  the  Shchelkin  spiral 
obstacle  array  is  shown  in  Fig.  4  together  with 
two  typical  pressure  traces:  one  corresponding 
to  a  pressure  measurement  close  to  the  igniter 
at  the  very  beginning  of  the  obstacle  array,  the 
other  roughly  two-thirds  of  the  way  along  the 
array.  One  can  immediately  note  that  the  pres¬ 
sure  front  steepens  with  increasing  distance 
along  the  obstacles  as  shown  on  the  oscilloscope 
traces;  this  is  qualitatively  similar  to  the  trend 
observed  by  Wagner  [2]  for  methane-air  flames. 
However,  in  the  present  case  of  the  hydrogen- 
air  flames,  the  pressure-time  records  exhibit 
strong  oscillations  over  and  above  the  baseline 
mean  pressure  changes.  This  is  due  to  the  inter¬ 
action  of  the  flow  with  the  periodically  spaced 
obstacles  and  is  compatible  with  the  wave  pattern 
shown  in  the  streak  records  reported  by  Wagner 
[2]  for  a  flame  propagating  in  a  regularly  spaced 
obstacle  array.  The  pressure  plotted  in  Fig.  4 
refers  to  the  maximum  mean  value.  The  time  to 
reach  this  maximum  mean  value  of  the  average 
pressure  ranges  from  10  ms  or  more  for  very  lean 
mixtures  (e.g„  H2  s  12%)  to  a  few  microseconds 
in  the  quasi-detonation  regime  (H2  ~  25%). 
For  very  lean  mixtures,  the  maximum  pressure 
obtained  is  the  same  for  all  the  transducer  loca¬ 
tions.  This  is  due  to  the  fact  that  the  flame  speed 
is  less  than  the  sound  speed  (~350  m/s),  which  is 
a  measure  of  the  rate  at  which  the  local  pres¬ 
sure  inhomogeneities  generated  by  the  flame  are 
dispersed  and  equalized.  For  H2  ~  13%,  where 
the  flame  speed  becomes  supersonic,  local  pres¬ 
sure  buildup  is  dispersed  at  speeds  of  the  order 
of  the  flame  speed.  Thus,  the  pressure  front 
closely  coincides  with  the  flame  front  and  the 
approach  to  steady  state  for  both  the  flame 
speed  and  the  maximum  pressure  become  identi¬ 
cal,  as  indicated  in  Fig.  4.  Beyond  the  obstacle 
region,  the  pressure  will  either  rise  further  or 
fall,  depending  on  whether  the  flame  or  quasi- 


detonation  undergoes  transition  to  a  normal 
C-J  detonation  or  decays  continuously  to  a  rela¬ 
tively  low  speed  deflagration  compatible  with  the 
obstacle-free  portion  of  the  tube. 

A  comparison  of  the  flame  developments  for 
the  three  obstacle  configurations  of  spiral  (BR  = 
0.44),  orifice  plate  (BR  =  0.44),  and  orifice  plate 
(BR  =  0.6)  can  be  seen  from  Figs.  2,  5,  and  6. 
In  terms  of  flame  acceleration  rates,  orifice  plates 
are  more  efficient  than  the  spiral,  for  the  same 
blockage  ratio,  even  though  this  may  not  be  a 
meaningful  comparison,  since  the  spiral  has  a 
different  geometry.  However,  for  the  same  orifice 
plate  configuration,  a  higher  blockage  ratio  gives 
rise  to  steady-state  conditions  much  earlier.  The 
steady-state  flame  speeds  and/or  quasi-detonation 
velocities  are  weakly  dependent  on  the  obstacle 


configuration  and  blockage  ratio.  However,  the 
first  sudden  jump  in  the  flame  speed  occuring 
at  about  13%  H2  is  common  to  all  three  ob¬ 
stacle  configurations,  indicating  that  the  kinetic 
explanation  is  a  plausible  one.  Furthermore,  the 
reacceleration  of  the  flame  to  detonation  after 
leaving  the  obstacle  section  at  the  critical  composi¬ 
tions  above  about  17%  H2  for  the  three  cases 
also  confirms  the  single-head  spin  limit  criterion 
discussed  earlier.  A  more  convenient  compari¬ 
son  of  the  steady-state  flame  speeds  for  various 
H2  concentrations  for  the  three  obstacle  con¬ 
figurations  is  given  in  Fig.  7.  Here,  the  first  velocity 
jump  in  the  flame  speed  at  about  13%  H2  for  all 
the  three  cases  as  well  as  subsequent  dependence 
of  the  flame  speed  on  blockage  ratio  for  different 
H2  concentrations  can  be  clearly  observed.  The 
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Fig.  7.  Variation  of  the  maximum  steady-state  flame  velocity  with  composition  of  the 
H2-aii  mixture  for  the  three  obstacle  array  configurations  studied. 


second  jump,  corresponding  to  the  transition  to 
the  quasi-detonation  regime,  is  also  common  to 
all  cases.  It  is  of  particular  interest  to  note  that 
the  flame  speed  just  prior  to  the  transition  is 
about  860  m/s  for  all  three  cases.  This  confirms 
the  notion  that  perhaps  a  critical  value  of  the 
deflagration  speed  is  required  before  transition 
from  deflagration  to  detonation  can  take  place. 
For  different  obstacle  configurations  and  blockage 
ratios,  this  critical  value  is  attained  at  different 
H2  concentrations.  For  example,  with  the  large 
blockage  ratio  of  0.6,  the  pressure  losses  are  suf¬ 
ficiently  severe  to  limit  the  steady-state  flame 
speed  to  a  value  below  the  critical  value  until 
almost  the  stoichiometric  composition  is  reached. 
However,  for  the  spiral  configuration,  where  the 
pressure  loss  are  much  less  severe,  the  flame 
accelerates  more  readily  and  the  critical  value 
required  for  transition  can  be  reached  at  a  lower 
H2  concentration  of  about  25%  H2.  It  is  possible 
that  for  very  dense  obstacle  arrays  with  sufficiently 


high  pressure  losses,  transition  to  detonation  may 
be  suppressed  entirely  while  propagating  in  the 
obstacle  field,  since  the  maximum  flame  speed 
attained  may  be  below  the  required  level.  Pre¬ 
sumably.  '-c  critical  value  of  the  flame  speed 
must  also  depend  on  the  H2  concentration  as  well 
as  other  parameters.  In  the  present  study,  about 
800  m/s  appears  to  be  the  minimum  critical 
value  of  flame  speed  required  for  transition  to 
detonation.  For  H2-air  mixtures  such  a  value 
of  flame  speed  corresponds  roughly  to  H2  con¬ 
centrations  around  the  stoichiometric  composi¬ 
tion. 

For  the  higher  blockage  ratio  of  0.6,  transition 
to  detonation  again  is  suppressed  for  rich  mixtures 
with  H2  ~  33%.  Thus,  for  this  particular  obstacle 
configuration,  transition  is  observed  only  in  the 
narrow  range  of  30%  ~  H2  ~  33%.  For  lower 
blockage  ratios  (i.e.,  BR  =  0.44),  transition  is 
observed  for  rich  mixtures  up  to  45%  H2.  The 
implication  here  is  that  eventually,  for  even  richer 
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mixtures  (i.e.,  H2  ~  45%),  transition  again  would 
not  be  possible  for  the  lower  blockage  ratio 
configurations. 

It  is  of  interest  to  note  that  in  the  range  25%  5 
H2  ~  45%  for  the  case  of  BR  =  0.44,  the  quasi¬ 
detonation  velocity  remains  more  or  less  constant 
with  composition,  even  though  the  normal  C-J 
detonation  velocity  shows  a  continuous  rise  with 
increasing  H2  concentration.  This  indicates  the 
strong  influence  of  pressures  losses  generated  by 
the  obstacles  on  the  propagation  of  these  quasi¬ 
detonations.  Increasing  the  detonation  velocity 
with  H2  concentration  increases  the  momentum 
losses,  thus  maintaining  a  more  or  less  constant 
value  of  the  quasi-detonation  velocity  itself. 


However,  for  different  obstacle  configurations 
and  blockage  ratios,  the  value  of  the  quasi-detona¬ 
tion  velocity  varies.  For  example,  for  the  spiral 
the  quasi-detonation  velocity  is  around  1800  m/s 
in  the  range  25%  <  H2  <  45%,  while  for  the 
orifice  plates  with  the  same  blockage  ratio  of 
0.44,  the  quasi-detonation  velocity  is  about  1600 
m/s.  This  indicates  that  the  pressure  losses  asso¬ 
ciated  with  the  orifice  plates  are  higher.  For  the 
orifice  plates  with  BR  =  0.6,  the  quasi-detonation 
velocity  is  only  about  1300  m/s. 

The  maximum  mean  pressure  measured  at  197 
cm  from  the  ignition  end  is  plotted  in  Fig.  8 
against  H2  concentration  for  all  three  obstacle 
configurations.  From  previous  figures,  we  note 
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that  steady -state  conditions  (either  turbulent 
deflagrations  or  quasi-detonations)  are  achieved 
at  this  location.  In  Fig.  8,  we  again  note  the  rapid 
increase  in  pressure  for  all  obstacle  configurations 
at  about  13%  H2.  Prior  to  transition  to  the  quasi¬ 
detonation  regime,  the  maximum  mean  pressure 
is  found  to  be  slightly  above  the  theoretical 
constant  volume  explosion  pressure.  Dependence 
of  the  pressure  on  obstacle  configuration  is  weak 
throughout.  In  the  range  of  compositions  where 
transition  from  high  speed  deflagration  to  quasi¬ 
detonation  (i.e.,  25%  ~  H2  5  30%)  occurs,  the 
pressure  rises  only  slightly  to  values  about  midway 
between  the  equilibrium  C-J  detonation  pressure 
and  the  constant  volume  explosion  pressure. 
Unlike  the  flame  speed,  which  shows  a  more 
distinct  jump  when  transition  to  quasi-detonation 
occurs,  the  pressure  shows  only  a  slight  increase. 
In  the  quasi-detonation  regime,  the  spiral  and 
orifice  plate  obstacles  with  the  same  blockage 
ratio  give  almost  identical  pressures  even  though 
the  detonation  velocity  shows  a  distinct  difference 
of  about  200  m/s.  For  the  higher  blockage  ratio 
plates  of  BR  =  0.6,  the  quasi-detonation  pressure 
is  also  practically  identical  to  the  other  configura¬ 
tions  (H2  ~  34%).  For  richer  mixtures  (i.e., 
H2  ~  34%),  where  Fig.  7  indicates  no  transition 
to  quasi-detonation  for  BR  =  0.6,  the  pressure  for 
the  high  speed  deflagration  regime  decreases 
only  slightly  from  the  corresponding  quasi-detona¬ 
tion  values.  In  fact,  it  may  be  concluded  that  as 
far  as  the  pressure  is  concerned,  little  difference 
is  observed  between  the  high  speed  deflagration 
and  the  quasi-detonation  regimes  when  the  block¬ 
age  ratio  or  pressure  losses  are  sufficiently  high. 
The  distinction  becomes  more  prominent  when 
the  pressure  losses  are  relatively  small  (e.g.,  the 
spiral)  and  the  two  regimes  can  be  separated  more 
clearly  when  transition  occurs. 

4.  CONCLUSIONS 

The  results  of  the  present  investigation  indicate 
that  both  chemistry  and  gas  dynamics  play  impor¬ 
tant  roles  in  the  control  of  freely  propagating 
turbulent  deflagrations.  For  mixtures  near  the 
limits  where  the  reaction  rates  are  low,  it  appears 
that  turbulent  quenching  is  responsible  for  the 


upper  bound  of  the  burning  rate  and  hence  of  the 
deflagration  speed.  In  other  words,  the  self-ac¬ 
celeration  process  will  bring  the  flame  up  to 
velocities  where  the  turbulent  mixing  rate  be¬ 
comes  of  the  same  order  of  magnitude  as  the 
reaction  rate.  Further  increase  in  turbulence  in¬ 
tensity  will  result  in  quenching  rather  than  an 
augmentation  of  the  combustion  processes. 
Equality  between  some  characteristic  eddy  life 
time  and  the  global  chemical  reaction  time  would 
be  an  appropriate  criterion  for  this  regime. 

For  sufficiently  short  reaction  times,  quenching 
will  not  set  in  until  the  displacement  flow  generated 
by  the  flame  starts  to  choke  gas  dynamically  from 
the  pressure  losses  due  to  the  obstacles.  Thus, 
frictional  choking  now  poses  the  upper  limit 
rather  than  the  chemical  reaction  rate,  and  the 
maximum  deflagration  speed  will  depend  on  the 
losses  as  well  as  the  power  output  of  the  combus¬ 
tion  processes  that  is  responsible  for  the  genera¬ 
tion  of  the  displacement  flow  itself.  Choking  or 
sonic  flow  conditions  present  the  criterion  for 
this  gas  dynamic  regime.  For  the  present  case  of 
H2-air  mixtures,  the  sudden  increase  in  reaction 
rate  at  around  13%  H2  leads  to  a  distinct  transi¬ 
tion  from  the  turbulent  quenching  to  the  gas 
dynamic  choking  regimes.  In  other  mixtures 
where  such  a  distinct  change  in  the  chemical 
reaction  rates  with  composition  is  absent,  it  would 
be  difficult  to  discern  the  two  regimes.  The  onset 
of  sonic  flow  conditions  would  perhaps  signify 
that  choking  is  present  and  that  the  gas  dynamic 
mechanism  begins  to  play  a  more  prominent  role. 

A  second  distinct  regime,  the  so-called  quasi¬ 
detonation  regime,  is  also  observed.  The  onset 
of  the  quasi-detonation  regime  is  marked  by  a 
sudden  jump  in  the  velocity  of  the  turbulent 
deflagration,  although  the  averaged  maximum 
pressure  does  not  demonstrate  such  a  distinct 
difference  between  the  two  regimes.  Depending 
on  the  pressure  loss,  the  quasi-detonation  velocities 
can  be  as  low  as  30%  of  the  normal  C-J  values, 
but  the  maximum  averaged  pressure  for  high 
speed  deflagrations  and  quasi-detonations  is 
about  the  same,  having  values  typically  midway 
between  the  constant  volume  explosion  pressure 
and  the  normal  C-J  detonation  pressure.  The 
pressure-time  histories  for  high  speed  turbulent 
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deflagrations  and  quasi-detonations  are  quite 
similar,  and  only  the  velocity  jump  marks  the 
transition  between  the  two  regimes.  Transition 
to  quasi-detonations  appears  to  require  a  minimum 
flame  speed  of  the  order  of  800  m/s.  In  other 
words,  for  a  given  obstacle  configuration,  the 
mixture  composition  in  which  transition  pccurs 
appears  to  be  governed  by  the  attainment  of  this 
minimum  flame  speed  of  about  800  m/s. 

Regarding  acceleration  to  steady  state,  it  is 
found  that  higher  pressure  losses  (i.e.,  larger 
blockage  ratio)  give  more  rapid  rate  of  flame 
acceleration.  This  is  compatible  with  the  notion 
that  pressure  losses  are  associated  with  the  ran¬ 
domization  of  the  mean  flow  into  turbulence, 
which  in  turn  augments  the  combustion  rates. 
For  very  lean  mixtures  (e.g.,  H2  as  10%),  the 
steady  state  is  achieved  very  rapidly,  suggesting 
that  turbulent  quenching  is  the  dominant  mech¬ 
anism  and  sets  in  immediately,  since  the  slow 
reaction  rate  does  not  require  high  turbulent 
intensities.  For  more  rapidly  burning  mixtures, 
the  flame  acceleration  to  steady  state  takes  cor¬ 
respondingly  longer,  since  it  takes  time  to  build 
up  a  sufficiently  high  gas  dynamic  flow  velocity 
to  produce  the  required  turbulence  intensity  for 
the  quenching  mechanism  to  be  effective.  For 
sufficiently  high  velocity,  frictional  choking 
may  set  in  and  present  another  mechanism  that 
limits  the  rate  of  buildup  of  the  mean  flow  velocity 
itself.  For  H2-air  mixtures,  the  very  rapid  increase 
in  the  kinetic  reaction  rates  at  around  13%  via 
the  OH  branching  step  results  in  a  very  rapid 
increase  in  the  turbulent  flame  speed  to  sonic 
velocity,  and  choking  then  becomes  the  dominant 
mechanism.  In  other  mixtures  with  a  different 
kinetic  mechanism,  there  may  not  be  a  sudden 
transition,  and  no  sharp  jump  in  velocity  would 
be  observed. 

Of  particular  interest  are  the  phenomena  down¬ 
stream  of  the  obstacle  section.  There  appears  to 
be  a  critical  H2  concentration  of  about  17%  H2 
above  which  the  deflagration  upon  leaving  the 
obstacle  filled  section  will  abruptly  undergo 
transition  to  a  normal  C-J  detonation.  Below 
this  critical  concentration,  the  deflagration  decays, 
since  there  are  no  more  obstacles  to  support  the 


required  turbulence  level  for  high  speed  combus¬ 
tion.  Corresponding  to  the  critical  composition 
of  17%  H2,  the  detonation  cell  size  is  about  156 
mm  and  agrees  closely  with  the  value  of  vD  for 
the  onset  of  single-head  spinning  detonation  in 
the  5  cm  tube  used.  This  is  in  accord  with  the 
criterion  that  the  onset  of  single-headed  spin 
detonations  represents  the  limiting  condition 
for  stable  detonation  propagation  in  a  given 
tube.  The  present  results  go  further  in  demonstrat¬ 
ing  that  transition  is  only  possible  in  a  given  tube 
if  the  mixture  composition  is  such  that  the  cell 
size  is  less  than  the  itD  value  for  single-head  spin. 
This  suggests  that  establishing  detonability  limits 
could  be  approached  by  carrying  out  experiments 
exploring  the  possibility  for  transition  itself. 

The  present  results,  which  demonstrate  the 
high  deflagration  speed  achieved  when  H2  ~  13%, 
suggest  that  the  lean  limit  for  H2-air  detonation 
could  be  in  the  neighborhood  of  13%  H2,  as 
suggested  by  Atkinson  et  al.  [7] .  In  fact,  our 
recent  experiments  [17]  in  a  30  cm  diam  tube 
demonstrate  that  13.8% H2  can  be  detonated  with 
a  corresponding  cell  size  of  about  60  cm.  For  13% 
H2 ,  with  an  estimated  cell  size  of  the  order  of  80 
cm,  the  minimum  tube  diameter  for  transition  or 
stable  single-head  detonation  would  be  of  the 
order  of  about  27  cm. 
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